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Abstract 



In this thesis I present recent work aimed at showing how currently competing 
theories of the early universe leave their imprint on the temperature anisotropics of 
the cosmic microwave background (CMB) radiation. 

After some preliminaries, where we review the current status of the field, we con- 
sider the three-point correlation function of the temperature anisotropics, as well as 
the inherent theoretical uncertainties associated with it, for which we derive explicit 
analytic formulae. 

These tools are of general validity and we apply them in the study of possible 
non-Gaussian features that may arise on large angular scales in the framework of both 
inflationary and topological defects models. 

In the case where we consider possible deviations of the CMB from Gaussian statis- 
tics within inflation, we develop a perturbative analysis for the study of spatial corre- 
lations in the infiaton field in the context of the stochastic approach to inflation. 

We also include an analysis of a particular geometry of the CMB three-point func- 
tion (the so-called 'collapsed' three-point function) in the case of post-recombination 
integrated effects, which arises generically whenever the mildly non-linear growth of 
perturbations is taken into account. 

We also devote a part of the thesis to the study of recently proposed analytic models 
for topological defects, and implement them in the analysis of both the CMB excess 
kurtosis (in the case of cosmic strings) and the CMB collapsed three-point function 
and skewness (in the case of textures) . 

Lastly, we present a study of the CMB anisotropics on the degree angular scale in 
the framework of the global texture scenario, and show the relevant features that arise; 
among these, the Doppler peaks. 
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Units and Conventions 



In this thesis we are concerned with the study of specific predictions for the cosmic 
microwave background (CMB) anisotropics as given by different early universe models. 
Thus, discussions will span a wide variety of length scales, ranging from the Planck 
length, passing through the grand unification scale (e.g., when infiation is supposed 
to have occurred and cosmological phase transitions might have led to the existence 
of topological defects), up to the size of the universe as a whole (as is the case when 
probing horizon-size perturbation scales, e.g., from the COBE-DMR maps).[] 

As in any branch of physics, it so happens that the chosen units and conventions 
are almost always dictated by the problem at hand (whether we are considering the 
microphysics of a cosmic string or the characteristic thickness of the wake of accreted 
matter that the string leaves behind due to its motion). It is not hard to imagine then 
that also here disparate units come into the game. Astronomers will not always feel 
comfortable with those units employed by particle physicists, and vice versa, and so it 
is worthwhile to spend some words in order to fix notation. 

The so-called natural units, namely h = c = ks = 1, will be employed, unless 
otherwise indicated, h is Planck's constant divided by 2tt, c is the speed of light, and 
ks is Boltzmann's constant. Thus, all dimensions can be expressed in terms of one 
energy-unit which is usually chosen as GeV = lO^eV, and so 

[Energy] = [Mass] = [Temperature] = [Length]"^ = [Time]^^ . 

Some conversion factors that will be useful in what follows are 



1 GeV 


= 1.60x10 ^ erg 


(Energy) 


1 GeV 


= 1.16x10^3 K 


(Temperature 


1 GeV 


= 1.78x10-24 g 


(Mass) 


1 GeV-i 


= 1.97x10"^^ cm 


(Length) 


1 GeV-i 


= 6.58x10-25 sec 


(Time) 



mp stands for the Planck mass, and associated quantities are given (both in natural 

^COBE stands for COsmic Background Explorer satellite, and DMR is the Differential Microwave 
Radiometer on board. 



VI 1 



and cgs units) by nip = 1.22 x lO^^GeV = 2.18 x IQ-^g, Ip = 8.2 x IQ-^OGeV-^ 
= 1.62 X 10~^^cm, and tp — 5.39 x 10~^sec. In these units Newton's constant is given 
by G = 6.67 x IQ-^cm^^g-^sec-^ = mp^. 

These fundamental units will be replaced by other, more suitable, ones when study- 
ing issues on large-scale structure formation. In that case it is better to employ astro- 
nomical units, which for historical reasons are based on solar system quantities. Thus, 
the parsec (the distance at which the Earth-Sun mean distance (1 a.u.) subtends one 
second of arc) is given by 1 pc = 3.26 light years. Of course, even this is much too 
small for describing the typical distances an astronomer has to deal with and thus 
many times distances are quoted in megaparsecs, 1 Mpc = 3.1 xlO^^cm. Regarding 
masses, the standard unit is given by the solar mass, Mq = 1.99 x lO^^g. 

It is worthwhile also recalling some cosmological parameters that will be used ex- 
tensively in this thesis, like for example the present Hubble time Hq^ — 3.09 x 10^'^h~^ 
sec = 9.78/i~^ Gyr, where h {little H) parameterises the uncertainty in the present 
value of the Hubble parameter Hq — 100 h km sec~^Mpc~^. Agreement with the pre- 
dictions of the hot big bang model plus direct observations yield the (conservative) 
constraint 0.4 < h < I. The Hubble distance is cH^'^ = 2997.9/i~^ Mpc, and the 
critical density pc = 2>HI/{S'kG) — 1.88/i^ x 10~^^ g cm"^. Lastly, the mean value (the 
monopole) of the CMB radiation temperature today is 2.726 K. 

Throughout we adopt the following conventions: Greek letters denote spacetime 
indices and run from to 3; spatial indices run from 1 to 3 and are given by Latin 
letters from the middle of the alphabet, while those from the beginning of the alphabet 
stand for group indices (unless stated otherwise). The metric signature is taken to be 

(+,-,-,-)• 

Many more conversion factors, as well as the values of fundamental constants and 
physical parameters (both astrophysical and cosmological) can be found in the appen- 
dices of the monograph by Kolb & Turner [1990]. 

Some abbreviations used in this thesis are the following. 

CMB = cosmic microwave background CDM = cold dark matter 
FRW = Fricdmann-Robcrtson-Walker HDM = hot dark matter 
GUT = grand unified theory 

and, lastly, COBE and DMR that we mentioned before. 
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Introduction 



In 1964 Penzias and Wilson discovered the Cosmic Microwave Background (CMB) ra- 
diation. Since then researchers have been studying the radiation background by all 
possible means, using ground, balloon, rocket, and satellite based experiments [Read- 
head Sz Lawrence, 1992]. This relic radiation is a remnant of an early hot phase of the 
universe, emitted when ionised hydrogen and electrons combined at a temperature of 
about 3000 K. Within the simplest models this occurs at a redshift z 1100 — 1300, 
although it is also possible that the hydrogen was reionised as recently as redshift 
z 100 [Jones & Wyse, 1985]. The CMB radiation is a picture of our universe when 
it was much smaller and hotter than it is today. 

It is well known that the CMB has a thermal (blackbody) spectrum, and is re- 
markably isotropic and uniform; in fact, after spurious effects arc removed from the 
maps, the level of anisotropies is smaller than a part in IC^. Only recently have 
experiments reliably detected such perturbations away from perfect isotropy. These 
perturbations were expected: a couple of years after the discovery of the CMB, Sachs 
& Wolfe [1967] showed how variations in the density of the cosmological fluid and 
gravitational wave perturbations result in CMB temperature fluctuations, even if the 
surface of last scattering was perfectly uniform in temperature. In the following years 
several other mechanisms for the generation of anisotropies, ranging all angular scales, 
were unveiled. 

Our aim in this thesis is to report a study of the structure of the CMB anisotropies 
and of the possible early universe models behind it. In this Introduction we briefly 
review the research done and we cite the references where the original parts of the 
thesis are published. More detailed accounts of the different topics treated here are 
given in the various Chapters below. 

In the first Chapter we include a brief account of the basics of what is considered 
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to be the 'standard cosmology', the why's and how's of its success as well as of its 
pressing problems. 

The second (rather lengthy) Chapter is devoted to an overview of the salient fea- 
tures of the two presently favourite theories of the early universe (or should we call 
them (just) 'ideas' or (pompously) 'scenarios'?), which were devised to complement 
the standard cosmology and render it more satisfactory, both from the theoretical 
and observational viewpoint. Inflationary models and 'seed' models with topological 
defects have been widely studied both in the past and recently, specially after the 
detection of the large-scale CMB anisotropics of cosmological origin that we call 
Their predictions have been confronted against a large bulk of observations and their 
performance is satisfactory, albeit not perfect. Both classes of models have their strong 
and weak points and we hope to have given a general flavour of them in Chapter 2 
(you can always consult [Linde, 1990] and [Vilenkin & Shellard, 1994] if not satisfied). 

The different sources of CMB anisotropics are then explained in detail in Chapter 
3, both on large as well as on small angular scales. We also give here a brief overview 
of the predictions from standard cold dark matter (CDM) plus inflationary scenarios 
and confront them with the present status of the observations on various scales. 

A quantitative account of the spectrum and the higher-order correlation functions 
of the CMB temperature anisotropics, together with the necessary introductory formu- 
lae is given in Chapter 4. Part of it, specially the analysis in §0, has been published 
in [Gangui, Lucchin, Matarrese & MoUerach, 1994]. 

Having given the general expressions for the spatial correlations of the anisotropics, 
we apply them in Chapter 5 to the study of the predictions from generalised models 
of inflation (i.e., a period of quasi-exponential expansion in the early universe). We 
use the stochastic approach to inflation [e.g., Starobinskii, 1986; Goncharov, Linde & 
Mukhanov, 1987], as the natural framework to self-consistently account for all second- 
order effects in the generation of scalar field fluctuations during inflation and their 
giving rise to curvature perturbations. We then account for the non-linear relation be- 



tween the inflaton fluctuation and the peculiar gravitational potential $, ending §|5Tl 



with the computation of the three-point correlation function for $. We then concen- 
trate on large angular scales, where the Sachs- Wolfe effect dominates the anisotropics, 
and compute (§ |5.2| ) the CMB temperature anisotropy bispectrum. From this, we cal- 
culate the collapsed three-point function (a particular geometry of the three-point 
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function defined in Chapter 4), and show its behaviour with the varying angle lag 
(§|5.2|). We specialise these results in §^]^ to a long list of interesting inflaton poten- 
tials. We confront these primordial predictions (coming from the non-linearities in the 
dynamics of the inflaton field) with the theoretical uncertainties in § |5.4| . These are 
given by the cosmic variance for the skewness, for which we show the dependence for 
a wide range of the primordial spectral index of density fluctuations. We end Chapter 
5 with a discussion of our results; what we report in this Chapter has been published 
in [Gangui, Lucchin, Matarrese & MoUerach, 1994] and [Gangui, 1994]. 

Even for the case of primordial Gaussian curvature fluctuations, the non-linear 
gravitational evolution gives rise to a non-vanishing three-point correlation function 
of the CMB. In Chapter 6 we calculate this contribution from the non-linear evolu- 
tion of Gaussian initial perturbations, as described by the Rees-Sciama (or integrated 
Sachs- Wolfe) effect. In § |6.2| we express the collapsed three-point function in terms of 
multipole amplitudes and we are then able to calculate its expectation value for any 
power spectrum and for any experimental setting on large angular scales. We also 
give an analytic expression for the rms collapsed three-point function (§ |6.3|) arising 
from the cosmic variance of a Gaussian fluctuation field. We apply our analysis to the 
case of COBEy-DMR in § |6.4| , where we also briefly discuss our results. The contents 
of this Chapter correspond to the study performed in [MoUerach, Gangui, Lucchin, & 
Matarrese, 1995]. 

The rest of the thesis is devoted to the study of possible signatures that topological 
defects may leave in the CMB radiation anisotropies. In Chapter 7 we introduce an 
analytic model for cosmic strings and calculate the excess CMB kurtosis that the 
network of strings produce on the relic photons after the era of recombination though 
the Kaiser-Stebbins effect (§ [7.2[ ). As we did in previous Chapters, here also we quantify 
a measure of the uncertainties inherent in large angular-scale analyses by computing 
the rms excess kurtosis as predicted by a Gaussian fluctuation field; we confront it 
with the mean value predicted from cosmic strings in § [7.3| . This Chapter contains 
essentially what was published in [Gangui & Perivolaropoulos, 1995]. 

In Chapter 8 we introduce another recently proposed analytic model for defects, 
but this time devised for treating textures. After a brief overview of the basic features 
of the model (§ |8.1| ) we calculate the multipole coefficients for the expansion of the 
CMB temperature anisotropies in spherical harmonics on the microwave sky (§ p.2|) . 
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At this point one may make use of the whole machinery developed in previous Chap- 
ters to compute texture-spot correlations (§8.^j|) as well as the angular bispectrum and, 
from this, also the collapsed three-point correlation function. We also estimate the 
theoretical uncertainties through the computation of the cosmic variance for the col- 
lapsed three-point function, for which we present an explicit analytic expression (§ ^.41) . 
The main aim of this study is to find out whether the predicted CMB non-Gaussian 
features may provide a useful tool to constrain the analytic texture model. The work 
presented in this Chapter is in progress. We provide, however, a preliminary discussion 
of the expected results in §8.5| . 

The last Chapter is devoted to the study of CMB anisotropics from textures, this 
time on small angular scales (~ 1°). As is well known, both inflation and defect 
models provide us with adequate mechanisms for the generation of large-scale structure 
and for producing the level of anisotropics that have been detected by COBE and 
other experimental settings on smaller scales. The small-scale angular power spectrum 
is most probably the first handle on the problem of discriminating between models, 
and this will be the focus of our study in Chapter 9. We first discuss ( §9.1| ) which 
contributions to the anisotropy effectively probe those perturbations with characteristic 
scale of order of the horizon distance at the era of recombination. In § |9.2| we set up 
the basic system of linearised perturbation equations describing the physics of our 
problem, as well as the global scalar (texture) field source terms. In § p73| we discuss 
the issue of initial conditions based on the causal nature of defects; we then compute the 
angular power spectrum and show the existence of the Doppler peaks (§ |9.4|) . Finally, 
we comment on the implications of our analysis in §^]5|. This Chapter is a somewhat 
slightly inflated version of what was reported in [Durrer, Gangui & Sakellariadou, 
1995], and recently submitted for publication. 

During my graduate studies here at SISSA I also did some work not directly related 
to the CMB radiation anisotropics [Castagnino, Gangui, Mazzitelli & Tkachev, 1993], 
a natural continuation of [Gangui, Mazzitelli & Castagnino, 1991] on the semiclassical 
limit of quantum cosmology; this research, however, will not be reported in this thesis. 



Chapter 1 

Standard cosmology: Overview 



1.1 The expanding universe 



One of the cornerstones of the standard cosmology is given by the large-scale homo- 
geneity and isotropy of our universe. This is based in the cosmological principle, a 
simple assumption requiring that the universe should look the same in every direction 
from every point in space, namely that our position is not preferred. 

The metric for a space with homogeneous and isotropic spatial sections is of the 
Friedmann-Robertson-Walker (FRW) form 



ds^ = dr - a\t) 



+ r^dO' + sin^ 9 



[l.l) 



where r, 9 and are spatial comoving coordinates, t is physical time measured by 
an observer at rest in the comoving frame (r, 6', 0) = const, and a{t) is the scale 
factor. After appropriate rescaling of the coordinates, k takes the values 1, —1 or for 
respectively constant positive, negative or zero three-space curvature. 

In 1929 Hubble announced a linear relation between the recessional velocities of 
nebulae and their distances from us. This is just one of the many kinematic effects that 
we may 'derive' from the metric (|1.1| ). If we take two comoving observers (or particles. 



if you want) separated by a distance d, the FRW metric tells us that this distance will 

grow in proportion to the scale factor, the recessional velocity given by f = Hd^ where 

H = a/a is the Hubble parameter (dots stand for t-derivatives) . Given the present 

uncertainties on its value|^, H is usually written as if = lOO/i km sec~^Mpc"^, with 

^The uncertainties are mainly due to the difRculty in measuring distances in astronomy [see, e.g., 
Rowan-Robinson, 1985]. 
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0.4 < h < 1. In addition to this, the expansion will make the measured wavelengths, 
Xobs, of the light emitted by stars, Xem, to redshift to lower frequencies by the usual 
Doppler effect: 1 + z = a{tobs) / o,{tem) = Xobs/ Km- This expression defines the redshift 
z, which according to observations is 2; > in the expanding universe. For objects 
with receding velocities much smaller than the velocity of light we have z ~ f/c and 
thus we can estimate their distances from us as d cz/H (assuming negligible peculiar 
velocities). 

Hence, we may calculate how far away nearby objects are (say, within our local 
group of galaxies). For instance , the Virgo clusteiQ has a redshift z ~ 0.004 which 
means a distance d ~ 12h^^ Mpc from us and a receding velocity v ~ 1150 km sec~^. 
Clearly distances as calculated above are valid for z << 1. When z ^ 1 Hubble's law 
is given by a somewhat more complicated formula [see, e.g., Kolb, 1994]. 

1.2 The Einstein field equations 

So far we have dealt just with the the properties of the FRW metric. For it to be 
an adequate representation of the line element, the Einstein field equations should be 
satisfied. These are 

Rf.u - \r9^u = SnGT^, , (1.2) 

where R^j^^, R, g^y and T^j, are the Riemann tensor, the Ricci scalar, the metric and 
the stress-energy-momentum tensor (for all the fields present), respectively^. 

In order to derive the dynamical evolution of the scale factor a(t), the form of 
T^i, must be specified. Consistent with the symmetries of the metric, the energy- 
momentum tensor must be diagonal with equal spatial components (isotropy). Thus 
T^u takes the perfect fluid form, characterised by a time-dependent energy density p{t) 
and pressure p{t) 

Tf,u = (p + p)u^u^ - pg^u ■ (1-3) 

^The case of the Andromeda nebula (M31) is different in two aspects: first, its distance d ^ 0.65 
Mpc can be determined by direct means and so the factor does not appear. Secondly, the 
negative velocity v ^ —270 km sec^^ {z ^ —0.0009) means that it is not receding from us but in fact 
approaching us (and thus its peculiar velocity wins, like in most of the other members of the local 
group of galaxies [Weinberg, 1972]. 

•^We are here not including the +A(7^^ on the right-hand side, corresponding to a cosmological 
constant (A) term. 
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Here, is the four-velocity of the comoving matter, with m° = 1 and -u* = 0. The 
local conservation equation for the energy-momentum tensor, T'^'l^ = yields 

p + 3H{p + p)=0, (1.4) 

where the second term corresponds to the dilution of p due to the expansion {H = a/a) 
and the third stands for the work done by the pressure of the fluid (first law). 

With regard to the equation of state for the fluid, we need to specify p = p{p). It 
is standard to assume the form p = wp and consider different types of components 
by choosing different values for w. In the case w = const we get p oc a~^^^~^'^'> (use 
Eq. ( |1.4D ). If the universe is flUed with pressureless non-relativistic matter ('dust') 
we are in the case where p « p and thus p oc a~^. Instead, for radiation, the ideal 



-4 



relativistic gas equation of state p = p/3 is the most adapt, and therefore^ p oc a' 
Another interesting equation of state is p = — p corresponding to p = const. This is the 
case of 'vacuum energy' and will be the relevant form of energy during the so-called 



inflationary epoch; we will give a brief review of inflation in § p.l| below. Note from the 
different ways in which dust and radiation scale with a{t) that there will be a moment 
in the history of the expanding universe when pressureless matter will dominate (if 
very early on it was radiation the dominant component - as is thought to be the case 
within the hot big bang). In fact, pmat/ Prad a{t). This moment is known as the time 
of 'equality' between matter and radiation, and is usually denoted as t^g. 

The main fleld equations for the FRW model with arbitrary equation of state (ar- 
bitrary matter content) are given by 

3H + 3H^ = -^8nG{p + 3p) , (1.5) 

namely, the Raychaudhuri equation for a perfect fluid (with vanishing shear and vor- 
ticity, no cosmological constant and where the fluid elements have geodesic motion), 
and the Friedmann equation 

3H^ + 3K = SttGp . (1.6) 
This last can be seen as the equation satisfled by the Ricci scalar = 6K of the 



^Consider for example photons: not only their density diminishes due to the growth of the volume 
(cx a^^), but the expansion also stretches their wavelength out, which corresponds to lowering their 
frequency, i.e., they redshift (hence the additional factor oc a~^). 
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three-dimensional space slices (again for vanishing shear and cosmological constant). [| 
The spatial curvature is written as K = k/a^ij:), with k constant^ (cf. Eq. (|0|)), and 
has time evolution given by K = —2HK. 

It is convenient at this point to define the density parameter Q = p/pc, where 
Pc = SH"^ / (SnG) is the (critical) density necessary for 'closing' the universe.[] In fact, 
from Eq. ( |1.6|) we get the relation to be satisfied by Q, 

^=1 + ^4^, (1.7) 

and so, for k = —1 (open three-hypersurfaces) we have < 1, while k = 1 (closed 
three-hypersurfaces) corresponds to f2 > 1 with the fiat k = case given by Q = 1 
(transition between open and closed t = const slices). 

From ( p..7|) we see that, when values other than the critical one Q = 1 are considered, 
we may write = K/ (Q — l), and from this = k/{Q — l), for k ^ 0. These relations 
apply for all times and all equations of state (including an eventual infiationary era) 
provided Q includes the energy density of all sources. We take p > 0, which is the 
case for ordinary matter and scalar fields. From the Raychaudhuri equation we may 
express the deceleration parameter q = —a/{aH'^) as follows 

q=^-(^- + w)n, (1.8) 

where we recall w = p/p. This shows that w = — | is a critical value, separating 
qualitatively different models. A period of evolution such that w < — | (g < 0), and 
hence when the usual energy inequalities are violated, is called 'inflationary'. 

From ( pTsD and the conservation equation (|L^) we get the evolution of the density 
parameter [Madsen & Ellis, 1988] 

3, 1,^ , ^ \ ^ X 

= -{w + -)Q{Q - 1) . (1.9) 
da a 6 



^The Friedmann equation (1.6) follows from the 0-0 Einstein equation. Combining Eq. ( |l.6| ) and 
the i-i Einstein equation we get (1.5). 

^If no point and no direction is preferred, the possible geometry of the universe becomes very 
simple, the curvature has to be the same everywhere (at a given time). 

^The present value of the critical density is pc — l.88h^ x 10^^^ g cm~'^, and taking into account 
the range of permitted values for h, this is ~ (3 — 12) x 10~^^ kg m~^ which in either case corresponds 
to a few hydrogen atoms per cubic meter. Just to compare, a 'really good' vacuum (in the laboratory) 
of 10~^ N/m^ at 300 K contains - 2 x 10" molecules per cubic meter. The universe seems to be 
empty indeed! 
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Note that ^7 = 1 and = are solutions regardless of the value of w. This and the 
previous equations imply that if > 1 at any time, then it will remain so for all times. 
The same applies for the open case: if < i7 < 1 at any given time, this will be true 
for all times. We will make use of Eq. ( |1.9D below, when we will come to consider the 
drawbacks of the standard cosmology in § |1.6i 



1.3 Thermal evolution and nucleosynthesis 

According to the standard hot big bang model, the early universe was well described 
by a state of thermal equilibrium. In fact, the interactions amongst the constituents of 
the primordial plasma should have been so effective that nuclear statistical equilibrium 
was established. This makes the study simple, mainly because the system may be fully 
described (neglecting chemical potentials for the time being) in terms of its temperature 
T. In a radiation dominated era the energy density and pressure are given by p = 3p = 
TT^^f^.T'^/SO, where the T-dependent g^, gives the effective number of distinct helicity 
states of bosons plus fermions with masses << T. For particles with masses much 
larger than T (in natural units) the density in equilibrium is suppressed exponentially 
[see, e.g., Weinberg, 1972; the notation however is borrowed from Kolb & Turner, 
1990]. 

With the expansion the density in each component diminishes and so it gets more 
difficult for the above mentioned effective interactions to keep 'working' as before. 
It is clear that there will be a point at which the interaction rate F of a certain 
species with other particles will fall below the characteristic expansion rate, given 
by H. Approximately at that moment that species is said to 'decouple' from the 
thermal fluid, and the temperature T^ec at which this happens is called the decoupling 
temperature^. 

Let us follow the history of the universe backwards in time. Very early on all 



^Massless neutrinos decouple when their interaction rate ~ G^T^ gets below H ~ T^/mp (Gp 
is the Fermi constant, characteristic strength of the weak interactions). This happens around T^^^ ~ 1 
MeV. Photons, on the other hand, depart from thermal equilibrium when free electron abundance is 
too low for maintaining Compton scattering equilibrium. This happens around T]^^ ~ 0.3 eV. (If you 



want to know why this is smaller than the binding energy of hydrogen, 13.6 eV, consult Ref. |232|) 



In the Chapters to come, by T^ec we will always mean Tj^^, namely the decoupling temperature for 
the CMB radiation, the 'messenger' that gives us direct information of how the universe was 400,000 
years after the bang. 'Direct' information of earlier times we cannot get (. . . if only we could detect 
neutrinos as easily as we can detect CMB photons . . .) 
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the matter in the universe was ionised and radiation was the dominant component. 
With the expansion, the ambient temperature cools down below ~ 13.6 eV and the 
recombination of electrons and protons to form hydrogen takes place, which diminishes 
the abundance of free electrons, making Compton scattering not so effective. This 
produces the decoupling of CMB radiation from matter and, assuming the universe is 
matter-dominated at this time, this occurs at about tdec ~ 5.6 x 10^'^{nh'^)~^/'^ sec. 
The corresponding redshift is about Zdec — HOO and the temperature T^ec — 3000 K 
= 0.3 eV. 

Going still backwards in time we reach the redshift z^q — 2 x W^flh^ when matter 
and radiation (i.e., very light particles) densities are comparable. By this time the 
universe was teg ^ 1.4 x 10^(n/i^)~^ years old, with a temperature Tgg ~ S.SQ/t^ eV 
and a density peq ~ 3.2 x 10-^^{QhY g/cm^. 

At even earlier times we reach densities and temperatures high enough for the 
synthesis of the lightest elements: when the age of universe was between 0.01 sec and 3 
minutes and its temperature was ar ound 10 - 0.1 McV the synthesis of D, '^Hc, ^Hc, ^Li 
took place. The calculation of the abundance of these elements from cosmological origin 
is one of the most useful probes of the standard hot big bang model (and certainly the 
earliest probe we can attain) [refer to Malaney & Mathews, 1993 for a recent review]. 

The outcome of primordial nucleosynthesis is very sensitive to the baryon to photon 
ratio 1] = ribar/nphoton, and agreement with all four observed abundances limit it in 
the range ~ (3 — 5) x 10^^°. Given the present value rj ^ 2.68 x lO'^Qsh'^ , this 
implies the important constraint 0.011 < fi^^^ < 0.019 [Walker et ah, 1991; Smith et 
al, 1993]. The consequences of this are well-known. Given the uncertainties on the 
Hubble parameter (buried in h) we derive 0.011 < Qb ^ 0.119. Recall that luminous 
matter contribute less than about 0.01 of critical density, hence there should be 'dark' 
matter (in particular, dark baryons). On the other hand, dynamical determinations of 
Q point towards ~ 0.15 — 0.2. This implies there should also exist 'non-baryonic' 
dark matter. 

The compatibility between nucleosynthesis predictions and the observed abun- 
dances is one of the successes of the hot big bang model and gives confidence that 
the standard cosmology provides an accurate accounting of the universe at least as 
early as 0.01 sec after the bang. 
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1.4 The CMB radiation and gravitational instabil- 
ity 

The CMB radiation provides another fundamental piece of evidence in favour of the hot 
beginning of our universe. After its discovery in 1965, the next feature that surprised 
people was its near-perfect black-body distribution (over more than three decades in 
wavelength A ~ 0.03-100 cm) with temperature ~ 2.7 K. Recently the COBE-FIRAS 
detector measured it to be 2.726 ± 0.01 K [Mather et al, 1994]. 

Once decoupled, the background radiation propagates freely, keeping its Planck 
spectrum and redshifting as T oc (1 + 2;). If the universe became reionised at a lower 
redshift (e.g., due to early star or quasar formation) then the 'last scattering surface' 
may be closer to u^. Once that we know the temperature of the relic radiation we 
may easily compute its number density and energy density [e.g., Kolb, 1994] 

~ 411 cm"^ ; ~ 4.71 x 10"^^g cm"^ . (1.10) 

In 1992 the COBE collaboration announced the discovery of the long sought-after 
anisotropics on angular scales ranging from about 7° through to 90°, at a magnitude 
of about 1 part in 10^ [Smoot et al, 1992]. One way of characterising the level of 
anisotropics is by the rms temperature variation, which the COBE team found to be 
30 ± 5/iK on a sky averaged with a 10° FWHM beam. The anisotropy of the CMB 
radiation will be the topic of Chapter |^ below. 

Let us move on now to the issue of large-scale structure formation. The favourite 

picture today is that of structures developing by gravitational instability from an initial 

spectrum of primordial density perturbations. One usually expands all quantities (like 

the density perturbation) in Fourier modes (we are working in fiat space). Once a 

particular mode becomes smaller than the horizon two competing effects will determine 

the future of the fluctuation. The dynamical time scale for gravitational collapse is 

Tgrav ~ (Gp)^^^^ and, uulcss we consider effectively pressureless fluids, there will be 

an analogue characteristic time of pressure response given by Tpress ~ A/fs, where A is 

the physical wavelength of the perturbation and Vg is the sound speed of the fluid. If 

Tgrav ^ 'tpress then prcssurc forccs cannot overcome the gravitational attraction and the 

^Throughout this thesis we will be considering scenarios where early reionisation does not take 
place [see, e.g., Sugiyama, Silk & Vittorio, 1993 for the effects that reionisation has on the CMB 
anisotropics] . 
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collapse is inevitable. This occurs for A > Aj ~ Vs/ {GpY^"^ . Hence, the Jeans length, 
A J, and associated mass, Mj ~ 7rpAj/6, define the scales above which structures 
become unstable to collapse .Q 

The subsequent evolution of fluctuations depend very much on the kind of matter 
that dominates after tgg- The flrst type of matter we may think of is, of course, baryonic 
dark matter. Photon difussion [Silk, 1968] plays a relevant role in this scenario, since it 
will dissipate small-scale fluctuations. Hence, large structures (in the form of pancakes) 
will form flrst; galaxies and smaller structure will be formed from fragmentation of 
these pancakes afterwards. This model, however, runs into problems since it cannot 
get the structure we now observe formed without generating too much anisotropics in 
the CMB radiation. This is mainly due to the fact that initial perturbations in the 
baryon component can only begin to grow after the era of recombination. Before that, 
baryons are not free to move through the radiation plasma to collapse. Hence, having 
'less' time to grow to a certain 'flxed' level, the amplitude of fluctuations at horizon 
crossing have to be much larger. As we said, this generates too much anisotropics 
[Primack, 1987]. 

If baryons will not help, cold dark matter (CDM) particles may be possible candi- 
dates; roughly, these are relics with very small internal velocity dispersion, namely, 
heavy particles that decouple early in the history of the universe and are non- 
relativistic by now[^ Other candidates are hot dark matter (HDM) particles: these are 
light particles (m < 100 eV) that decouple late and are still relativistic when galactic 
scales cross the horizon. Example of this are massive neutrinos with irty a few tens of 
eV. 

In a CDM universe and due to the small growth of the perturbations that takes place 
between horizon crossing and the time of equality (for scales less than Aeg — 13(fi/i^)~^ 
Mpc) the density contrast 5p/ p increases going to smaller scales. Then the flrst objects 
to form are of sub-galactic size leading to the 'botton-up' hierarchical scenario. These 
small-mass systems are subsequently clustered into larger systems that become non- 
linear afterwards. The hierarchical clustering begins with masses of order the baryon 

^°We will define the mass associated to a given perturbation scale in the next section. 

^^The possible exception being the axion, the Goldstone boson arising from a U{\) global Peccei- 
Quinn symmetry breaking phase transition, that acquires a small mass below Aqcd ~ few MeV, due 
to instanton efi'ects. In many simple models the axion has never been in thermal equilibrium [see, 
e.g., Turner, 1990 and Raffelt, 1990 for reviews]. 
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Jeans mass at recombination {Mj{trec) ~ IO^Mq) and continues until the present. 

The situation is very different in the case of HDM. As we mentioned before, imagine 
we have neutrinos of mass rrii, ~ a few tens of eV. They become non-relativistic ap- 
proximately at teq- However, small-scale fluctuations are prevented from growing due 
to 'free-streaming', namely, the high thermal velocities endowed by neutrinos make 
them simply stream away from the overdense regions; in so doing they erase the fluc- 
tuations. We can define an 'effective' Jeans length to this effect, which we will call the 
free-streaming length, Xf^; the corresponding mass is given by Mfg ~ IO^^Mq. Hence, 
in order for perturbations (in a HDM scenario) to survive free-streaming and collapse 
to form a bound structure, the scale has to be that corresponding to superclusters. 
Galaxies and smaller structure form by fragmentation in this 'top-down' picture. 

1.5 The perturbation spectrum at horizon crossing 

It is standard in any treatment of the evolution of cosmic structures to assume a 
statistically homogeneous and isotropic density field p(x), and its fluctuations to be 
the seed needed by gravity for the subsequent clumping of matter. 

These fluctuations are given by the density contrast 6{x) = (p(x) — (p))/(p), and 
are calculated as the departures of the density from the mean value (p). By definition 
the mean value of 6 over the statistical ensemble is zero. However, its mean squared 
value is not, and it is called the variance cx^, representing a key quantity in the analysis. 

It is straightforward to give an expression for cr^ in terms of the power spectral 
density function of the field 6 (or, more friendly, power spectrum) P{k), 



The variance gives no information of the spatial structure of the density field. It 
however depends on time due to the evolution of the Fourier modes 5k oc P^/^(fc), and 
therefore encodes useful information on the amplitude of the perturbations. 

The second equality of Eq. ( |1 . 1 Ij ) defines V{k) = P(fc)fc^/(27r^), representing the 
contribution to the variance per unit logarithmic interval in k: this quantity lends itself 
well for comparison of large-scale galaxy clustering [Coles & Lucchin, 1995]. 

The variance as given above suffers one main drawback: it contains no information 
on the relative contribution of the different modes. It may even yield infrared or 
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ultraviolet divergencies, depending on the form of P{k). In practice what people do 
is to introduce some kind of resolution scale, say r, in the form of a window function, 
which acts as a filter smearing information on the modes smaller than r. The mean 
mass contained inside a spherical volume ~ is (M) ~ (p)r^. One can then define 
the mass variance in this volume as 

< = {[-Jfj ) = P{k)W\kr)k'dk. (1.12) 

The second equality follows after some straightforward steps [Kolb & Turner, 1990]. 
(Tm is the rms mass fluctuation on the given scale r and thus depends on this scale, 
and through it on the mass M. As already mentioned, the window function W'^{kr) 
makes the dominant contribution to o"m to come from wavelengths greater than r.Q 

Spurious effects (leading to incorrect results for ctm) due to boundaries in the win- 
dow function may be prevented by taking a Gaussian ansatz W'^{kr) ~ exp(— fc^r^). If 
we also consider a featureless power spectrum P{k) oc A;", with n the so-called spectral 
index of density perturbations, we easily get ctm oc r^i'^+^)/'^^ with the amplitude of aM 
posing no problems provided we take n > —3. 

Let us write this result for the rms mass fluctuation as ctm (t) oc M-("+3)/6, where 
the t-dependence is to emphasize the fact that we are computing the fluctuation at a 
particular time t, in contrast to the one at horizon crossing time, as we will see below. 
Of particular importance is the choice n = 1, and we can see heuristically why this 
is so: The perturbation in the metric, as given by the gravitational potential $, may 
be expressed as $(r) ~ G 6M/r ~ G 6p{r) ~ Gpa^ir"^ oc M"("+^)/^M^/^, and for 
n = 1 this result is independent of the mass scale M, hence the name scale invariant 
spectrum [Peebles & Yu, 1970; Harrison, 1970; Zel'dovich, 1972]. 

It is useful to define the mass associated to a given perturbation A as the total mass 
contained within a sphere of radius A/2, i.e., M = -IpA^. According to this, we may 
write the horizon mass in cold particles as Mh = ^PcdmH~^, and given that the density 
in the cold component scales as pcdm oc (1 + z)^ and that H^^ oc a^(t) oc (1 + z)~'^ 
{H^^ oc a'^/^(t) oc (1 + z)~^^'^ ) during radiation (matter) domination, we get Mh oc 
(1 + z)-^ {Mh oc (1 + z)~^/2 ) for z > z^g {z < z^g and z » Q-'^). 

The horizon crossing time (or redshift zh{M)) of a mass scale M is commonly 

^^This is a generic property of any W'^{kr), regardless of the specific form, i.e., top-hat, Gaussian, 
etc. 
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defined as that time or redshift at which M coincides with the mass inside the horizon, 
Mh- Thus we have Mh{zh{M)) = M. Prom our discussion of the last paragraph 
we easily find the dependence on the scale M of the horizon crossing redshift, namely 

zh{M) oc [zh{M) oc for zh{M) > Zeq {zh{M) < Zeq and Zh{M) » 

Let us turn back now to the rms mass fiuctuation ctm- Any perturbation on super- 
horizon scales grows purely kinematically from the time t of its generation until the 
time it enters the particle horizon, and its mass fiuctuation is given (at least in the 
liner regime) by ctm oc t oc (1 + ^)~^ (ctm oc t^/^ oc (1 + z)~^) before (after) the time of 
equivalence. Then, we have (for zh{M) > Zeq) 



while for zh{M) < z^q we have 

aM{zH{M)) ~ aM{t) 



1 + Zeq\ 'fl + ZH{My 



1 + z{t) I \ l + z, 



eq 



RecaUing now that the fiuctuation at its origin t was given by (Tm(^) oc M^^^+^^Z^ 
we find in both cases that the rms mass fiuctuation at the time of horizon crossing 
is given by aM{zH{M)) oc M^^"'^^/^. This shows that mass fiuctuations with the 
Harrison-Zel'dovich spectrum {n — 1) cross the horizon with amplitudes independent 
of the particular scale. 

1.6 Drawbacks of the standard model 

In previous sections we have briefiy reviewed the standard hot big bang cosmology 
and emphasised its remarkable successes in explaining a host of observational facts: 
among these, the dynamical nature (the expansion) of our universe, the origin of the 
cosmic microwave background from the decoupling between matter and radiation as 
a relic of the initial hot thermal phase. It also provides a natural framework for 
understanding how the large-scale structure developed, and for the accurate prediction 
of light elements abundance produced during cosmological nucleosynthesis. 
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In this section we will focus, instead, on its shortcomings (or at least, some of them). 
These are seen not as inconsistencies of the model but just issues that the model cannot 
explain, when we extrapolate its highly accurate predictions back in time, beyond, say, 
the era of nucleosynthesis and before. We leave for § ^.1.3| the discussion on how the 
inflationary scenario yields a well-defined, albeit somewhat speculative, solution to 
these drawbacks, based on early universe microphysics. We are not discussing here the 
'cosmological constant' problem, namely, why the present value of A (or equivalently, 
the present energy of the vacuum) is so small compared to any other physical scalej^; 
this mainly because it keeps on being a problem even after the advent of inflation: 
in fact inflation makes use of the virtues of a vacuum dominated period in the early 
universe. While offering solutions to the problems listed below, inflation sheds no light 
on the problem of the cosmological constant [see Weinberg, 1989 for an authoritative 
review; also Carrol et al, 1992 and Ng, 1992 for more recent accounts]. 

The horizon (or large— scale smoothness) problem 

The relic CMB radiation we detect today was emitted at the time of recombi- 
nation. It is uniform to better than a part in 10^, which implies that the universe 
on the largest scales (greater than, say, 100 h^^ Mpc) must have been very smooth, 
since otherwise larger density inhomogeneities would have produced a higher level of 
anisotropics, which is not observed. The existence of particle horizons in standard cos- 
mology precludes microphysic events from explaining this observed smoothness. The 
causal horizon at the time of the last scattering subtends an apparent angle of order 
2 degrees today, and yet the radiation we receive from all directions of the sky have 
the same features. If there was no correlation between distant regions, how then very 
distant (causally disconnected) spots of the sky got in agreement to produce the same 
radiation features, e.g., the same level of anisotropics and temperature? This is just 
one way of formulating the horizon problem. 

The fiatness/oldness problem 

^•^In fact, observational limits impose p\ = A/{8ttG) to be smaller than the critical density, and 
so \pa\ < lO""** GeV^, in natural units, corresponding to |A| < 10~^^ cm~^. Alternatively, if we 
construct m\ = [\pA\i^/c)^]^^^ which has dimensions of mass, the above limits constrain m\ < 10"'^^ 
eV. These unexplained very small values are undesired in cosmology, hence the A-problem. 
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We saw in § O how Eq. ( [1.9| ) implies that the density parameter f2 does not remain 



constant while the universe expands, but instead evolves away from 1. Assuming the 
standard evolution according to the big bang model extrapolated to very early times 
and given that observations indicate that Q is very close to 1 today, we conclude 
that it should have been much closer to 1 in the past. Going to the Planck time 
we get |i7(10~^^sec) — 1| ~ (9(10~^°), and even at the time of nucleosynthesis we get 
|fi(lsec) — 1| ~ 0(10^^^). These very small numbers are but one aspect of the flatness 
problem. 

This implies that the universe was very close to critical in the past and that 
the radius of curvature of the universe acurv = H^^/\H — 1|^/^ = a/\/k was much 
much greater than the Hubble radius H^^. Were this not the case, i.e., suppose 
[H-^/a^urv? = {\k\/a^)/{87TGp/3) ~ 0(1) at the Planck time, then, if closed {k > 0) 
the universe would have collapsed after a few Planck times (clearly not verified - our 
universe is older than this), while if open {k < 0) it would have become curvature dom- 
inated, with the scale factor going like a oc t (coasting), and cooling down so quickly 
that it would have acquired the present CMB radiation temperature of ~ 3 K at the 
age of 10~^^ sec, which clearly is at variance with the age of the universe inferred from 
observations [Kolb & Turner, 1990]. This is just another face of the same problem, 
namely the difficulty in answering the question of 'why' the universe is so old. 

The unwanted relics problem 

Baryogenesis is an example of the virtues of getting together cosmology and grand 
unified theories (GUTs). However, the overproduction of unwanted relics, arising from 
phase transitions in the early universe, destroys this 'friendship'. The overproduction 
traces to the smallness of the horizon at very early times. Defects (e.g., magnetic 
monopoles) are produced at an abundance of about 1 per horizon volume (cf. the 
Kibble mechanism, see § p.2.2| ). This yields a monopole to photon ration of order 



{Tout I'fnp)^ and a present fimonop far in excess of 1, clearly intolerable cosmologically 
speaking; see §p.2.4| for details. 



We will see in § |2.2| below that the simplest GUTs generically predict the formation 
of topological defects, and why many of these defects are a disaster to cosmology; in 
particular the low-energy standard model of particle physics cannot be reached (as 
the last step of a chain of phase transitions from a larger symmetry group) without 
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producing local monopoles [Preskil, 1979]. The standard cosmology has no means of 
ridding the universe of these overproduced rehcs; hence the problem. 



Chapter 2 

Theories of the early universe 



A particularly interesting cosmological issue is the origin of structure in the universe. 
This structure is believed to have emerged from the growth of primordial matter- 
density fluctuations amplified by gravity. The link of cosmology to particle physics 
theories has led to the generation of two classes of theories which attempt to provide 
physically motivated solutions to this problem of the origin of structure in the universe. 

The first class of theories are those based on a mechanism called inRation; ac- 
cording to it, the very early universe underwent a brief epoch of extraordinarily rapid 
expansion. Primordial ripples in all forms of matter-energy perturbations at that time 
were enormously amplified and stretched to cosmological sizes and, after the action of 
gravity, became the large-scale structure that we see today. The initial idea that an 
early epoch of accelerated expansion would have interesting implications for cosmology 
is due to Guth [1981]. 

According to the second class of theories, those based on topological defects, pri- 
mordial fluctuations were produced by a superposition of seeds made of localised dis- 
tributions of energy density trapped during a symmetry breaking phase transition in 
the early universe. This idea was first proposed by Kibble [1976], and has been fully 
worked out since then by many authors. 

This Chapter is devoted to a general review of these two competing models of large- 
scale structure formation and CMB anisotropy generation. The Chapter is divided into 
two 'big' sections, each of which deals with one of these scenarios. Let us begin first 
(in alphabetical order) with inflationary models. 
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2.1 Inflation 

As we have discussed in Chapter |I] the standard big bang cosmology is remarkably 
successful. It explains the expansion of the universe, the origin of the CMB radiation 
and it allows us to follow with good accuracy the development of our universe from 
the time of nucleosynthesis (around 1 sec after the bang) up to the present time (~ 15 
Gyr). However, it also presents some shortcomings, namely, the flatness/oldness, the 
horizon/large-scale smoothness, the cosmological constant, and the unwanted relics 
problems. 

Until the 1980 's the was seemingly no foreseeable solution to these problems. With 
the advent of the inflationary scenario the way of thinking the early universe changed 
drastically. The scenario, which is based on causal microphysical events that might 
have occurred at times ~ 10~^^ sec in the history of the universe (and well after the 
Planck era ~ 10~^^ sec), offers a framework within which it is possible to explain some 
of the above mentioned problems. 

In the following sections we give a brief account of inflation. First we review a 
bit of the history leading to the inflationary idea, namely that it was worthwhile to 
study a period of exponential expansion in the early evolution of our universe. Then 
we concentrate in the basic facts related to the dynamics of the inflaton field, responsi- 
ble for driving the quasi-exponential expansion. We then show how inflation explains 
the handful of cosmological facts with which the standard model alone cannot cope. 
From this we go on to one of the nicest features that inflation predicts, namely, the 
origin of density and metric perturbations, whose consequences for structure formation 
and gravitational wave generation the reader may already appreciate. After a succinct 
explanation on how these perturbations evolve, we move on to a brief survey of infla- 
tionary models, trying to put them in perspective and explaining what features make 
them attractive (or why they proved unsuitable) as a realisation of the scenario. We 
finish the section with a short account of the stochastic approach to inflation. The 
results of this will be helpful in understanding the generation of space correlations in 
the primordial inflaton field that, on horizon entry, will lead to non-Gaussian features 
in the CMB radiation on large angular scales (cf. Chapter |^ below). 
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2.1.1 The paradigm: some history 

Soon after entering the subject one realises that currently there are many inflation- 
ary universe models, and some of them involve very different underlying high energy 
physics. Of course, something in common among them all is the existence of an early 
stage of exponential (or quasi-exponential) expansion while the universe was domi- 
nated by vacuum energy^ and filled just by almost homogeneous fields and nearly no 
other form of energy. This feature, shared by all simple models, is what sometimes 
goes under the name of the inflationary paradigm. The end of inflation is signaled by 
the decay of the vacuum energy into lighter particles, the interaction of which with 
one another leads to a state of hot thermal equilibrium. From that point onwards the 
universe is well described by the hot big bang model. 

The inflationary models (while not yet carrying this name) have a long history that 
traces back to the papers by Hoyle, Gold, Bondi, Sato and others [Lindley, 1985, cited in 



Ref. ||112|| , see also Gliner, 1965, 1970], where the possible existence of an accelerated 
expansion stage was first envisaged. A few years later Linde [1974, 1979] came up 
with the idea that homogeneous classical scalar fields (which appear in virtually all 
GUTs) could play the role of an unstable vacuum state, whose decay may give rise to 
enormous entropy production and heat the universe. Soon afterwards it was realised 
that quantum corrections in the theory of gravity also led to an exponential expansion: 
the Starobinskii model [Starobinskii, 1979]. 

But the major breakthrough came with Guth [1981] paper, where the real power of 
inflation for resolving the shortcomings of the standard model was spelled out. In spite 
of the fact that his original model did not succeed (essentially due to the extremely 
inhomogeneous universe that was produced after the transition - see § |2.1.6| below) 



it laid down the main idea and it took just months for people to propose the 'new' 
(or 'slow-rollover') scenario where the drawbacks of Guth's model were solved [Linde, 
1982a; Albrech & Steinhardt, 1982]. However, in the search for simple and powerful 
models, there is little doubt that the first prize goes to Linde's chaotic models [Linde, 
1983a]. In these models the inflaton is there just to implement inflation, not being part 
of any unified theory. Moreover, no special potential is required: an ordinary A0'^ or 
even a mass autointeraction term will do the job. There is neither phase transition nor 



^In most models the way of implementing the idea of a vacuum energy dominated universe is by 
assuming a scalar field whose initial state was displaced from its true vacuum (lowest energy) state. 



22 



Chapter 2. THEORIES OF THE EARLY UNIVERSE 



spontaneous symmetry breaking process involved, and the only requirement is that the 
inflaton field be displaced from the minimum of its potential initially. Different regions 
of the universe have arbitrarily different ('chaotically' distributed) initial values for 0. 

2.1.2 Scalar field dynamics 

The basic idea of inflation is that there was an epoch early in the history of the universe 
when potential, or vacuum, energy was the dominant component of the energy density. 
The usual way to realise an accelerated expansion is by means of forms of energy other 
than ordinary matter or radiation. Early enough energy cannot be described in term 
of an ideal gas but it should be described in terms of quantum flelds. Usually a scalar 
field is implemented to drive this infiationary era and in the present section we will 
consider its dynamics in detail. 

Let us study now the relevant equations for a homogeneous scalar field with 
effective potential V{(f)) in the framework of a FRW model, in presence of radiation, 
Pr(t). The classical equations of motion [Ellis, 1991] are given by the Friedmann 
equation 

H' + K^ ^SnG[V{(l>) + + ^ (2.1) 
the Raychaudhuri equation 

H + H' = ^8t:G[V{<P) - 0^ - p,] , (2.2) 

the energy conservation equation for 

f^[l^' + V{cP)]^-3H^'-S , (2.3) 

which is equivalent to the equation of motion for (p, and the energy conservation for 
radiation 

^Or = -^Hpr + 6 , (2.4) 

at 

where 5 accounts for the creation of ultra-relativistic particles (radiation) due to the 
time variation of the scalar field. S may be expressed as [Kolb & Turner, 1990] 

S^r^, (2.5) 

where the characteristic time F"^ for particle creation depends upon the interactions of 
(f) with other fields. As we are mainly interested in the phase of adiabatic evolution of 
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the scalar field, where F << 3H and thus particle creation processes are not operative, 
we will neglect the 6 term in what follows. When arriving at the reheating phase (at 
the very end of inflation) the damping of the coherent oscillations of the scalar field 
will lead to the creation of light particles, which, after thermalisation, will heat the 
universe to a temperature appropriate for continuing its evolution in the radiation era. 

It is a common practice, when dealing with the above equations, to assume a 
particular kind of rolling condition for the scalar field. The slow-rolling approximation 
101 << 3if0 yields SHc/) + dV/dcj) = and this describes the evolution of during 
inflation. On the other hand at the end of inflation the situation reverts itself and we 
have + dV/d(j) + TcfP' ~ 0. This signals the fast-rolling evolution during reheating. 

Although these are the approximations that one commonly finds in the literature, 
we will here follow the analysis of Ellis [1991] and study the dynamics of the field 
exactly. Combining Eqs. ( |2.1| ), (p^) (now without radiation density contribution) and 
assuming 7^ we have 

V^(0) = {MG)-\H + + 2K] , (2.6) 

and 

02(t) = {A^G)-^\-H + K] . (2.7) 

To specify a model we need only proceed as follows. First, choose a value for k (cf. 
Eq. ( |1.7|) ) and the initial value for 0. Secondly we specify the time-dependent scale 
factor a{t) (such that the right-hand side of Eq. (p.7|) is positive) and compute from 



this the Hubble parameter and its derivative. Use then ( p.7|) to get 0(t), integrate it to 
find 0(t) and finally invert it to get t(0). Finally, plug this into ( p.6| ) to get V{(j)). On 
the other hand, if = 0, equation ( |2.3|) tells us that dV/dcp ~ and thus the potential 
is flat. In summary, in all cases it is possible to find (in principle) the form of the 
inflationary effective potential satisfying the classical inflationary equations (without 
any assumption on the way 'rolls') with the desired functional form for the scale 
factor a(t), the chosen curvature k and initial condition for the scalar field. 

A specific example was given by Lucchin & Matarrese [1985]. Taking flat spatial 
sections = they chose a power-law behaviour for the scale factor a{t) cx t^, with 
constant p > 1 [Abbott & Wise, 1984b]. They studied the evolution of the system of 
classical equations from an initial time ti with 0(ti) = 0j 7^ and assumed (as it is 
usually done) that after a brief period of inflation the radiation contribution to the 
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total kinetic and potential energy is depressed, and so we may neglect the contribution 
of pr to the relevant equations. With this ansatz for the scale factor we readily get 

P t-\ (2.8) 



which gives 



0(t)^</),±^/^ln(t/t,) . (2.9) 



Inserting the solution with the plus sign into ( ^.1[ ) with K = 0, and making use of 
we get 

It should be said that this potential is to be considered as an approximation (for the 
evolution of the infiaton until before the time of reheating) of a more complex potential. 
In particular the featureless exponential potential cannot provide an oscillatory end to 
inflation. Power-law inflation, however, is of particular interest because exact analytic 
solutions exist both for the dynamics of inflation and for the density perturbations 
generated. 

Slow— rolling 

As we mentioned above, in a universe dominated by a homogeneous scalar field (j), 
minimally coupled to gravity, the equation of motion is given by ( p.3| ), which we may 
write 

4) + 3H^ + dV/d(f) = . (2.11) 

The energy density and pressure (neglecting radiation or other matter-energy compo- 
nents) are given by 

p^l^(0) + i02 ; p^_r(0) + l02^ (2.12) 

We then see that provided the field rolls slowly, i.e., 0^ < V, we will have 3p < —p 
and, from equation (|2.2| ) or equivalently a = — 47rG'(p + 3p)a/3, we have a > 0, namely, 
an accelerated expansion. 

In most of the usually considered models of inflation there are three conditions 
that are satisfied. The first is a statement about the solution of the field equation, and 
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says that the motion of the field in overdamped, namely that the force term dV/d(j) 
balances the friction term 3H(j), and so 

^ 1 (2.13) 



3H dcj) 

The following two conditions are statements about the form of the potential. The first 
one constrains the steepness of the potential (or better, its squared) 



ml (V 



-77 \ «1> (2-14) 



where ' = djdc^). This implies that the condition 0^ < V is well satisfied, and ~ 
87rG'V(</))/3 which implies that the Hubble parameter is slowly varying. The third 
condition that is usually satisfied is 

(2.15) 

and is independent of the other two conditions [Liddle & Lyth, 1993]. 
2.1.3 Problems inflation can solve 

We have seen before that the standard model, although extremely successful in its 
predictions, left some open questions (cf. §1.6| ). In this section we review some com- 
mon features of infiationary cosmologies and make the way to understand how these 
problems may be solved. 

Let us first consider the horizon problem. We saw in the previous section that we 
may define the accelerated expansion of the universe by the condition a > 0. We may 
calculate now 

which, together with a > implies d{H~^/ a)/dt < 0, namely that the comoving Hubble 
length rn = H^^/a is decreasing during any accelerated expansion period. This is the 
basic mechanism through which infiation solves the horizon problem. 

At the beginning of the infiationary period the comoving Hubble length is large and 
therefore perturbations on very large scales (like our present horizon scale, and beyond) 
are generated causally within r^- With the accelerated expansion rn decreases to such 
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an extent that its subsequent increase during the radiation and matter eras following 
inflation is not enough to give it back to the length it had before the inflationary epoch. 

While vh decreases during inflation all those perturbation scales that are fixed 
(in comoving coordinates) are effectively seen as 'exiting' the Hubble radius^. After 
inflation ends, during the radiation and matter eras, vh grows again and the Hubble 
radius stretches beyond some of these perturbation scales: they are 'entering' the 
Hubble radius. However, by that time these scales have already had time to get causally 
connected before exiting H^^ and, e.g., produce the same level of CMB anisotropies 
in the sky. 

With regards to the flatness problem, we have seen in §|L6| that the evolution 
equation satisfied by Q makes it evolve away from unity (cf. Eq. (|1.9|) ). However 
this happens for a standard perfect fluid equation of state, with p/p = w > —1/3 
(thus, satisfying the strong energy condition). During inflation the expansion is in 
general driven by a slow-rolling scalar field, and in this case the previous condition on 
w does not apply. As we can see from Eqs. ( |2.12| ), during a scalar field dominated 
universe we have w ~ — 1 and for this value of w equation ( p..9|) tells us that Q will 
rapidly approach unity. Provided sufficient inflation is achieved (necessary for solving 
the horizon problem, for example), — 1| at the end of inflation is dynamically driven 
to a value small enough to render Q ~ 0.2 today, as it is indeed estimated from, e.g., 
rotation-curve measurements in spiral galaxies, and other dynamical determinations 

of fi. g 

Finally, inflation solves the unwanted relics problem, the essential reason being the 
following: since the patch of the universe that we observe today was once inside a 
causally connected region, presumably of the order of the correlation length ^ (see 
§ |2.2.2| below), the Higgs field could have been aligned throughout the patch (as this 

^One many times finds statements in the literature about perturbations exiting and entering the 
horizon. Clearly, once the particle horizon encompasses a given distance scale, there is causal com- 
munication and homogenisation on that scale. Therefore one hardly sees how this scale can ever 'exit' 
the horizon and loose the causal contact it has achieved before. In order to avoid many headaches 
one should speak in terms of the Hubble radius, H~^, quantity that remains nearly constant during 
an inflationary phase (however, I am aware that this thesis might be plagued with sentences where 
the word horizon (in place of the more correct 'Hubble radius') is employed; my apologies) [see Ellis 
& Rothman, 1993 for a pedagogical review]. 

•^One comment is in order: many times one finds statements saying that inflation solves the flatness 
problem too well, thus predicting a ft value too close to 1 today. Recently, however, there has been 
a growing interest in models of inflation with ^ 1 [see, e.g., Amendola, Baccigalupi, & Occhionero, 
1995 and references therein]. 
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is, in fact, the lowest energy configuration). Thus, being no domains with conflicting 
Higgs orientations, the production of defects is grossly suppressed, and we expect of 
order 1 or less topological defects in our present Hubble radius. In other words, the 
huge expansion produced by an early inflationary era in the history of the universe 
dilutes the abundance of (the otherwise overproduced) magnetic monopoles or any 
other cosmological 'pollutant' relief. 

2.1.4 Generation of density perturbations 

One of the most important problems of cosmology is the problem of the origin of the 
primordial density inhomogeneities that gave rise to the clumpy universe where we live. 
This is closely connected with the issue of initial conditions. Before the advent of the 
inflationary scenario there was virtually no idea of which processes were the responsible 
for the formation of the large-scale structure we see today. It was possible that galaxies 
were originated by vortex perturbations of the metric [Ozernoi & Chernin, 1967], from 
sudden events like the explosion of stars [Ostriker & Cowie, 1980], or even from the 
formation of black holes [Carr & Rees, 1984]. 

In the inflationary cosmology this situation changed. First of all, it was understood 
that all those perturbations present before inflation were rendered irrelevant for galaxy 
formation, since inflation washes out all initial inhomogeneities. Further, a finite pe- 
riod of accelerated expansion of the universe naturally explains that perturbations on 
cosmologically interesting scales had their origin inside the Hubble radius at some point 
in the inflationary phase. 

The analysis of the linear evolution of density perturbations is normally performed 
by means of a Fourier expansion in modes. Depending on the physical wavelength of 
the modes, relative to the Hubble radius, the evolution splits into two qualitatively 
different regimes: for perturbations of size smaller than H~^, microphysical processes 
(such as pressure support, quantum mechanical effects, etc) can act and alter their 
evolution. On the contrary, when the typical size of the perturbation is beyond the 
Hubble radius, their amplitude remains essentially constant (due to the large friction 
term 3H(f) in the equation of motion) and the only effect upon them is a conformal 

^This also implies, however, that any primordial baryon number density will be diluted away 
and therefore a sufficient temperature from reheating as well as baryon-number and CP-violating 
interactions will be required after inflation [Kolb & Turner, 1990]. 
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stretching in their wavelength due to the expansion of the universe. 

Inflation has the means to produce scalar (density) and tensor (gravitational waves) 
perturbations on cosmologically interesting scales. The dynamics of producing density 
perturbations involves the quantum mechanical fluctuations of a scalar field in a nearly 
de Sitter space. Since the couplings of the inflaton are necessarily weak, e.g., in order 
not to generate too much CMB anisotropics, its quantum fluctuations can be estimated 
by the vacuum fluctuations of a free field in a de Sitter space [see Gibbons & Hawking, 
1977; Bunch & Davies, 1978]. In what follows we consider the Lagrangian 

in the background metric ds'^ = dt^ — exp(2i7t)[(ix^ + dy'^ + dz^]. To keep things 
simple we consider a massive free field, V^(0) = m?(fp' /2. From (|2.17| ) we may derive 



the equation of motion, and then quantise the system using the usual equal time 
commutation relations. After Fourier decomposing the field 0(x, t) we get the equation 
of motion for a particular mode as follows 

[dl + mdt + k'^ exp{-2H{t - Q) - m2]0fc(t) = , (2.18) 

where is the time at which inflation began. Upon defining z = kH~^ exp{—H{t — t^)) 
and i> = (^^ — m'^H^'^^ , equation (|2.18|) may be written 



This has the form of a Bessel equation and the solutions are 

Mt)ocz'/'Hl'''\z) , (2.20) 

where HI)^\ H'^'> stand for Hankel functions [Bunch & Davies, 1978]. The most general 
solution is a linear combination with coefficients Ci, C2 satisfying |cip — |c2p = 1. 
Different values for these constants lead to different vacuum states for the quantum 
theory. The mode 0fc(t) describes the evolution of a perturbation of wavelength ~ 
k~^exp{Ht). For very early times this wavelength is much smaller than and 
at such short distances de Sitter and Minkowski spaces are indistinguishable. The 
short wavelength limit corresponds to large z's and, given that in this limit we have 
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Hl^'^^ oc exp(—{+)ikAt), we see that the choice which corresponds to positive frequency 
modes in the flat space hmit is given by (ci, C2) = (1,0). Upon normahsation one gets 

Mt) = l\[^^M-lH{t - U))H^\z) . (2.21) 

The spectrum of fluctuations of the scalar field is ~ (|0^(t)|). Using the solution 
( p.21| ) and taking the limit of large time {z « 1) and iv? « (and thus we have 
z/ ~ I — \^{m/HY) we get the variance of the scalar field perturbation given by 

" exp -^^(t-t,)U d{\nk)(^\ . (2.22) 



The upper limit of the integral is fixed by the last (the smallest) wavelength that 
crosses the horizon at time t during the inflationary expansion. The lower limit takes 
into account that inflation starts at time t^, and therefore there is a first (maximum) 
wavelength that crosses the horizon when inflation began. From Eq. ( p.22| ) we easily 
see that the contribution to \5(f)\ per logarithmic interval of k is given by 



and, in the limit iv? « it reduces to well-known result V^^'^{k) ^ 11/271. A 
measurement of the (p^s will yield random phases, and the distribution of the moduli 
will have a dispersion ~ (l^^fcP) = H"^/ (2A;^). Accordingly, the spectrum of the inflaton 

_____ -j^ /2 

(cf. Eq. ( p^.23| )) a few times after the horizon exit is given by (k) ~ H(,x/2n 
[Vilenkin & Ford, 1982; Linde, 1982b; Starobinskii, 1982], where Hex stands for the 
(slowly varying) Hubble parameter at horizon exit. 

2.1.5 Evolution of fluctuations 

In this section we follow the recent review by Liddle & Lyth [1993]. The spectrum of 
the density contrast after matter domination may be written 

aH 



Vs=(^] T\k)5l{k) , (2.24) 



where T{k) is the linear transfer function which becomes equal to 1 on very large 
scales, meaning that very little evolution is suffered by those large scales that entered 
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the horizon recently, during the matter dominated eraQ. The quantity Sjj{k) specifies 
the initial spectrum of density perturbations at horizon entry (hence the subindex h) 
and is related to the spectrum of the initial curvature perturbations TZ by 

SUk) = ^Vn{k) . (2.25) 

5lj is exactly equal to the value of Vs on horizon entry on scales much larger than 100 
h-^ Mpc, roughly the size of the horizon at the time of equality between matter and 
radiation, and is approximately equal to it on smaller scales. 

The standard assumption is that d'jj oc corresponding to a density spectrum 

P oc with spectral index n. The standard choice n = 1 yields the 'fiat' spectrum, 
proposed by Peebles & Yu [1970], Harrison [1970] and Zel'dovich [1970] as being the 
only one leading to small perturbations on all scales^. 

The curvature perturbation TZ of comoving hypersurfaces is given in term of the 
perturbations of the infiaton field by [Lyth, 1985; Sasaki, 1986] 

7^ = ^S(f) . (2.26) 

The spectrum of 71 is then given by 

-P'^" = S-V'/ ■ (2.27) 



Now, given that the curvature perturbation TZ is constant after horizon exit [Liddle & 
Lyth, 1993], thus V^"^ remains constant even though H/4> and may vary sepa- 
rately. As long as the scale is far outside the horizon we have 

= JS^ , (2.28) 
27r(/)^„. 



^The correct transfer function fit comes only after numerically solving the relevant perturbation 
equations, with initial conditions specifying the relative abundance of the different matter components, 
and taking into account the free streaming of neutrinos, the diffusion of radiation on scales smaller 
than Silk's scale, etc. We will not go into the details here; these are reviewed in, e.g., [Efstathiou, 
1990]. 

® Compared with the standard n — 1 spectrum, tilted models with n < 1 lead to relatively less 
power on small scales and have been advocated many times, in particular to alleviate the too large 
pairwise galaxy velocity dispersion on scales of order 1 Mpc. On the other hand, recent analyses of 
the COBE-DMR data on very large angular scales tend to indicate a 'blue' spectra {n > 1). However 
this may be partially due to the (large-scale) 'tail' of the so-called Doppler peaks that, e.g., CDM 
plus inflation models predict to appear on the degree scale. See Chapter O. 



2.1. INFLATION 31 



where ex subindex stands for quantities evaluated at horizon exit (cf. the last paragraph 
of §2.1.4| ). Using Eqs. ( |2.25|) and ( |2.13| ) we find the amplitude of the adiabatic density 



perturbation at horizon crossing [Lyth, 1985] 



SUk) ^ I^Ci , (2.29, 



where e^x is the small parameter given by (|2.14| ) at horizon exit. One therefore has 



a perturbation amplitude on re-entry that is determined by the conditions just as it 
left the horizon, which is physically reasonable. In general one does not know the 
values of V^x and eex so one cannot use them to predict the spectral form of the 
perturbations; in fact, within generic infiationary models, these values depend on the 
particular length-scale of the perturbation. In an exactly exponential infiation these 
two parameters are constant (independent of the particular scale exiting the horizon 
during the infiationary era) and therefore the perturbation amplitude on horizon entry 
is constant and independent of the scale, leading to the 'fiat' spectrum we mentioned 
above. 

In general the spectral index of density perturbations depends on the considered 
scale. However, if this dependence is weak (at least within a cosmologically interesting 
range of scales) we may define n — 1 = dlnSfj/dllnk), which is implicitly based in 
the power-law dependence of Sjj. Differentiating equation ( p^.29D , and using the slow- 



rolling conditions on x and e, Liddle & Lyth [1992] derive the result 

n = 1 + 2xex - 6eex (2.30) 

where now the ex subindex refers, for definiteness, to the moment when the observable 
universe leaves the horizon. 

Analogously, for tensor perturbations (gravitational waves) we can define ut = 
d In Vgw / d{lii k) and, following a similar analysis, we getQ 

riT = -2eex . (2.31) 

We thus see that in order to have significant deviation from the fiat n = 1 and = 
spectra we necessarily need to violate the slow-rolling conditions e << 1 and |x| << 1. 



^More accurate results fo r n an d nr are given in [Stewart & Lyth, 1993] and [Kolb & Vadas, 1994]. 



We will use their results in §5.3.2 below in the framework of some interesting inflationary models. 
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Nevertheless, we will see below that within many of the currently used inflationary 
models, values for n below 1 can be achieved. Models leading to n > 1 have recently 
been studied by MoUerach, Matarrese & Lucchin [1994]. 

2.1.6 An overview of models 

We here briefly review a few inflationary models; the list is by no means exhaustive. 
We will say more about some of them, as well as introduce a couple of other popu- 
lar models, in Chapter ||, when we will study non-Gaussian features (e.g., the CMB 
skewness) predicted by primordial non-linearities in the evolution of the inflaton in the 
framework of these models. We want to emphasise that although these models solve 
the cosmological problems outlined in § p.l.3| there is as of 1995 no convincing way to 
realise the scenario [see Brandenberger, 1995; Turner, 1995 for recent reviews]. 

'Old' inflation 

The 'old' inflationary model [Guth, 1981; Guth &; Tye, 1980] is based on a scalar 
field theory undergoing a first-order phase transition [Kirzhnits & Linde, 1976]. In 
this model the universe has an initial expansion state with very high temperature, 
leading to a symmetry restoration in the early universe. The effective potential Vt(0) 
has a local (metastable) minimum at = even at a very low temperature. The 
potential also has a global (true) minimum at some other value, say, (p = a where the 
potential vanishes (in order to avoid a large cosmological constant at present time). In 
old inflation the crucial feature is a barrier in the potential separating the symmetric 
high-temperature minimum from the low-temperature true vacuum. As a result of 
the presence of the barrier the universe remains in a supercooled vacuum state = 
for a long time. The energy-momentum tensor of such a state rapidly becomes equal 
to Tn„ = gnuV{0) (all other form of matter-energy rapidly redshift) and the universe 
expands exponentially (inflates) until the moment when the false vacuum decays. As 
the temperature decreases with the expansion the scalar field configuration, trapped 
in the false vacuum, becomes metastable and eventually decays to the global minimum 
by quantum tunnelling through the barrier. This process leads to the nucleation of 
bubbles of true vacuum, and these bubbles expand at the velocity of light converting 
false vacuum to true. 
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Reheating of the universe occurs due to bubble-wall collisions. Sufficient inflation 
was never a real concern; the problem with this classical picture is in the termination 
of the false-vacuum phase, usually referred to as the 'graceful exit' problem. For 
successful inflation it is necessary to convert the vacuum energy to radiation. This is 
accomplished through the collision of the bubbles. 

Guth himself [Guth, 1981; see also Guth & Weinberg, 1983] realised that the sce- 
nario had a serious problem, in that the typical radius of a bubble today would be 
much smaller than our observable horizon. Thus the model predicted extremely large 
inhomogeneities inside the Hubble radius, in contradiction with the observed CMB 
radiation isotropy. The way out was thought to be in the percolation of the generated 
bubbles, which would then homogenise in a region larger than our present horizon. In 
these collisions the energy density tied up in the bubble-walls may be converted to 
entropy, and in order to have a graceful exit there must be many collisions. But, it 
so happens that the exponential expansion of the background space overwhelms the 
bubble growth (the volume inside the bubbles expands only with a power-law). This 
prevents percolation and the 'graceful exit' problem from being solved. 



'New' inflation 



Soon after the original (and unsuccessful) model laid down the main idea of the 
convenience of an early era of accelerated expansion, inflation was revived by the 
realisation that it was possible to have an inflationary scenario without recourse to 
a strongly ffist-order phase transition. Linde [1982a] and independently Albrecht & 
Steinhardt [1982] put forwards a modified scenario, the 'new' inflationary model. The 
starting point is a scalar field theory with a double well 'mexican hat' potential which 
undergoes a second-order phase transition. V{(f)) is symmetric and = is a local 
maximum of the zero temperature potential. As in old inflation, here also, finite 
temperature effects are the responsible for confining (p to lay near the maximum at 
temperatures greater than the critical one, Tc. For temperatures below Tc thermal 
fluctuations make = unstable and the field evolves towards one of the global 
minima, (p = icr? ruled by the equation of motion (cf. (|2 . 1 1|) ) 



+ 3H^ + dV/d(j) = 



(2.32) 



34 



Chapter 2. THEORIES OF THE EARLY UNIVERSE 



The transition proceeds now by spinoidal decomposition^ and hence t) will be 
homogeneous within one correlation length. 

If the potential near the maximum is flat enough the term can be neglected 
and the scalar field will undergo a period of slow-rolling. The field has both 'kinetic' 
and 'potential' energy; however, under the slow-roll hypothesis the velocity of the 
Higgs field will be slow and the potential energy will dominate, driving the accelerated 
expansion of the universe. The phase transition occurs gradually and if significant 
inflation takes place huge regions of homogeneous space will be produced; we would 
be living today deep inside one of these regions. 

There is no graceful exit problem in the new inflationary model. Since the spinoidal 
decomposition domains are established before the onset of inflation, any boundary walls 
will be inflated away outside our present Hubble radius. Thus, our observable universe 
should contain less than one topological defect produced in the transition. This is good 
news for the gauge monopole problem, but also bad news for global defects and local 
cosmic strings. 

Reheating occurs in this model as soon as the slow-roll evolution of the breaks 
down. The scalar field will begin its accelerated motion towards one of the absolute 
minima of the potential. In so doing it overshoots and then starts oscillating about 
it. The amplitude of this oscillation is damped by the expansion of the universe and 
predominantly by the couplings that has with other fields. These couplings lead 
to the decay of into lighter relativistic particles that after thermalisation reheat 
the universe. However, constraints from the level of density perturbations and CMB 
anisotropics generated by this model force the couplings between and other fields 
to be extremely small. This is not only viewed as a 'fine tuning' problem difficult to 
explain by any underlying particle theory, but also leads to a low reheat temperature, 
which may make baryogenesis problematical^ This is why the scalar field responsible 
for the spontaneous symmetry breaking in the minimal SU{5) grand unified theory 
with a Coleman-Weinberg potential was discarded as a possible candidate for the 
inflaton [Bardeen, Steinhardt & Turner 1983; Brandenberger & Kahn, 1984]. 



^We will be more precise about this concept in §§2.2.1-2.2.2 below [see also Mermin, 1979]. 
^Recently there has been an interest in reconsidering reheating scenarios. See e.g. Refs. |l|, pl2| , 



109| , 186, Needless to say much of the power of inflation relays in being able to find a way to 



re-establish the conditions necessary for a smooth matching with the standard hot big bang model, 
which is so successful in its predictions after, say, 1 sec (nucleosynthesis) in the history of the universe. 
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Chaotic inflation with a polynomial potential 

A model that successfully implements inflation and that can occur for very general 
types of scalar field potentials was proposed in [Linde, 1983a]. Here the scalar field 
is not part of any grand unified theory and its presence is just to drive an era of 
infiation. No specific thermal equilibrium state previous to infiation is required, unlike 
what happened in 'new' inflation; in fact this was considered as a weak point, since 
very small coupling constants were required in order that the generated density and 
CMB perturbations be in agreement with observations, and thus the achievement of 
thermal equilibrium which would drive the inflaton to the minimum of its potential 
was hard to justify as a generic initial condition. 

Lindc's idea solves this problem by envisioning that different initial values for the 
inflaton fleld will be attained in different regions of the universe (these initial values 
being 'chaotically' distributed). In those Hubble volumes where the conditions for 
an adequate period of inffation are satisfled the problems of the standard big bang 
cosmology might be solved, and they yield a universe like ours. This view leads to a 
very inhomogeneous universe, albeit on scales much larger than our present horizon. 

Let us now explain the basic features of the model and see how the fleld (j)[x, t) could 
be distributed in the early universe. We know that the upper limit for an accurate 
classical description of spacetime is given by the Planck scale (times near to tp ~ 10^^^ 
sec). At that epoch the effective potential is given with an accuracy of 0{m%), due 
to the uncertainty principle. Thus a typical distribution of the fleld in a classical 
universe, which emerges from the spacetime 'foam' at that time, is given by [Linde, 



Consider, for simplicity, the V{(f)) = A0''/4 model, with A << 1. Let us focus on 
those domains of the universe where the fleld was initially homogeneous (on a scale 
larger than H^^) in the sense that {di,(f))^ < V{(f)), and sufficiently large, > mp. 



p ~ V{(f)) and the Hubble parameter H ~ {8nV{(f))/3)^ /mp ~ A^/^0^/mp is slowly 
varying. The equation of motion can be expressed as 

+ v^eTrA— = -A0^ . (2.33) 
nip 

Assuming slow-rolling (condition valid for 0^ >> mp/(67r)) we get 0(t) = 
(j)inexp {—J -^mp t), with 0j„ the initial value. Furthermore, this domain expands 



1990] V{(j)) < m% and ((9^0)^ < m%. 



The equation of motion in such a domain is + "iHcf) 



\(f)^. The energy density is 
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with a scale factor 



a(t) = ttin exp 



H{t)dt 



TT 



— {ct>l-<P\t)) 



m. 



(2.34) 



This is a quasi-exponential expansion for (p » mp and by the time the slow-rolling 
condition breaks down the domain will have expanded by a factor ~ exp (vr^^^/mp). 
For ^ A Amp the universe expands more than times, enough for solving the 
horizon and flatness problems[3. The reheating picture is similar to the previous models 
of inflation: for < mp/3 the scalar field oscillates about the minimum of the potential, 
and the vacuum energy stored is converted into radiation. However, there is still 'fine 
tuning' in the sense that, e.g., an adequate level of CMB anisotropics constrains A to 
take a very small value. 

inflation 

We have seen that within chaotic models of inflation the scalar field associated 
to the inflaton had not a direct connection with a grand unified theory. Its presence 
was justified just by the 'need' to have an accelerated era in the early universe. Here 
instead we will see how an inflaton field naturally arises in the case when we minimally 
extend General Relativity. We will concentrate in an inflationary model motivated by 
a modification of Einstein gravity. The Lagrangian for gravity is|^ 

C = -{16nG)-^[R + i?V(6M2)] . (2.35) 

The second (oc R^) part is what modifies the usual Einstein-Hilbert first term. The 



justification of the ansatz (|2.35|) can be found in [Cardoso & Ovrut, 1993] in the 
framework of supergravity theories. After a conformal transformation [Whitt, 1984] 
the R^ term may be eliminated and we are left with the standard-gravity Lagrangian, 
but interacting with an additional scalar field (identified with R) given by 

= (3/(167rG))^/^ In (l + R/{3M^)) , (2.36) 



where the new scalar field must satisfy the potential [see, e.g., Liddle & Lyth, 1993 
for a review, and Salopek, Bond & Bardeen, 1989 for an analysis of the fluctuation 

^°Note that (pm ^ A.Amp does not contradict ¥{(()) < mp, provided we take A small enough, in 
agreement with many reasonable theories. 

^^The first model of inflation, studied by Starobinskii [1980], was of this form. 



2.1. INFLATION 37 



spectrum] 

Vi<f>) = TJTT^ - (-(167rG/3)^/V)l . (2.37) 



327rG 

In the regime (p > mp this potential produces the slow-rollover evolution of the scalar 
field, and thus infiation is achieved. Approximately when the field reaches the value 
<pend ~ (SttG)^^/^ iufiatiou ends and the field oscillates. For a scale of the order of 
our present horizon the small parameters of Eqs. (|2.14| ) and (|2.15| ) can be calculated 
yielding e ~ 10~^ and x ~ —0.02, and thus a spectral index of density perturbations 
n ~ 0.96. The tensor mode contribution to the anisotropics is negligible in this model. 

Power— law infiation with an exponential potential 

Exponential infiation is by no means a necessary condition to solve the problems of 
the standard big bang cosmology. In principle any type of accelerated expansion may 
solve, say, the horizon and fiatness problems, and we will here consider the case in which 
the scale factor is given by a{t) oc t^, with p > 1 [Abbott & Wise, 1984b]. Lucchin & 
Matarrese [1985] were the first to study the cosmological-background equations with 
a scale factor of this form, and found exact solutions yielding an infiaton potential of 
the form (cf. § p.l.2| ) 



V{<j)) = V,exp[-2J ^ . (2.38) 

\ \ p nip ) 

This potential may also find a motivation in the context of extended infiation [see, 
e.g.. La & Steinhardt, 1989]. Furthermore, the featureless exponential potential does 
not provide a way to end the infiationary stage, say, by oscillations of the field, and 
so suitable modifications or even bubble nucleation mechanisms (e.g., in extended 
infiation) should be invoked. Nevertheless, the slow-roll conditions are satisfied if 
p >> 1 and we get e = x/2 = 1/p independent of the scale. The spectral indices are 
nx — — 2/p and r;,~l — 2/p. The tensor-to-scalar ratio, R, of the contributions to 
the variance of the CMB anisotropics is i? ~ 12/p. 

It is beyond the scope of the present thesis to give a detailed account of all the 
existent models that have been proposed in the literature in the recent past. A few 



more of them we will discuss in 3. 1-5.3.2 in connection with the CMB skewness 



^^Exact analytic solutions for both the inflaton dynamics and for the density perturbations exist 
in this model. The exact spectral indices are ut = — 2/(p— 1) and n = 1 — 2/(p— 1), in very good 
agreement with the slow-roll results [see, e.g., Lyth & Stewart, 1992]. 
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they predict. Some of these models have been strongly constrained by observations 
and only by contorting somewhat unnaturally can they still be alive. However, 'many' 
others resist and this is seen by many authors as a weakness of the whole idea. It 
would be desirable to clean up a bit the 'closet of models' and leave there just the best 
candidates we have. 

2.1.7 Stochastic approach to inflation 

In § |2.1.2| we have studied the dynamics of the inflaton field assuming it homogeneous 
and obeying a classical Klein-Gordon equation. Quantum fluctuations in the field are 
stretched to super-horizon scales and behave very much like classical perturbations^. 
However, the quantum nature of the scalar field leads to very interesting physics as 
well, and the stochastic approach is the natural framework to study it. The basic idea 
is due to Vilenkin [1983] and Starobinskii [1986]. The approach consists in splitting 
the scalar field in momentum space into long and short wavelength modes. 

Starting from the Heisenberg operator equation of motion for the scalar field, the 
evolution of the long wavelength part satisfies a classical, but stochastic, equation of 
motion. The quantum effects, in the form of short wavelength modes, build up a noise 
term, as we will explain below. The field satisfies Vi^V^^ + dV/d(j) = 0. The scalar 
field may be written = 0iong + 0short, or more explicitly 

<P = 0iong(x, t) + J Q{k~ eaH) (akMY-'^^'' + 4€(^)e*'"') , (2-39) 

where (^iong(^; ^) contains only modes such that k « aH; aj,, are the usual creation 
and annihilation operators, e is a constant much smaller than 1, and G stands for the 
step function. Usually inflation takes place in regions where the scalar field potential 
is not very steep and thus the short wavelength part satisfies the massless equation 

^•^Classicality of long wavelength quantum fluctuations in de Sitter space has been established in 
[Guth & Pi, 1985] in a mode-by-mode treatment. This is the standard picture but, although it has 
been widely used, it is not widely understood; see [Mijic, 1994]. The origin of classical perturbations 
as a consequence of quantum fluctuations is a subtle subject. The quantum to classical transition 
may be studied in the context of quantum cosmology, by means of the Wheeler-DeWitt equation for 
the wave function of the universe in a minisuperspace. For the achievement of the classical behaviour, 
the system under study has to show not only a strong correlation between coordinates and momenta 
(i.e., behave in a classical way), but also all signs of quantum interference should disappear. This 
last leads to the notion of decoherence of the density matrix [Zurek, 1982; see also Halliwell, 1989; 
Padmanabhan, 1989; Gangui et al., 1991; Castagnino et al., 1993, and references therein]. 
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VjyV^^short = 0. Thus the short wavelength modes can be taken as 



From this Starobinskii [1986] derives the equation of motion for the long wavelength 
'coarse-grained' part which, in the slow-rolling approximation, reads 

1 dV 

01ong(x,t) = + fiS,t) , (2.41) 



where the spatial gradient of 0iong has been neglected since it is subdominant for modes 
k « aH (and hence the evolution of the coarse-grained field can be followed in each 
domain independently), and where 

/(^' ^) = TTW^ I "^'^ ~ {cik<Pke-''^-' + ai<Ple'^-^) , (2.42) 



(27r)3/2 

where ks = eaH stands for the inverse of the coarse-grained domain radius. The 
correlation function for /, built up by all high frequency modes, may be computed 
from this last equation and is given by 

{f{x,t)f{x,t')) = ^6{t-t') , (2.43) 

and has a white noise spectrum in time. 

According to this picture the universe is described by the value of the field in dif- 
ferent coarse-grained domains of comoving size kj^ ^ {eaH)~^. With the inflationary 
expansion (whereby the comoving radius of these domains decrease) an initial domain 
of radius {eaiHi)~^ gets divided into 0{e^) subdomains of size ~ (£02-^2)"^, after 
one Hubble time-step such that ^2 — ^1 + (since the temporal rate of change of 
is smaller than H). The magnitude of the coarse-grained field in these smaller 
domains in determined by both the classical (convective) and the noise (stochastic) 
parts appearing in the right-hand side of Eq. ( p.41j ) .p^ The classical force 'pushes' 
the field down the potential in all the subdomains, whereas the stochastic force acts 
with different strengths]^ and arbitrary direction in different subdomains. This process 
repeats itself in all subsequent Hubble time-steps. 

^^The reader may have noted that we have already dropped the label long of the coarse-grained 
field. 

^^The typical magnitud e of t he change produced in the field in one Hubble time is given by the 
variance ~ H/2t:, cf. Eq. ( |2.43D . 
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The result of this process is that in the new subdomains the field takes on different 
(uncorrelated) values: in those where the classical 'drag' overwhelms the stochastic 
'push' the field goes down the potential and infiation will eventually end, yielding a 
domain like our present universe. The fact that our present causal region was, at the 
time of infiation, divided into many coarse-grained subdomains with slightly different 
values for explains the existence of energy-density perturbations (and eventually also 
the structure originated from them). However, in those domains where the opposite 
happens (namely, where {3H'^)''^dV/d(f) < H/2'k) there will be regions where the 
stochastic force is greater than the classical one and points upwards the potential. 
These domains will keep on inflating and, as they expand much faster than their non- 
inflating counterparts, they flU the majority of the universe (even if at each time more 
and more regions end their accelerated expansion). This picture leads to a 'stationary' 
universe, where there could have been no beginning and would be no end, with some 
domains inflating and others ending inflation all the time, forever: Linde's eternal 
inflation picture [Linde, 1990]. 

At this point one could attempt to study the probability distribution P(0, t) for 
the coarse-grained fleld by means of the Fokker-Planck equation. By doing this we 
are concentrating on one particular domain, say with comoving coordinate Let us 
consider, for deflniteness, chaotic power-law inflation with an exponential potential. 
By suitably redeflning a new coarse-grained fleld 0' = 0'(0 — (pel)) as a non-linear 
function of the difference between the old and its classical value[^, the analysis 
allows one to derive the dispersion of 0' around its mean value (0'). This dispersion 
depends on the scale in question, namely, on the time a particular scale left the horizon 
during the inflationary expansion. It turns out that for scales smaller than our present 
horizon (relevant for structure formation) the actual value of 0' is very close to (0'), 
and consequently the distribution of the fleld is a narrow Gaussian to a very good 
approximation and is peaked at the classical trajectory 0^/. Instead, if we look at 
scales much larger than our horizon, the dispersion of 0', non-linearly related to the 
coarse-grained fleld 0, can be of the same order as its mean value. Therefore, the 
linearisation of this non-linear 0' = 0'(0 — 0c/) relation does not apply anymore and, 
as a consequence of this, a highly non-Gaussian distribution for is predicted on very 



""^^see e.g., [Matarrese et al, 1989] for a detailed study of the probability distribution function P in 
the case of many relevant chaotic models of inflation. 



17. 



Pel 



is the solution of Eq. ( 2.41 ) when the noise term is 'switched off' 
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large scales [see MoUeracli et al, 1991 for details]. 

In the previous paragraph we considered correlations of the coarse-grained field 
at different times but at the same point. Now we will concentrate also on spatial 
correlations which could be derived from the Langevin-type equation ( [2.41 ). Because 
of the time dependence of the Hubble parameter it has been proposed that a more 
fundamental time variable is given by a oc ln(a) = / H{t)dt [Starobinskii, 1986; see 
also Chapter ^. From the definition of the noise (cf. Eq. (|2.42|) ) we can compute 

1 r°° 

(/(fi, ai)/(x2, aa)) = tt^ dkk^ 5 {k - ks{xi, ai)) S {k - ks{x2, 02)) 

ZTT^ JO 

dk dk 

x-7^(xi,ai)— ^(f2,a2)0fe(ai)0fc(«2)io(fc|^i - ^2!) , (2.44) 
da da 

which yields the correlation in the noise for spatially separated points. Now we may 
compute the spatial correlation for the fluctuations in the field, defined as 6(f){x, a) = 
0(x, a) — 0d(a)- Given that (^d is solution of equation ( p.41| ) (now written down in 
terms of the new time variable a) without noise term, we have 



6(f){x, a) = I f{x, a')da' , (2-45) 

and from this 

(fi,a)5(/)(x2,a')) = / dai I rfa2(/(fi, ai)/(x2, 02)) . (2.46) 



Using ( |2.44| ), ( |2.46| ) may be cast as 
{64>{xi,a)6(f){x2,a)) = — / dkk (l)k{(yk)(pl{ak)jo{k\xi-X2\) , (2.47) 



27r2 Jks{o) 

where a^ corresponds to the time when k = eaH. From Eq. ( p2.40[ ), in the case of short 
wavelength modes, we have (j)k{oik) — iH {(j){ak)) exp(2e)/(v/2p/^), and Eq. ( |2.47] ) may 
be cast as 

]^ i"fcs(min(a,o')) dk 

(50(fi, a)(50(f2, a')) = TTT-yi / -r^(<^(^i' afc))^(0(^2, afc))jo(fc|^i - ^2!) • 

(2.48) 

This is a function of the value the Hubble parameter takes on when the scale k~^ is 
equal to the coarse-grained radius at each of the points. 

We saw before that if we want to concentrate on scales of the order of, and smaller 
than, our present horizon the distribution of values of the coarse-grained field turns 
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out to be highly peaked around the classical trajectory 0c«- Then, for these scales it is a 
good approximation to take for H the value it would have were the field just described 
by (pel- To be specific, let us consider again the case of an exponential potential V{(j)) oc 
exp(— AVSttG^), where = 2/p (cf. Eq. ( p^.lUD ). In this case the classical solution 
of Eq. (|2.41|) is 4>d{c() = 4>o — A(87rG')^^/^a, and so we have H{(f)ci{a)) oc exp(— A^q;/2) 
and also k = eao-f^o exp[(l — A^/2)afc], where Hq = H{a = 0), etc. We finally get 
[MoUerach et al, 1991] 

(2.49) 

We see from this that in the case p » 1, namely when A — > (when power-law 
inflation approaches de Sitter inflation) we recover the scale invariant spectrum of 
fluctuations with amplitude Hq/2tt. However, for an arbitrary value of A in the range 
< A < -\/2 we have a spectral index of density fluctuations n = 1 — 2A^/ (2 — A^) < 1. 
This spectral index 'tilts' the spectrum, with the net effect of transferring more power 
on large scales, as is the benchmark of power-law inflation [Lucchin & Matarrese, 
1985]. 

We will come to use extensively the stochastic approach to inflation reviewed in 
this section again in Chapter |^. Then we will perform a second-order perturbative 
expansion around the classical solution to the Langevin-type equation for the coarse- 
grained scalar field, the ultimate aim being the computation of higher-order correlation 
functions for the perturbations in the inflaton. Upon horizon re-entry of the relevant 
perturbation scales, these correlations will be passed on to the peculiar gravitational 
potential and, via the Sachs- Wolfe effect, to the CMB temperature anisotropics. 



2.2 Topological defects 

A central concept of particle physics theories attempting to unify all the fundamental 
interactions is the concept of symmetry breaking. As the universe expanded and cooled 
down, first the gravitational interaction, and subsequently all other known forces: the 
strong, the weak and the electromagnetic force, would have begun adopting their own 
identities. In the context of the standard hot big bang theory the spontaneous breaking 
of fundamental symmetries is realised as a phase transition in the early universe. Such 
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phase transitions have several exciting cosmological consequences and thus provide an 
important hnk between particle physics and cosmology. 

There are several symmetries which are expected to break down in the course of 
time. In each of these transitions the space-time gets 'oriented' by the presence of 
a hypothetical force field called the 'Higgs field', pervading all the space. This field 
orientation signals the transition from a state of higher symmetry to a final state where 
the system under consideration obeys a smaller group of symmetry rules. As an every- 
day analogy we may consider the transition from liquid water to ice; the formation of 
the crystal structure ice (where water molecules are arranged in a well defined lattice), 
breaks the symmetry possessed when the system was in the higher temperature liquid 
phase, when every direction in the system was equivalent. In the same way, it is 
precisely the orientation in the Higgs field which breaks the highly symmetric state 
between particles and forces. 

Kibble [1976] made use of a model in which the phase transition proceeds by the 
formation of uncorrelated domains that subsequently coalesce, leaving behind relics in 
the form of defects. Such relic 'flaws' are unique examples of incredible amounts of 
energy and this feature attracted the minds of many cosmologists. In the expanding 
universe, widely separated regions in space have not had enough time to 'communicate' 
amongst themselves and are therefore not correlated, due to a lack of causal contact. 
It is therefore natural to suppose that different regions ended up having arbitrary 
orientations of the Higgs field and that, when they merged together, it was hard for 
domains with very different preferred directions to adjust themselves and fit smoothly. 

Most interestingly, unlike other proposed mechanism for the generation of cosmolog- 
ical observable features, topological defects can be reproduced in the laboratory! This 
is possible thanks to 'scaling' properties with which these defects are endowed. In fact, 
when all relevant lengths are uniformly scaled down, experimentalist have within their 
reach a physically equivalent situation as that in the early universe but with all the ad- 
vantages of a manageable laboratory experiment. In 1985, Zurek [1985] proposed using 
the transition from normal to superfluid ^He to test the Kibble mechanism. His idea 
was that after a rapid quench, defects would form throughout the system. By studying 
their formation, interesting hints for cosmology could be discovered]^. Initially the su- 

^^Note however that the analogy is not complete; there are some differences between the two classes 
of systems. In condensed matter systems the dynamics is in general non-relativistic and friction 
dominated, whereas in flat spacetime defect motions are highly relativistic, even in the case of being 
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perfluid transition of helium turned out to be hard to deal with; temperatures of order 
~2K were necessary and very extreme laboratory conditions were required. However, 
people managed to perform similar experimental tests by resorting to particular or- 
ganic compounds called 'hquid crystals' [de Gennes, 1974]. These substances undergo 
transitions at temperatures ranging from 10° to 200° centigrade and generate structure 
easily detectable with the naked eye or with a microscope [Chuang et al, 1991; Bowick 
et aJ., 1994]. After an initial transient, the created network of defects evolves in a 
self-similar manner and, although dynamically changing sizes, always looks the same 
(in a statistical sense). 

Recently McClintock and co-workers [Hendry et al, 1994; see also Yurke, 1994] have 
succeeded in carrying out Zurek's original idea. By making a fast adiabatic expansion 
through the critical density in "^Hc they were able to observe copious production of 
quantised vortices (the condensed-matter analogue of cosmic strings). These different 
tests (both with liquid crystals and helium) provide a kind of 'experimental confir- 
mation' of cosmological topological defect theory, increasing the credibility in these 
ideas. 

These 'non-conventional' theoretical models for generating large-scale structure 
have been very well studied and highly detailed simulations have been performed to 
test and confront their predictions against a large bulk of observational material [see 
Vilenkin & Shellard, 1994 for an updated review]. Not only can they account for a host 
of astrophysical data, but they also do it in a way so transparent and economical that 
they seriously challenge the so far 'standard' inflationary model. One severe testing 
ground for these models is given by the cosmic microwave background anisotropy test. 
According to the big bang theory, the extremely high energies reigning during the flrst 
instants of time in the early universe favoured the close interaction between matter and 
radiation. Perturbations in this primordial plasma remained imprinted in the radiation 
that decoupled from matter around a hundred thousand years after the bang. These 
signatures in the relic radiation are of key importance in trying to unveil the actual 
mechanism behind the origin of the large-scale structure surrounding us today. Any 
proposed 'seed' for structure formation should not produce too much anisotropy in the 
CMB radiation in order to be consistent with current experiments. This turns out to 
be one of the most stringent tests on the models. 

in presence of a surrounding plasma that may eventually damp somewhat its dynamics. 
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Cosmic strings can account for the formation of large-scale filaments and sheets, 
observed galaxy formation rates and galactic magnetic fields. They also generate pe- 
cuhar velocities on large scales, and are consistent with the statistical properties of the 
CMB anisotropies measured by COBE-DMR on large angular scales. Texture simula- 
tions are able to reproduce galaxy-galaxy correlation functions, get clusters of galaxies 
significantly clustered on the correct scales and reproduce large-scale streaming mo- 
tions as observed. Moreover, this scenario predicts the existence of density peaks that 
are much more numerous than expected from the standard infiationary mechanism of 
generation of structure. This in turn leads to earlier galaxy formation and could help 
in explaining the existence of very old objects, like quasars. 

2.2.1 Phase transitions and finite temperature field theory 

Phase transitions are known to occur in the early universe. Examples are the quark 
to hadron (confinement) transition, which QCD predicts at an energy around 1 GeV, 
and the electroweak phase transition at about 250 GcV. Within grand unified theories, 
aiming to describe the physics beyond the standard model, other phase transitions are 
predicted to occur at energies of order 10^^ GeV; during these, the Higgs field tends to 
fall towards the minima of its potential while the overall temperature of the universe 
decreases as a consequence of the expansion. 

A familiar theory to make a bit more quantitative the above considerations is the 
A|0|** theory, 

^-lm' + lml\<f>\'-^\cf>\\ (2.50) 

with ml > 0. The second and third terms on the right hand side yield the usual 
'mexican hat' potential for the complex scalar field. For energies much larger than 
the critical temperature, Tc, the fields are in the so-called 'false' vacuum: a highly 
symmetric state characterised by a vacuum expectation value (|0|) = 0. But when 
energies decrease the symmetry is spontaneously broken: a new 'true' vacuum develops 
and the scalar field rolls down the potential and sits onto one of the degenerate new 
minima. In this situation the vacuum expectation value becomes = Guiq/X. 

Research done in the 1970's in finite-temperature field theory [Weinberg, 1974; 
Dolan & Jackiw, 1974; Kirzhnits & Linde, 1974] has led to the result that the 
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temperature-dependent effective potentiai can be written down as 

VV(|0|) = -lm2(T)|0p + A|0|4 (2.51) 

witli T2 = 24m2/A, m2(T) = ml{l - T^T^), and {\(j)\f = Qm'^{T)/X. We easily see 
tliat wlien T approaclies from below the symmetry is restored, and again we fiave 
(101) = 0. In condensed-matter jargon, tlie transition described above is second-order, 
or also first-order proceeding by spinoidal decomposition (caused by a rapid quench 
in the system) [Mermin, 1979] .[^ 

2.2.2 The Kibble mechanism 

The model described in the last subsection is an example in which the transition may 
be second-order. As we saw, for temperatures much larger than the critical one the 
vacuum expectation value of the scalar field vanishes at all points of space, whereas 
for T < Tc it evolves smoothly in time towards a non vanishing (|0|). Both thermal 
and quantum fluctuations influence the new value taken by (|0|) and therefore it has 
no reasons to be uniform in space. This leads to the existence of domains wherein the 
(|0(x)|) is coherent and regions where it is not. The consequences of this fact are the 
subject of this section. 

Phase transitions can also be first-order proceeding via bubble nucleation. At very 
high energies the symmetry breaking potential has (|0|) = as the only vacuum state. 
When the temperature goes down to a set of vacua, degenerate to the previous one, 
develops. However this time the transition is not smooth as before, for a potential 
barrier separates the old (false) and the new (true) vacua. Provided the barrier at 
this small temperature is high enough, compared to the thermal energy present in the 
system, the field will remain trapped in the false vacuum state even for small (< Tc) 
temperatures. Classically, this is the complete picture. However, quantum tunnelling 
effects can liberate the field from the old vacuum state, at least in some regions of 
space: there is a probability per unit time and volume in space that at a point x 
a bubble of true vacuum will nucleate. The result is thus the formation of bubbles 

^^In a first-order phase transition the order parameter (e.g., (|<^|) in our case) is not continuous. It 
may proceed by bubble nucleation [Callan & Coleman, 1977; Linde, 1983b] or by spinoidal decompo- 
sition [Langer, 1992]. Phase transitions can also be continuous second-order processes. The 'order' 
depends sensitively on the ratio of the coupling constants appearing in the Lagrangian. 
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of true vacuum with the value of the field in each bubble being independent of the 
value of the field in all other bubbles. This leads again to the formation of domains 
where the fields are correlated, whereas no correlation exits between fields belonging 
to different domains. Then, after creation the bubble will expand at the speed of light 
surrounded by a 'sea' of false vacuum domains. As opposed to second-order phase 
transitions, here the nucleation process is extremely inhomogeneous and is not 

a continuous function of time. 

Let us turn now to the study of correlation lengths and their role in the formation 
of topological defects. One important feature in determining the size of the domains 
where is coherent is given by the spatial correlation of the field (p. Simple 

field theoretic considerations [see, e.g., Ref. |^] for long wavelength fluctuations of 
lead to different functional behaviours for the correlation function G{r) = (0(ri)0(r2)), 
where we noted r = |ri — r2|. What is found depends radically on the whether the 
wanted correlation is computed between points in space separated by a distance r much 
smaller or much larger than a characteristic length = m{T) ~ known as 

the correlation length. We have 

^exp(-|) r»e 
G{r) ~ { (2.52) 



•2 



h r«i. 

This tells us that domains of size ^ ~ arise where the field is correlated. On 
the other hand, well beyond ^ no correlations exist and thus points separated apart by 
r >> ^ will belong to domains with in principle arbitrarily different orientations of the 
Higgs field. This in turn leads, after the merging of these domains in a cosmological 
setting, to the existence of defects, where field configurations fail to match smoothly. 

However, when T ^ Tc we have m ^ and so ^ oo, suggesting perhaps 
that for all points of space the field becomes correlated. This fact clearly violates 
causality. The existence of particle horizons in cosmological models constrains 
microphysical interactions over distances beyond this causal domain. Therefore we get 
an upper bound to the correlation length as ^ < ~ t. 

The general feature of the existence of uncorrelated domains has become known as 
the Kibble mechanism [Kibble, 1976] and it seems to be generic to most types of phase 
transitions. 
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2.2.3 A survey of topological defects 

Different models for the Higgs field lead to the formation of a whole variety of topo- 
logical defects, with very different characteristics and dimensions. Some of the pro- 
posed theories have symmetry breaking patterns leading to the formation of 'domain 
walls' (mirror reflection discrete symmetry): incredibly thin planar surfaces trapping 
enormous concentrations of mass-energy which separate domains of conflicting field 
orientations, similar to two-dimensional sheet-like structures found in ferromagnets. 
Within other theories, cosmological fields get distributed in such a way that the old 
(symmetric) phase gets confined into a finite region of space surrounded completely 
by the new (non-symmetric) phase. This situation leads to the generation of defects 
with linear geometry called 'cosmic strings'. Theoretical reasons suggest these strings 
(vortex lines) do not have any loose ends in order that the two phases not get mixed 
up. This leaves infinite strings and closed loops as the only possible alternatives for 
these defects to manifest themselves in the early universe^. 

With a bit more abstraction scientists have even conceived other (semi) topologi- 
cal defects, called 'textures'. These are conceptually simple objects, yet, it is not so 
easy to imagine them for they are just global field configurations living on a three- 
sphere vacuum manifold (the minima of the effective potential energy), whose non 
linear evolution perturbs spacetime. Turok [1989] was the first to realise that many 
unified theories predicted the existence of peculiar Higgs field configurations known as 
(texture) knots, and that these could be of potential interest for cosmology. Several 
features make these defects interesting. In contrast to domain walls and cosmic strings, 
textures have no core and thus the energy is more evenly distributed over space. Sec- 
ondly, they are unstable to collapse and it is precisely this last feature which makes 
these objects cosmologically relevant, for this instability makes texture knots shrink 
to a microscopic size, unwind and radiate away all their energy. In so doing, they 
generate a gravitational field that perturbs the surrounding matter in a way which can 
seed structure formation. 

Let us now explore the conditions for the existence of topological defects. It is 

^°'Monopole' is another possible topological defect; we defer its discussion to the next subsection. 
Cosmic strings bounded by monopoles is yet another possibility in GUT phase transitions of the 
kind, e.g., G^KxC/(l)— ^K. The first transition yields monopoles carrying a magnetic charge 
of the U{1) gauge field, while in the second transition the magnetic field in squeezed into flux tubes 
connecting monopoles and antimonopolcs [Langacker & Pi, f980]. 
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widely accepted that the final goal of particle physics is to provide a unified gauge 
theory comprising strong, weak and electromagnetic interactions (and some day may 
be also gravitation, if we really want to go beyond the Planck scale) . This unified theory 
is to describe the physics at very high temperatures, when the age of the universe was 
slightly bigger than the Planck time. At this stage, the universe was in a state with 
the highest possible symmetry, described by a symmetry group G, and the Lagrangian 
modehng the system of all possible particles present should be invariant under the 
action of the elements of G. 

As wc explained before, the form of the finite temperature effective potential of 
the system is subject to variations during the cooling down evolution of the universe. 
This leads to a chain of phase transitions whereby some of the symmetries present 
in the beginning are not present anymore at lower temperatures. The first of these 
transitions may be described as G^H, where now H stands for the new (smaller) 
unbroken symmetry group ruling the system. This chain of symmetry breakdowns 
eventually ends up with SU{3) x SU{2) x f/(l), the symmetry group underlying the 
'standard model' of particle physics, as it should be. 

The smaller group H contains elements g of the original group that, when acting on 
the field 0, they leave it with its expectation value, roughly, g0 = 0. However, we are 
interested in the vacuum manifold, as it is this the one that tells us something about 
the state of minimum energy, where the fields most probably sit down. In topology 
theory the manifold of degenerate vacuum states is identified with the cosct space A4 = 
G/H. Field configurations that arc left invariant when acted upon by the elements of 
A4 are precisely those that minimise the free energy. 

The importance of the study of the vacuum manifold lies in the fact that it is pre- 
cisely the topology of what determines the type of defect that will arise. Homotopy 
theory tells us how to map Ai into physical space in a non-trivial way, and what en- 
suing defect will be produced. For instance, the existence of non contractable loops in 
Ai is the requisite for the formation of cosmic strings. In formal language this comes 
about whenever we have the first homotopy group 7ri(A1) ^ 1, where 1 corresponds to 
the trivial group. If the vacuum manifold is disconnected we then have 7io{M.) ^ 1, and 
domain walls are predicted to form in the boundary of these regions where the field (f) 
is away from the minimum of the potential. Analogously, if 7r2(A1) 7^ 1 it follows that 
the vacuum manifold contains non contractable two-spheres, and the ensuing defect is 
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a monopole. Textures arise when Ai contains non contractable three-spheres and in 
this case it is the third homotopy group, n-s^M), the one that is non triviaL 

2.2.4 Local and global monopoles and domain walls 

Generically topological defects will be produced if the conditions for their existence 
are met. Then for example if the unbroken group H contains a disconnected part, 
like an explicit U{1) factor (something that is quite common in many phase transition 
schemes discussed in the literature), monopoles will be left as relics of the transition. 
This is due to the fundamental theorem on the second homotopy group of coset spaces 
[Mermin, 1979], which states that for a simply-connected covering group G we have0 



with Ho being the component of the unbroken group connected to the identity. Then 

we see that since monopoles are associated with unshrinkable surfaces in G/H, the 

previous equation implies their existence if H is multiply-connected. The reader may 

guess what the consequences are for GUT phase transitions: in grand unified theories 

a semi-simple gauge group G is broken in several stages down to H = SU (3) x SU{1). 

Since in this case 7ri(H) = Z, the integers, we have 712 (G/H) 7^ 1 and therefore gauge 

monopole solutions exist [Preskill, 1979]. 

Monopoles are yet another example of stable topological defects. Their formation 

stems from the fact that the vacuum expectation value of the symmetry breaking Higgs 

field has random orientations {{4>'^) pointing in different directions in group space) on 

scales greater than the horizon. One expects therefore to have a probability of order 

unity that a monopole configuration will result after the phase transition (cf. the 

Kibble mechanism). Thus, about one monopole per Hubble volume should arise and 

we have for the number density rimonop ~ ~ T^/rrip, where is the critical 

temperature, when the transition occurs. We also know the entropy density at this 

temperature, s ~ Tj^, and so the monopole to entropy ratio is rimonop/s — 100(Tc/mp)^. 

In the absence of non-adiabatic processes after monopole creation this constant ratio 

^^The isomorfism between two groups is noted as =. Note that by using the theorem we therefore 
can reduce the computation of 1^2 for a coset space to the computation of tti for a group. A word of 
warning: the focus in this thesis is on the physics and the mathematically-oriented reader should bear 
this in mind, especially when we will become a bit sloppy with the notation. In case this happens, 
consult the book [Steenrod, 1951] for a clear exposition of these matters. 




(2.53) 
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determines their present abundance. For the typical value Tc ~ 10^^ GeV we have 
nmonop/ s ~ 10^^^. This estimate leads to a present flmonoph"^ — 10^^, for the superheavy 
monopoles rrimonop — 10^^ GeV that are createcQ. This value contradicts standard 
cosmology and the presently most attractive way out seems to be to allow for an early 
period of inflation: the massive entropy production will hence lead to an exponential 
decrease of the initial nmonopis ratio, yielding flmonop consistent with observations.^ 
In summary, the broad-brush picture one has in mind is that of a mechanism that 
could solve the monopole problem by 'weeping' these unwanted relics out of our sight, 
to scales much bigger than the one that will eventually become our present horizon 
today. 

Note that these arguments do not apply for global monopoles as these (in the 
absence of gauge fields) possess long-range forces that lead to a decrease of their 
number in comoving coordinates. The large attractive force between global monopoles 
and antimonopoles leads to a high annihilation probability and hence monopole over- 
production does not take place. Simulations performed by Bennett & Rhie [1990] 
showed that global monopole evolution rapidly settles into a scale invariant regime 
with only a few monopoles per horizon volume at all times. 

Given that global monopoles do not represent a danger for cosmology one may 
proceed in studying their observable consequences. The gravitational fields of global 
monopoles may lead to matter clustering and CMB anisotropics. Given an average 
number of monopoles per horizon of ~ 4, Bennett & Rhie [1990] estimate a scale 
invariant spectrum of fluctuations {5p/ p)h ~ SOGr^^ at horizon crossing]^. In a subse- 
quent paper they simulate the large-scale CMB anisotropics and, upon normalisation 
with COBE-DMR, they get roughly Grf ~ 6 x 10~^ in agreement with a GUT energy 
scale fj [Bennett & Rhie, 1993]. 

^^These are the actual figures for a gauge SU{5) GUT second-order phase transition. PreskiU 
[1979] has shown that in this case monopole antimonopole annihilation is not effective to reduce 
their abundance. Guth & Weinberg [1983] did the case for a first-order phase transition and drew 
qualitatively similar conclusions regarding the excess of monopoles. 

^•^The inflationary expansion reaches an end in the so-called reheating process, when the enormous 
vacuum energy driving inflation is transferred to coherent oscillations of the inflaton field. These 
oscillations will in turn be damped by the creation of light particles whose final fate is to thermalise 
and reheat the universe (cf. § ^.1.6| ). 

^"^The spectrum of density fluctuations on smaller scales has also been computed. They normalise 
the spectrum at 8h~^ Mpc and agreement with observations lead them to assume that galaxies are 
clustered more strongly than the overall mass density, this implying a 'biasing' of a few [see Bennett, 
Rhie & Weinberg, 1993 for details]. 
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Let us concentrate now on domain walls, and briefly try to show why they are not 
welcome in any cosmological context (at least in the simple version we here consider - 
there is always room for more complicated (and contrived) models). If the symmetry 
breaking pattern is appropriate at least one domain wall per horizon volume will be 
formed. The mass per unit surface of these two-dimensional objects is given by ~ 
_\i/2^3^ where A as usual is the coupling constant in the symmetry breaking potential 
for the Higgs field. Domain walls are generally horizon-sized and therefore their mass 
is given by ~ \^/'^ff'H~'^ . This implies a mass energy density roughly given by pow ~ 
7ft~^ and we may readily see now how the problem arises: the critical density goes 
as Peril ~ which implies VLj:)y/{t) ~ {'fj/mp)'^fjt. Taking a typical GUT value for 
fj we get flr)w{t ~ 10^'^^ sec) ~ 1 already at the time of the phase transition. It is 
not hard to imagine that today this will be at variance with observations; in fact we 
get Qowit ~ lO^^sec) ~ 10^^. This indicates that models where domain walls are 
produced are tightly constrained, and the general feeling is that it is best to avoid 
them altogether [see Kolb & Turner, 1990 for further details]. 



2.2.5 Are defects inflated away? 

It is important to realise the relevance that the Kibble's mechanism has for cosmology; 
nearly every sensible grand unified theory (with its own symmetry breaking pattern) 



predicts the existence of defects. We have seen in § 2.1.3 how an early era of inflation 
helped in getting rid of the unwanted relics. One could well wonder if the very same 
Higgs field responsible for breaking the symmetry would not be the same one respon- 
sible for driving an era of inflation, thereby diluting the density of the relic defects. 
This would get rid not only of (the unwanted) monopoles and domain walls but also 
of any other (cosmologically appealing) defect. Let us sketch why this actually does 
not occur (cf. |3^). Take first the symmetry breaking potential of Eq. ( |2.51|) at zero 
temperature and add to it a harmless ^-independent term 3m^/(2A). This will not 
affect the dynamics at all. Then we are led to 



V{<P) = ^[r-v') , (2.54) 

with f] = (6m^/A)^/^ the symmetry breaking energy scale, and where for the present 
heuristic digression we just took a real Higgs field. Consider now the equation of 
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motion for (f), 

A 

~ — — = —0^ - m^4> -m^0 , (2.55) 
0(p 3! 

for << fj very near the false vacuum of the effective mexican hat potential and 
where, for simplicity, the expansion of the universe and possible interactions of (f) with 
other fields were neglected. The typical time scale of the solution is r ~ m~^. For an 
inflationary epoch to be effective we need r >> H~^, i.e., a sufficiently large number of 
e-folds of slow-roUing solution. Note, however, that after some e-folds of exponential 
expansion the curvature term in the Priedmann equation becomes subdominant and 
we have H"^ ~ SnG V{0)/3 ~ (27rm^/3)(77/mp)^. So, unless fj > mp, which seems 
unlikely for a GUT phase transition, we are led to r << and therefore the amount 
of inflation is not enough for getting rid of the defects generated during the transition 
by hiding them well beyond our present horizon. 

2.2.6 Cosmic strings 

Cosmic strings are without any doubt the topological defect most thoroughly studied, 
both in cosmology and solid-state physics (vortices). The canonical example, describ- 
ing flux tubes in superconductors, is given by the Lagrangian 

^ = -\f,.F^'' + \\D,ct>\' - ^ (I0r - fl'Y ' (2-56) 

with Ffjti, — d[fj^Atj-\, where Aj, is the gauge field and the covariant derivative is — 
+ leA^, with e the gauge coupling constant. This Lagrangian is invariant under 
the action of the Abelian group G = f/(l), and the spontaneous breakdown of the 
symmetry leads to a vacuum manifold M. that is a circle, S^, i.e., the potential is 
minimised for = fjexp{i6), with arbitrary < ^ < 27r. Each possible value of 6 
corresponds to a particular 'direction' in the field space. 

Now, as we have seen earlier, due to the overall cooling down of the universe, 
there will be regions where the scalar field rolls down to different vacuum states. The 
choice of vacuum is totally independent for regions separated apart by one correlation 
length or more, thus leading to the formation of domains of size ^ ~ fj~^. When these 
domains coalesce they give rise to edges in the interface. If we now draw a imaginary 
circle around one of these edges and the angle 9 varies by 27i then by contracting this 
loop we reach a point where we cannot go any further without leaving the manifold 
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Ai. This is a small region where the variable 6 is not defined and, by continuity, the 
field should be = 0. In order to minimise the spatial gradient energy these small 
regions line up and form a line-like defect called cosmic string. 

The width of the string is roughly m^^ ~ (V^^/)~^, rn^ being the Higgs mass. 
The string mass per unit length, or tension, is /i ~ r^^. This means that for GUT 
cosmic strings, where fj ~ 10^^ GeV, we have Gfi ~ 10~^. We will see below that the 
dimensionless combination G/x, present in all signatures due to strings, is of the right 
order of magnitude for rendering these defects cosmologically interesting. 

There is an important difference between global and gauge (or local) cosmic strings: 
local strings have their energy confined mainly in a thin core, due to the presence of 
gauge fields that cancel the gradients of the field outside it. Also these gauge fields 
make it possible for the string to have a quantised magnetic fiux along the core. On 
the other hand, if the string was generated from the breakdown of a global symmetry 
there are no gauge fields, just Goldstone bosons, which, being massless, give rise to 
long-range forces. No gauge fields can compensate the gradients of 4> this time and 
therefore there is an infinite string mass per unit length. 

Just to get a rough idea of the kind of models studied in the literature, consider 
the case G = 5*0(10) that is broken to H = SU{5) x Z2. For this pattern we have 
ni(A^) = Z2, which is clearly non trivial and therefore cosmic strings are formed 
[Kibble et ai., 1982] .[] 

2.2.7 Observational features from cosmic strings 

Let us finish this brief account of cosmic strings by providing just a quick description 
of their remarkable cosmological features. 

Many of the proposed observational tests for the existence of cosmic strings are 
based on their gravitational interactions. In fact, the gravitational field around a 
straight static string is very unusual [Vilenkin, 1981]. A test particle in its vicinity 
feels no Newtonian attraction; however the existence of a deficit angle A = SirGfi 
makes the topology of space around the string that of a cone. Two particles initially at 
rest while the string is far away, will suddenly begin moving towards each other after 
the string has passed between them. Their head-on velocities will be proportional 

^^In the analysis one uses the fundamental theorem stating that, for a simply-connected Lie group 
G breaking down to H, we have 7ri(G/H) = 7ro(H); see [Hilton, 1953]. 
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to A. Hence, the moving string will built up a wake of particles behind it that may 
eventually form the 'seed' for accreting more matter into sheet -like structures. 

Also, the peculiar topology around the string makes it act as a cylindric gravita- 
tional lens that may produce double images of distant light sources, e.g., quasars. The 
angle between the two images produced by a typical GUT string would be oc Gfi and 
of order of a few seconds of arc, independent of the impact parameter and with no 
relative magnification between the images. 

The situation gets even more interesting when we allow the string to have small- 
scale structure, called wiggles, as in fact simulations indicate. The wiggles modify 
its effective mass per unit length, /i, and also built up a Newtonian attractive term 
in the velocity boost inflicted on nearby test particles. These wiggles also produce 
inhomogeneities in the wake of accreting matter and may lead to the fragmentation of 
the structure. The 'top-down' scenario of structure formation thus follows naturally 
in a universe with strings. 

If cosmic strings really exist, the anisotropics in the CMB they produce would have 
a characteristic signature. One such effect comes about due to the Doppler shift that 
background radiation suffers when a string intersects the line of sight. The conical 
topology around the string will produce a differential redshift of photons passing on 
different sides of it, ^ ~ SnGjlvy, with 7 = (1 — t>^)~^/^ the Lorentz factor and v 
the velocity of the moving string. This 'stringy' signature was first studied by Kaiser 
& Stebbins [1984] and Gott [1985]. We will return to it in Chapter |^ where we will 
also give a more detailed description of this effect, and implement it within an analytic 
model for the study on the CMB excess kurtosis parameter. 



2.2.8 Global textures 

Whenever a global non-Abelian symmetry is spontaneously and completely broken 
(hopefully at a grand unification scale), global defects called textures are generated. 
Theories where this global symmetry is only partially broken do not lead to global 
textures, but instead to global monopoles and non-topological textures. As we al- 
ready mentioned global monopoles do not suffer the same constraints as their gauge 
counterparts (cf. § |2.2.4[ ): essentially, having no associated gauge fields they do not 



possess magnetic charge and therefore cannot be accelerated by the magnetic field in 
our galaxy, unlike the gauge monopoles. Thus there will be no 'evaporation' flux of 
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monopoles from the halo of the galaxy and the Parker limit [Parker, 1970] does not 
apply. On the other hand, non-topological textures are a generalisation that allows 
the broken subgroup H to contain non-Abelian factors. It is then possible to have 
TTs trivial as in, e.g., 5*0(5) — > 5*0(4) broken by a vector, for which case we have 
Ai = 5^, the four-sphere [Turok, 1989]. Having explained this, let us concentrate in 
global topological textures from now on. 

Textures, unlike monopoles or cosmic strings, are not well localised in space. This 
is due to the fact that the field remains in the vacuum everywhere, in contrast to 
what happens for other defects, where the field leaves the vacuum manifold precisely 
where the defect core is. Since textures do not possess a core, all the energy of the field 
configuration is in the form of field gradients. This fact is what makes them interesting 
objects only when coming from global theories: the presence of gauge fields could 
(by a suitable reorientation) compensate the gradients of (p and yield -0/^0 = 0, hence 
canceling out (gauging away) the energy of the configuration^ 

One feature endowed by textures that really makes these defects peculiar is their 
being unstable to collapse. The initial field configuration is set at the phase transition, 
when develops a nonzero vacuum expectation value. lives in the vacuum manifold 
Ai and winds around in a non-trivial way on scales greater than the correlation 
length, ^ < t. The evolution is determined by the nonlinear dynamics of 0. When 
the typical size of the defect becomes of the order of the horizon, it collapses on itself. 
The collapse continues until eventually the size of the defect becomes of the order 
of fj~^, and at that point the energy in gradients is large enough to raise the field 
from its vacuum state. This makes the defect unwind, leaving behind a trivial field 
configuration. As a result ^ grows to about the horizon scale, and then keeps growing 
with it. As still larger scales come across the horizon knots are constantly formed, 
since the field points in different directions on Ai in different Hubble volumes. This 
is the scaling regime for textures, and when it holds simulations show that one should 
expect to find of order 0.04 unwinding collapses per horizon volume per Hubble time 
[Turok, 1989]. However, unwinding events are not the most frequent feature [Borrill et 

^^This does not imply, however, that the classical dynamics of a gauge texture is trivial. The 
evolution of the (j)-Afj, system will be determined by the competing tendencies of the global field to 
unwind and of the gauge field to compensate the gradients. The result depends on the characteristic 
size L of the texture: in the range m'^^ « L « vtl^ ^ (e??)"^ the behaviour of the gauge texture 
resembles that of the global texture, as it should, since in the limit niA very small (e 0) the gauge 
texture turns into a global one [Turok & Zadrozny, 1990]. 
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al, 1994], and when one considers random field configurations without an unwinding 
event the number raises to about 1 collapse per horizon volume per Hubble time. We 
will be using these results from the simulations in Chapter ^ when we will consider an 
analytic texture model. 

During the radiation era, and when the correlation length is already growing with 
the Hubble radius, the texture field has energy density ptexture ~ (V0)^ ~ fj^/H'"^, 
and remains a fixed fraction of the total density pc ~ t"^ yielding VLtexture ~ SGfj'^. 
Thus we do not need to worry about textures dominating the universe. 

Needless to say, the texture collapse generates perturbations in the metric of space- 
time. These in turn will affect the photon geodesies leading to CMB anisotropies, the 
clearest possible signature to probe the existence of these exotic objects being the ap- 
pearance of hot and cold spots in the microwave maps. Due to the scaling behaviour 
mentioned above, the density fluctuations induced by textures on any scale at hori- 
zon crossing are given by {dpj p)h ~ ^G'tf. CMB temperature anisotropies will be of 
the same amplitude. Numerically-simulated maps, with patterns smoothed over 10° 
angular scales, by Bennett & Rhie [1993] yield, upon normalisation to the COBEr- 
DMR data, a dimensionless value Gff ~ 10~^, in good agreement with a GUT phase 
transition energy scale. 

2.2.9 Evolution of global textures 

We mentioned earlier that the breakdown of any non-Abelian global symmetry led 
to the formation of textures. The simplest possible example involves the breakdown 
of a global SU{2) by a complex doublet 0", where the latter may be expressed as a 
four-component scalar field, i.e., a = 1...4. We may write the Lagrangian of the 
theory much in the same way as was done in Eq. (|2.56|) , but now we drop the gauge 
fields (thus the covariant derivatives become partial derivatives). Let us take the 
symmetry breaking potential as follows, = ^ — v'^)'^- The situation in which 
a global SU{2) in broken by a complex doublet with this potential V is equivalent to 
the theory where 50 (4) is broken by a four-component vector to S0{3), by making 
(j)'^ take on a vacuum expectation value. We then have the vacuum manifold given 
by S0{4:)/ S0{3) = S^, namely, a three-sphere with 0'^0a = fj"^- As 773(6'^) 7^ 1 (in 
fact, 773(5*^) = Z) we see we will have non-trivial solutions of the field (p"" and global 
textures will arise. 
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As usual, variation of the action with respect to the field 0*^ yields the equation of 
motion 

0^" + 2^/-VV = -a^^, (2.57) 

where primes denote derivatives with respect to conformal time and V is computed in 
comoving coordinates. When the symmetry in broken three of the initially four degrees 
of freedom go into massless Goldstone bosons associated with the three directions 
tangential to the vacuum three-sphere. The 'radial' massive mode that remains {m^ ~ 
\/Xfj) will not be excited, provided we concentrate on length scales much larger than 

To solve for the dynamics of the field (p^, two different approaches have been im- 
plemented in the literature. The first one faces directly the full equation ( |2.57|) , trying 



to solve it numerically. The alternative to this exploits the fact that, at temperatures 
smaller than Tc, the field is constrained to live in the true vacuum. By implementing 
this fact via a Lagrange multiplieiQ we get 



V^V.^" = - ; <f = f,\ (2.58) 



with the covariant derivative operator. Eq. (|2.58| ) represents a non-linear sigma 



model for the interaction of the three massless modes [Rajaraman, 1982]. This last 
approach is only valid when probing length scales larger than the inverse of the mass 
m^^. As we mentioned before, when this condition is not met the gradients of the field 
are strong enough to make it leave the vacuum manifold and unwind. 

The approach (cf. Eqs. ( p.58| )) is suitable for analytic inspection. In fact, an exact 



fiat space solution was found assuming a spherically symmetric ansatz. This solu- 
tion represents the collapse and subsequent conversion of a texture knot into massless 
Goldstone bosons, and is known as the spherically symmetric self-similar (SSSS) exact 
unwinding solution. We will say no more here with regard to the solution, but just 
refer the interested reader to the original articles [see, e.g., Turok & Spergel, 1990; 
Notzold ,1991]. In what follows we will be mainly interested in considering the results 
of the numerical simulations from Eq. (|2.57|) , the reason being that these simulations 
take full account also of the energy stored in gradients of the field, and not just in the 



^^In fact, in the action the couphng constant A of the 'mexican hat' potential is interpreted as the 
Lagrange multipher. 
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unwinding events. We will apply the results of the simulations performed in [Durrer 
& Zhou, 1995] in Chapter ^ below. 
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Chapter 3 

CMB anisotropies 



The standard method of characterising the CMB fluctuation spectrum is by means 
of a multipole moment expansion: current experiments measure the distribution of 
temperatures in the sky, AT/T(7), where 7 points in a certain direction. Similar to 
what one does in Fourier analysis, namely, the study of a function by means of its 
expansion in plane waves, in our case we expand any given function deflned on the 
sky (like AT/T) in spherical harmonics, Y^. Since we are only able to measure mean 
values for the stochastic distribution of temperature fluctuations, and since physical 
quantities should not depend on the particular direction we are looking at (isotropy 
assumption on very large scales), we see that at the end of the day only the multipole 
index i will be relevant 

The temperature autocorrelation function tells us about possible anisotropies in 
the radiation coming from points in the sky separated by an angle 'd, and is given by 

AT AT 1 
{C^m = (^(7)^(7')) = ^ E(2^ + l)C.P.(cos^) (3.1) 

where (■) means averaging on the sky and 7 ■ 7' = cos'd. 

The quantity i{£ + l)Ce is usually referred to as the angular power spectrum. With 
this deflnition, an initial distribution of density fluctuations with n = 1 produces a flat 
(^-independent) spectrum on large angular scales {i small), namely, a 'plateau' in the 
i{i + l)Ci vs. £ plot. Notice that both inflation and topological defect models lead to 
approximately scale invariant spectra on large scales; therefore, observations on large 

^We are necessarily being a bit loose at this stage in trying to give a brief introduction. A detailed 
account of these matters will be given in Chapter 
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scales alone cannot by themselves select one, and not the other, of these competing 
models. 

Interestingly enough, going to smaller scales it turns out that power spectra pre- 
dicted from inflationary and defect models are not similar: the Doppler peaks expected 
on angular scales ~ 1° might be a helpful tool to discriminate between them; see Chap- 



On even smaller scales the details of the process of recombination become impor- 
tant: the last scattering surface has a finite thickness, which is due to the noninstanta- 
neous recombination of hydrogen, and to the residual ionisation level remaining. The 
finite width of this shell smooths out fluctuations on scales ~ 8fi^^/^ arcmin or smaller. 

The detection and detailed analysis of temperature anisotropies on different angular 
scales provide a unique measure of the primordial density fluctuations from which 
large-scale structure evolved, as well as a direct probe of the as yet uncertain values of 
the cosmological parameters. In the following sections we will enumerate the different 
sources of temperature anisotropies as well as the physical processes responsible for 
their suppression, together with a brief explanation of the physics behind them. 

3.1 Contributions to ^ 

Several are the authoritative reviews on this subject [e.g., Efstathiou, 1990; White et 
al., 1994], and we refer the reader to them for more detailed presentations of the items 
treated here. In this section we will briefly account for the main effects that may leave 
their imprint in the CMB temperature anisotropies. 

Dipole anisotropy 

The dipole anisotropy was for a long time, until the COBE detection, the only mea- 
sured temperature variation. Both 'intrinsic' and 'extrinsic' contributions are present 
in the maps and that is why it is source of large uncertainties. Essentially, the extrinsic 
part is 



where 6 is the angle between the direction of motion and that of observation, and 
is mainly due to our peculiar motion with respect to the comoving rest frame of the 
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CMB radiation. The intrinsic dipole on the other hand is expected to be two orders 
of magnitude smaller and of comparable amplitude to the quadrupole.0 

This dipole anisotropy was detected and from its trace in the ^ maps both v and 
the direction of our peculiar motion can be determined. The best estimate to date is 
from COBE which gives a speed of 365 ± 18 km sec~^ [Smoot et al, 1991]. When this 
velocity is decomposed and i) the motion of the Earth around the Sun, ii) the motion of 
the Sun in the galaxy and, finally, iii) the motion of the Milky Way towards Andromeda 
in the Local Group (LG) are subtracted, the motion of our LG may be determined to 
be ~ 600 km sec~^ towards the direction (in galactic coordinates) i ~ 270°, b ~ 3O°.0 

Sachs— Wolfe effect 

The second effect, and the most relevant one on very large scales, is the Sachs- 
Wolfe effect [Sachs & Wolfe, 1967]. When focusing on these perturbations we are 
facing rather 'simple' physics: since the scales involved are much larger than the size 
of the horizon at last scattering the analysis is to a great extent transparent to the 
microphysics of recombination, and consequently also to some of the still unknown 
cosmo logical parameters, like Qb, ^cdm, ^hdm, ^a, h, and ionisation history.^ 

The physics behind Sachs-Wolfe fluctuations is straightforward: matter pertur- 
bations distort space, and therefore also distort geodesies of photons that were last 
scattered at the recombination era. The net effect can be expressed as ^ ~ 
Since the perturbations in the gravitational potential $ are time independent in the 
linear regime, we only need to worry about $ on the last scattering surface. Blueshift 
and redshift cancellation effects ensure that distortions to the CMB photons on their 
way to us are unimportant .0 



^On the same footing, the amplitude of the extrinsic quadrupole is a factor ^ ^ that of the extrinsic 
dipole. 

•^For a comprehensive introduction to these matters consult Padmanabhan [1993]. 

''On the other hand results for these scales are sensitive to other parameters such as relative (am- 
plitude) contributions to ^ from energy density (scalar metric), gravitational wave (tensor metric), 
isocurvature scalar and other primordial fluctuations, and to the shape of the initial fluctuation spectra 
(as parameterised by the spectral indexes of the different components). 

^Note however that we expect perturbations to go non-linear in the end, and therefore intervening 
clumps will undoubtedly play a role - see below. 
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Intrinsic fluctuations 

In this paragraph we will consider fluctuations intrinsic to the radiation fleld it- 
self on the last scattering surface. These vary according to whether the primordial 
fluctuations are adiabatic or isocurvature in nature. 

Adiabatic perturbations are those for which the densities in matter and radiation 
(and by the equivalence principle also all other components of the energy density) are 
perturbed so that the entropy per baryon is the same as in the unperturbed state.0 

Since the radiation fleld energy density goes like prad oc (and remembering that 
the density contrast is deflned as 5 = we have ^ = \6rad = l^mat, where the last 
equality follows from the adiabatic condition. Now, since the CMB photons start off 
their trip during matter domination, one has 6p = Pmat^mat + Prad^rad — P^mat- Thus 
we get ^ = 1^ for adiabatic perturbations. 

On the other hand, isocurvature perturbations do not affect the total energy density 
and are therefore characterised hj 6p = [see, e.g., Kolb & Turner, 1990 for a review]. 
They correspond to perturbations in the form of the local equation of state. Let us 
consider a two-fluid system, e.g., a non-relativistic CDM component plus radiation. 
Then, we have p = rricdmncdm + prad and thus ^ = -\{pcdm/ prad){5ncdm/ncdm)- We 
call s the entropy density (s oc is proportional to the number density of relativistic 
particles) and ricdm the number density of CDM particles, thus s/ricdm is a measure of 
the entropy per particle. 

It turns out to be convenient to work with the following perturbation in the number 
density 6{ncdm/ s) / {ucdm/ s) = S{ncdm)/{ncdm) - 3^, and with this we get 

AT _ I 5{ncdm/s) f jPcdm/Prad) \ o\ 
T ~ 4 (n^dm/s) \{l + lpcdm/prad)) ' 

Now, at late times the universe becomes matter dominated pcdm >> Prad, and thus 
~ —\6{ncdm/ s) / {ncdm/ s). When a given mode becomes smaller than the horizon, 
microphysical processes play a role and fluctuations in the local pressure can generate 
fluctuations in the energy density. From that moment onwards the difference between 

^For an initial state to be adiabatic (or isoentropic) it is necessary that, for any two types of 
components, say A and B, the ratio of the number density of the A particles to the B particles, 
nA/nB be a constant independent of position. The origin of this terminology draws back to the 
universe with just radiation and plasma content, wherein nrad/ns is a measure of the entropy per 
baryon. 
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adiabatic and isocurvature is irrelevant. Further, the perturbation in the number 
density becomes a density perturbation of the same amphtude S^ricdm/ s) / {ricdm/ s) ~ 
5p/p yielding finally ^ ~ —^6p/p. 

Infiationary models are normally invoked as examples of a scenario for the genera- 
tion of adiabatic perturbations, by means of quasi-de Sitter fiuctuations in the infiaton. 
On the other hand, cosmic strings provide an example of isocurvature perturbations 
(of course, not the only one, but an interesting one): during the phase transition when 
these defects are produced, some of the radiation energy density goes into the forma- 
tion of the network of strings. Thus, the total energy density remains unaltered while 
a change in the effective equation of state of the matter takes place in those regions 
where the defects arise (we will say more on strings and their infiuence on the CMB 
anisotropies in §7.1|). 

Acoustic oscillations and the Doppler peaks 

On smaller scales, say ~ 1°, one is effectively probing those perturbations with char- 
acteristic scales of order of the horizon distance, or smaller than it, at recombination. 
Microphysical processes play definitely an important role now. 

Baryons and photons are still tightly coupled through Compton scattering and may 
be studied as a single fiuid. Oscillations in a given mode of the density perturbations 
are expected when pressure forces dominate over gravity. While in the case of CDM 
the fiuid is pressureless and thus always gravity wins, in generic fiuids the behaviour 
of an arbitrary mode will depend strongly on the sound speed in the particular fiuid. 
This latter in turn determines the Jeans length Aj, which is precisely that perturbation 
wavelength for which the pressure and gravity forces exactly balance (refer to § [1.4|) : 
near recombination we have Aj ~ 27iCsH~^, with Cg ~ Once the wavelength 

of a given mode falls below the Jeans length, pressure forces take over and therefore 
acoustic oscillations appear. 

These oscillations arise both in the density and velocity perturbations, and together 
lead to the existence of peaks in the i{i+l)Ce vs. i plot as remnants of the time at which 
matter and radiation decoupled. Even though the terminology might be somewhat 
misleading, these features in the angular power spectrum are known under the name 
of 'Doppler peaks', although no Doppler effect is taking place since the components 
are oscillating together and little or no velocity difference exists between them. 
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Sunyaev-Zel'dovich fluctuations 

Non primeval in origin but at the root of spectral distortions detected in a few 
reach clusters, Sunyaev-Zel'dovich fluctuations have been very well studied [see, e.g., 
Zel'dovich & Sunyaev, 1969; Sunyaev & Zel'dovich, 1970]. A generic feature in struc- 
ture formation scenarios is the existence of large regions of tenuous hot gas, left behind 
in the process of clumping and subsequent fragmentation of matter. The Sunyaev- 
Zel'dovich effect arises as a consequence of the CMB photons passing through these 
clouds of hot electrons, where they are (inverse Compton) scattered and get distorted. 
The effect may be expressed as an equivalent temperature perturbation as a function 
of the wavelength, and on the Ray leigh- Jeans side yields 



where ks is the Boltzmann constant, and and rrie refer to the electron's temperature 
and mass, respectively. 

Secondeiry anisotropies: the Vishniac effect 

If the binding energy of the gaseous components of newly bound structures is 
thermalised, then temperature fluctuations may result on sub-arcminute scales. We 
will see below that primary fluctuations might be greatly suppressed on these scales 
due to the smearing associated with the finite thickness of the last scattering surface, 
and thus newly-generated secondary anisotropies may dominate. While in principle 
there are several sources in second-order perturbation theory, there is a contribution, 
coming from a product of velocities and densities, that dominates over all others [e.g., 
Hu et al, 1994; Dodclson & Jubas, 1995]; this is known as the Vishniac effect [Vishniac, 
1987]. This contribution turns out to be negligible for standard ionisation histories; it 
is however important whenever we allow for late reionisation, as infall of baryons into 
the CDM perturbation wells plus the growth of the perturbation in the intervening 
period makes this a relevant source of anisotropies. 
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3.2 Suppression of anisotropies 

Width of the last scattering surface 

We explained before (e.g., in § |1.3| ) how the expansion of the universe makes the 
mean energy density gradually to decrease, and that this produces the decoupling of 
radiation from matter (due to the lower abundance of partners to interact in order to 
keep the equilibrium). Now, it may well happen that the last scattering surface was 
in fact not sudden, i.e., that it took some time for the majority of photons to freely 
stream away from their interactions with matter. This would be the most plausible 
situation and thus we see that this surface actually has a 'width', which means that 
every time we look in some direction on the sky we are in fact detecting a bunch of 
photons that were emitted from different points along the line of sight. This implies a 
kind of averaging of the fluctuations, and a destructive interference between the modes 
with wavelengths smaller than the thickness of the last scattering surface, with the 
result of a net 'washing out' of perturbations on very small angular scales, of order 
8n~^/^ arcmin. 

Silk damping 

When focusing on very small scales it becomes apparent that baryons and radiation 
are in reality imperfectly coupled. Photons possess a mean free path Ac in the baryons 
due to Compton scattering. Due to the frequent collisions, suffered by a large part of 
the photons, they randomly walk through the baryons with the result that hot and cold 
spots get mixed. Fluctuations only remain in a relatively small unscattered fraction of 
them, causing a nearly exponential decrease in amplitude as the characteristic scale of 
the perturbation becomes smaller than the diffusion length ~ \/N\c- Photons tend 
to stream away of the otherwise growing perturbations and in so doing they smooth 
the distribution of fluctuations (photon diffusion). This effect ends up by suppressing 
the level of anisotropies on very small scales, and is known as the Silk damping [Silk, 
1968]. 

Both the flnite thickness of the last scattering surface and Silk damping are the 
most eflicient mechanisms for the suppression of anisotropies on very small angular 
scales. We may see this in Figure |3.1| below; these effects are the responsible for the 
suppression of the angular power spectrum for ^ > 1000. 
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Early reionisation 

At last scattering the mean free path of the photons goes like oc {xenf,)~^. Ac 
increases due to the decrease in the ionisation fraction Xe during recombination [Hu, 
Sugiyama & Silk, 1995]. From this we see that the diffusion scale is dependent on 
the ionisation history (as well as on the number density in baryons, nt) and therefore 
a delay in last scattering, as might happen if early reionisation occurs, may lead to 
an extended photon diffusion with the resulting destruction of anisotropics, specially 
at the degree scale. The acoustic Doppler peaks are the dominant feature on these 
small scales and the net effect of reionisation would be a reduction (if not a complete 
suppression) of them. On the other hand, for large angular scales we do not expect the 
ionisation history to play any relevant role. Therefore, the normalisation of the CMB 
spectrum on large scales is safe from these uncertainties. 

3.3 Adiabatic fluctuations as source of ^ 

The dominant contribution to large scale fluctuations comes from the Sachs-Wolfe 
effect. Sachs & Wolfe [1967] in their classic paper on linear perturbations around 
an Einstein-de Sitter background universe, derived a general formula for the CMB 
anisotropy associated with scalar and metric perturbations. In this perturbed uni- 
verse model each photon is assumed to have been emitted at conformal time T]r (at 
recombination) and detected by us at r/o (present time).[] 

The temperature fluctuation, to first order in the perturbing quantity, may be 
expressed as 



where x^^rj) = •y^^rjo — rf) is the unperturbed photon geodesic, and = (1,7)) with 7 

pointing along the line of sight. In the above equation a uniform source was assumed 

(thus, no intrinsic fluctuations so far - see below), and we also used the correspondence 

/loo = 2$ relating the metric perturbation to the gravitational potential. 

The three terms in Eq. (|3.5|) are identified as (i) the gravitational redshift (that 

may also depend on the position) due to the difference between $'s at recombination 

''In Ref . | |1 76| , recombination was taken as instantaneous. We have seen that the actual width of 
the last scattering surface is not relevant but for very small scales. 
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and now; (ii) Doppler shifts, due to peculiar motions of both emitter and receiver; 
and (iii) a term sensible to time dependences of the metric perturbation. In practice, 
one is not interested in anisotropies produced by peculiar motions of the observer 
which create a pure dipole anisotropy on the sky, whereas we will see in Chapter 
H that perturbations in the peculiar velocity of the fluid components is key for the 
development of the Doppler peaks. The time dependence of hf^i, is relevant for models 
with gravitational waves, and therefore the integral term plays an important role. 
Similarly, for A-dominated or non-flat models in general, (iii) is relevant, as well. 

In the Einstein-de Sitter model and if perturbations in the metric arise only due 
to density perturbations, then the integrand is essentially 48^ /dt]. If we consider 
just growing mode density perturbations in the linear regime, $ is constant in time| 
and therefore the integral vanishes. On the other hand, there have been scenarios 
proposed (such as those where the occurrence of a 'late-time phase transition' greatly 
influences the total temperature anisotropies) where it is in fact the integrated term 
the interesting one. Also, this term will become relevant when studying the (mildly) 
non-linear growth of perturbations and their influence on CMB non-Gaussian features, 
cf. Chapter |] below. 



At this point we should stress that the result of Eq. ( |3.5| ) has not (yet) the usual 
form that one usually meets when considering the Sachs- Wolfe effect, even in the 
absence of tensor metric fluctuations. In fact we still need to take into account the 
intrinsic anisotropies in the radiation fleld, cf. § p.l| . For adiabatic perturbations, 
we found ^ ~ |^ &nd, recalling that an overdensity generates a larger Newtonian 
potential (^ ~ 2$), we see there is a partial cancellation, yielding the familiar result 

~ 1$, where we assumed $ = at the time of observation, and where the Doppler 
and integrated terms in (|3.5|) were not explicitly written down. 



With regards to ^ coming from isocurvature fluctuations, we have seen in §3.1| 
how a positive density perturbation (and the corresponding $) is compensated by a 
negative fluctuation in the temperature of the CMB radiation. Thus, this perturbation 



AT 
T 

with the corresponding adiabatic models, for the same perturbation in the density. We 



adds to the Sachs- Wolfe one, and this yields relatively large as to be compared 



®This applies for perfect fluids in a flat space regardless of whether the universe is dominated by 
matter or radiation. It may be most easily seen from the Poisson equation, and remembering the 
fact that the density contrast grows with the scale factor as 5 cx a (cx a^) for matter (radiation) 
domination. 
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will not go into further details but just refer the interested reader to the reviews in 
[White, Scott & Silk, 1994] and [Stebbins, 1993] for nice accounts of these matters. 

3.4 Non— linear effects 

We have seen how linear perturbations in a flat, matter dominated FRW universe led 
Sachs and Wolfe to derive their well-known formula. Now, what would happen if we 
considered a different cosmo logical model and/or we let the initially small perturba- 
tions go non-linear and form large-scale density inhomogeneities? Rees & Sciama 
[1968] were the first to note that these evolving mass concentrations would make the 
gravitational potential $ vary with time, even for a fixed space point. Thus, radiation 
travelling through these intervening clumps will be perturbed and the integral in Eq. 
( p. 51) will depend on what lies along the photon's path. 

Recently, using ray tracing techniques through N-body CDM simulations, Tuluie 
& Laguna [1995] showed that two effects related to the variations of the gravitational 
potential are relevant on small scales. The studied effects are the intrinsic evolution 
of $ (e.g., changes in the depth of the potential wells) during the time the photons 
take to transverse the protostructure, and the change in the gravitational potential as 
the large-scale structures formed (like clusters) move across the microwave sky. They 
find that both these effects can be of similar amplitude, and amount to a level of 
anisotropics of order 10"'' < ^ < 10~^ on scales of ~ 1°. 

We will revisit integrated effects due to a varying gravitational potential during the 
mildly non-linear growth of perturbations in Chapter ^ below, when we will calculate 
the post-recombination contribution to the CMB three-point temperature correlation 
function. 



3.5 Overview 

Although we have not introduced so far the necessary notation and machinery for 
giving a 'quantitative' account of the many contributions to the CMB temperature 
anisotropics on the various angular scales, let us at least qualitatively depict what 
currently is though to be the 'standard model' for the fluctuations in the relic radiation. 



In Figure |3.1| we show the generic prediction for the angular power spectrum from 
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Multipole Moments - / 



Figure 3.1: Angular power spectra predicted within inflationary cosmologies [from 
Steinliardt, 1995]. Tie spectra are for CDM and cosmological parameters h = 
0.5 , fifi = 0.05, with cosmological constant A = 0. 
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standard inflation plus CDM on aU relevant angular scales, namely, from the very 
large {i small) to the smallest scales (large €s) currently probed by observations. In 
this figure, the upper solid curve comes from a primordial (power-law) spectrum of 
density fluctuations, with spectral index n = 1 on large angular scales (those scales 
which suffered little or no evolution at all) and no gravitational wave contribution. 
Vertical hashmarks give an idea of the theoretical uncertainties inherent to the CMB 
anisotropics which go under the name 'cosmic variance' and, as we see from the plot, 
are specially relevant on large angular scales. The lower curve (dashed) is for n — 0.85 
and shows the apparent 'tilt' in the otherwise flat Sachs-Wolfe plateau; this fact leads 
to spectra with more power on large scale (as to be compared with the 'flat' Harrison- 
Zel'dovich spectrum n — 1). This curve also shows how the corresponding spectrum 
gets modifled when the amplitude on large scales is built up by both scalar and tensor 
modes (50-50 mixture at the quadrupole level in the plot). For scales well within the 
horizon at recombination (as is the case for £ > 150) gravitational waves redshift away 
and thus the overall amplitude gets equally reduced (as to be compared with the case 
with solely scalar contribution to the anisotropics). 
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Figure 3.2: Amplitude of the fluctuations in different experiments, as a function of the 
scale (multipole i ~ 9~^) [from Scott et al, 1995]. Qfiat is the best-fitting amplitude 
of a flat power spectrum, given by i{i + l)Ci = (247r/5)(Q//(it/^o)^- 



While on large angular scales the observational results are well settled, on smaller 



scales the situation is not as clear. In Figure we show the present status regarding 
observations. As it is clear from this figure there is some indication around scales 
i ~ 200 that the level of anisotropies grows, hinting for the existence of the Doppler 
peaks [Scott, Silk & White, 1995]. However, it is also apparent from the not always 
compatible data points and vertical error bars that presently nothing can be definitely 
concluded. More, independent pieces of observations are thus needed to corroborate 
the present experimental results and to fill in the gaps. 
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Chapter 4 
CMB statistics 



4.1 Dealing with spherical harmonics 

In the following sections we will compute correlations of temperature anisotropics. The 
pattern of these anisotropics is distributed on the CMB sky two-sphere and, as such, 
it will prove useful to expand the temperature fluctuation field ^ (as we may do with 
any function that lives on a sphere) in spherical harmonics. 

It is conventional to define the spherical harmonics in terms of the associated Leg- 
endre functions as follows 



with £ and m integers fulfiUing the relations £ >0 and |m| < £. 

As in the expression above we will often use the angles 6*, or the unit vector 
7 = (sin^cos0, sin^sin0,cos^) interchangeably; when working with temperature 
anisotropics 7 will be taken as the line-of-sight direction on the sky. 

Being the product of e*"**^, which form a complete set of orthogonal functions in 
the index m on the interval < < 27r, and Pp{cos9), which form a similar set in 
the index £ for fixed m on the interval — 1 < cos ^ < 1, it is clear that the Y/" will form 
an analogous set on the surface of the unit sphere in the two indices i, m. 

Under complex conjugation they change as Yf"^ = {—l)"^Yp* . The orthonormality 
condition is 

J dn^Yr{l)Yl!'\l) = hl'^ram' (4.2) 



7.^ 
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where the differential sohd angle is given by d^l^ = sinO d9 d(j). These properties, 
together with the completeness relation, allow one to expand any square-integrable 
function on the unit sphere, e.g.. A, as 

oo I 

A(7) = E E <Yni) (4.3) 

£=0m=-£ 

where the complex quantities are given by 

aT = J dn^Vriimi)- (4.4) 

We will be interested in the situation in which A is actually the CMB tem- 
perature fluctuation, and in this case we will call the multipole moments. These 
are relevant quantities as in general predictions of a theory are expressed in terms of 
predictions for the a^. 

Also useful in this work is the spherical harmonics addition theorem 

47r ^ 

^.(cos^) = -— E >^r*(^i,0i)>^r(^2,02) (4.5) 

m=—i 

This theorem applies to any spherical triangle on the unit sphere, with sides di, 62 
and 6*, and with the angle |0i — ^2! opposite to 6*, so that (being 71 = (6*1,01) and 
72 = (6*2, 02)) we have 71-72 = cos6' = cos 6*1 cos 6*2 + sin 6*1 sin ^2 cos |0i — 02|. 

The spherical harmonics as defined above are connected with the bases of the even 
(or single-valued) representations of the rotation group on the sphere.^ 

We will see below that, when dealing with correlations, it will prove convenient to 
define the rotationally symmetric quantities = Y,m=_i I^T^P, which we will call the 
multipole coefEcients]^ 

That the above defined quantities are rotationally symmetric can be seen from the 
fact that a suitable average over the m index was performed. Given the expansion in 
spherical harmonics of any function 

/(7i) = EE<^r(7i), (4.6) 

£ m 

^The odd (double- valued) representations arise for half-integer £ and the corresponding bases are 
called spinors. 

^It is standard in the literature to work in terms of ensemble averages, in order to get theoretical 
predictions from different models. In that case one computes the quantities Ce = {Qj)/{2£+ 1) and 
it is actually to these C^'s that one assigns the name multipole coefficients. 
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one could wonder how the multipole moments change for a rotated version of the 
same function, say written in terms of the same Y£™(7i)'s. That is, we express 

fil2) = J2iJ2m ^''''^^a™^/"(72) and, by recalhng the relation 

Yr{i2) = T.DLmyr'iii) (4.7) 

m' 

we get = J2eJ2m' '^'^^'^'^aT'^e^' ill) where we defined 

(-^-^af = Y.DLm^"''^ciT. (4.8) 

m 

This means that after rotation the new are just a linear combination of the old 
ones, with i index fixed. We see then that a rotation on the sphere will not change the 
multipole coefficients Qj, which depend only on 

Now going back to the defining expression for the Qj, and baring in mind what we 
discussed above, we see that no matter how we rotate a given function the different 
i contributions to it are independent and do not mix: say, if we focus on the dipole 
{i = 1) contribution (the linear part), and we rotate, it will never pick up quadratic 
{£ = 2), cubic {£ = 3) or any other £ 7^ 1 terms. 

4.2 The CMB two— point correlation function 

Our main aim below will be to compute the CMB three-point temperature correlation 
function. The very existence of non-vanishing higher-order correlation functions on 
large angular scales is due to possible departures from a Gaussian behaviour of the 
primordial peculiar gravitational potential (cf. Chapter ^). Also, even for the case of 
primordial Gaussian curvature fluctuations, the mildly non-linear gravitational evolu- 
tion gives rise to a non-vanishing three-point correlation function and we will see in 
Chapter ^ below what predictions we get for this integrated effect. 

As a flrst step we will deflne the connected two- and three-point correlation func- 
tions of the CMB temperature anisotropics as measured by a given observer, i.e., 
on a single microwave sky. Let us then deflne the temperature fluctuation fleld 
^(x;7) = (T(x; 7) — Tq (x) ) /Tq (x) , where x specifles the position of the observer. 



3 



The so-called D-functions are in fact {2£ + 1) x {2i + 1) unitary matrices D^, and they form a 



(2£ + l)-dimensional irreducible representation of the groups SO{3) and SU{2). See Ref.|152] for a 
clear exposition of their properties. 
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the unit vector 7 points in a given direction from x and To (a;) = / ^^T{x;^) repre- 
sents the mean temperature of the CMB measured by that observer (i.e., the monopole 
term) . 

The angular two-point correlation function C2{x;a) measured by an observer 
placed in x is the average product of temperature fluctuations in two directions 71 
and 72 whose angular separation is a; this can be written as 

C,{x;a)^ /^/^<^(7i-72-cosa)^(f;7i)^(f;72). (4.9) 
As well known, in the limit a ^ one recovers the CMB variance 6*2(0;) = 

47r I. T 

AT, 



/ -T-^[^{x;^)]'^. Expanding the temperature fluctuation in spherical harmonics 



(^;7) = E E aT{x)WeYr{^), (4.10) 
e=i m=-e 

and writing the Dirac delta function as a completeness relation for Legendre polyno- 
mials Pi 

00 

5(cos^i - COS ^2) = Y.{i+l)Pi{cosei)Pi{cos92), (4.11) 
e=o 

we easily arrive at the expression 

C2{x;a) = i^E^^(cos«)g|(f)>V,2 , (4.12) 

where Qj = EL=-eK?- 

It is standard to compare theory and experiments through the a^. However real 
measurement of temperature anisotropics are limited by finite resolution and specific 
strategies are employed by the many current experimental settings. How to account 
in our definitions for this modeling of the real world? It has become conventional to 
include this piece of dirty physics by the use of a hltcr, limiting the range in i (and 
hence, the angular scale) to which the experiment is sensitive, in the form of a window 
function. In the previous expressions We represents in fact the window function of the 
specific experiment. Setting Wq = Wi = automatically accounts for both monopole 
and dipole subtraction. 

Let us remind the reader that the subtraction techniques employed in current 
anisotropy experiments do not allow us to get a value for the monopole. Moreover, the 
dipole contribution present in the maps (essentially ^ ~ ^ cos^, where 9 is the angle 
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between the direction of motion and that of observation - see §3.1| )p| is considered to 
be mainly due to our pecuhar motion with respect to the comoving rest frame of the 
CMB radiation. 

We will often concentrate below on large angular scales. In the case of the COBE- 
DMR experiment the window function for £ > 2 is Wi ^ exp — + l)a^ , where 
cr is the dispersion of the antenna-beam profile, which measures the angular response 
of the detector [e.g., Wright et al, 1992]. In some cases the quadrupole term is also 
subtracted from the maps [e.g., Smoot et al, 1992]; in this case we also set W2 = 0.| 

4.3 The CMB three— point correlation function 

The analogous expression for the angular three-point correlation function is obtained 
by taking the average product of temperature fluctuations in three directions 71, 72 
and 73 with fixed angular separations a (between 71 and 72), (3 (between 72 and 
73) and 7 (between 71 and 73); these angles have to satisfy the obvious inequalities 
|tt — 7I < /5 < a + 7- One then has 



C^ix; a, (3,'y) = / — — / / (icos'(9i25(cos'i9i2 — cosa) / (icost923 

J Att Jo 2n J-i J~i 

/■I AT AT AT 

x5(cos?923 - COS/9) / rfcos^9i3(5(cos79i3 - C0S7)— ^(f;7i)— ^(f;72)^7r(^;73), 
J-i 111 

(4.13) 

where cos 'dai3 = 7a ■ 7/3 and (fi2 is the azimuthal angle of 72 on the plane orthogonal 
to 7i. The above relation can be rewritten in a form analogous to Eq.( |4.9|) , namely 



C,{x; a, (3, 7) = N{a, f^^l) J ^ J ^ J " 72 - cos a) 

AT AT AT 

X(5(72-73 -cos/3)5(7i-73-cos7) — (f;7i) — (f;72) — (f;73) , (4.14) 



^This contributes with a i = 1 term in the expansion in spherical harmonics. 

^When the first analyses of the DMR maps were performed it was noted that the quadrupole was 
also easy target of second-order Doppler effects, and these could contaminate the former with a rms 
kinematic quadrupole of order 1.3 . This, together with the high intrinsic theoretical uncertainty 
(cosmic variance, see below) of the lower order multipoles, led the team to subtract, not only the 
mean and dipole, but also the quadrupole from the maps. It was afterwards recognised that this in 
fact makes little difference in the maximum likelihood analyses done for estimating the spectral index 
n and the Qrms-ps normalisation [e.g., Gorski et al, 1994], and it is standard in current literature 
to state both fittings, with and without 1 = 2 contribution. 
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where iV(a,/?, 7) = ^/l — cos^o; —cos'^P —008^7 +2 cos a cos /? C0S7.Q Setting a = 
/3 = 7 = in these general expressions one obtains the CMB skewness C^lx) = 
/ ^^[^(x; 7)]^. Also useful are the egm/atera/ three-point correlation function [e.g., 
Falk et al, 1993] and the collapsed one [e.g., Hinshaw et ai., 1994 and references 
therein], corresponding to the choices a = /5 = 7, and a = 7, /3 = 0, respectively. 
Alternative statistical estimators, more suited to discriminate bumpy non-Gaussian 
signatures in noisy data, have been recently introduced by Graham et al. [1993]. In 
all the above formulas, full-sky coverage was assumed, for simplicity. The effects of 
partial sky coverage on some of the statistical quantities considered here are discussed 
in detail by Scott, Srednicki & White [1994] (see also Ref.^). 

Following the procedure used above for C2{x;a), we can rewrite the three-point 
function in the form 

xE E P,(cosa)P,(cos/5)PKcos7)HX"HSr"^'^r/r'^^ 

where the coefficients UJ^'^f^'' = / dfi^F™^* (7)^^72 (7)1^^7(7), which can be easily 
expressed in terms of Clebsch-Gordan coefficients [e.g., Messiah 1976], are only non- 
zero if the indices £j, rrii {i = 1,2,3,) fulfill the relations: \ij — ik\ < < \^j + h\, 
^1 + ^2 + ^3 = even and mi = 1712 + m^. We may occasionally in what follows make use 
of a related quantity UTXiT'' = (-1)™'^"^™™' ^^i^h has the further advantage of 
being invariant under circular permutation of its columns (^™'^ . In fact, this is apparent 
from the following expression in terms of the '3j' symbol^ 

^ ((2£i + 1)(2£2 + 1)(2£3 + (i^ {X 1 %) . (4.16) 

The collapsed three-point function measured by the observer in x reads 

^1/2,^3 rni,m2,m3 

which, for a = 0, gives a useful expression for the skewness. 

® To show that the latter expression is correctly normalised one can expand the delta functions 
in Legendre polynomials ( [t.ll[ ), use the spherical harmonics addition theorem ([4.51) , and note the 
relation 2,(2£ + l)P,(x)Pf (y)^^(^) = |(1 - - - + 2xyz)-^l^ (see, e.g., jl7§). 

''This expresses, for the particular case of the spherical harmonics, a more general theorem, called 
the Wigner-EcJcart theorem. 
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4.4 Predictions of a theory: ensemble averages 

To obtain definite predictions for tlie statistics described above, one needs to exploit 
tlie random nature of tlie multipole moments a^. Tliese coefficients tell us about the 
underlying theory for the generation of temperature fluctuations and, depending on 
the case, will be taken as following a Gaussian statistics (in which case predictions 
are fully specified by giving the two-point function for the a™'s) or, more generally, 
a non-Gaussian one (and in this case also the bispectrum will be needed,^ as we will 
show below). In either case, the a^'s are stochastic random variables of position, with 
zero mean and rotationally invariant variance depending only on i: (a™(x)) = and 

(a^Hf)a^^*(f)) = 5,,e,5^,„.Ai , (4-18) 

where brackets (■) stand for an average over the ensemble of possible universes. 

Within standard inflationary models, fluctuations are Gaussian and therefore know- 
ing the Cis is enough for characterising the angular fluctuation spectrum. The C^'s 
however are not trivial to compute, specially on small scales (large i) where they en- 
code crucial information regarding the cosmological parameters. On the other hand, 
on large scales where the Sachs- Wolfe contribution to the anisotropy is likely to domi- 
nate the signal, we have ^{x; 7) = |$(x + ro7), with tq = 2/Ho the horizon distance 
and Hq the Hubble constant, and the multipole coefficients are given by 

Ce oc / dkk'^PJk)jj(kro) , (4.19) 
Jo 

where P<i>{k) is the gravitational potential power-spectrum and je is the £-th order 
spherical Bessel function.^ 

For large scales we can make the approximation P<s,{k) oc fc"""^, where n corre- 
sponds to the primordial index of density fluctuations (e.g., n = 1 is the Zel'dovich, 
scale invariant case), in which case [e.g.. Bond & Efstathiou, 1987; Fabbri, Lucchin & 
Matarrese, 1987] 

^ Q^T{i + n/2-l/2)T{9/2-n/2) 



5 r{i + 5/2-n/2)r{3/2 + n/2)'- 



(4.20) 



^Of course, for a complete description we would need all higher-order correlations (o™^ ■ • ■ • ■ ■)• 
However, in what follows we concentrate in the first of these combinations, namely, the bispectrum. 

^The possible infrared divergence of this expression for £ = has no practical effect on observable 
quantities, since the monopole is always removed. 
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with Q = (Q2)V2 ^J^g 

rms quadrupole, which is simply related to the quantity Qrms-ps 
defined by [Smoot et ai., 1992]: Q = \f^Qrms-PS /Tq. 

Notice that for the quadrupole we get Ci=2 = and this depends on the 

spectral index n (i.e., the normalisation of the fluctuation spectrum depends on the 
form of the spectrum). If we instead take = 1 we get £{£ + 1)Q = 6Q^/5, which is 
independent of the multipole index i. This is often referred to as a 'flat' spectrum, since 
gravitational potential fluctuations (as well as the amplitude of the density contrast at 
horizon crossing) are independent of the scale.|^ 

Roughly speaking the amplitude for Ci is determined by fluctuations on angular 
scales 9 n/i radians.]^ As we mentioned in Chapter |^, a plot of i{i + l)Ce is usually 
referred to as the angular power spectrum. With this deflnition an exact n = 1 initial 
spectrum, assuming the relevant scales suffered virtually no evolution after becoming 
smaller than the horizon (as is the case for very large scales), produces a flat spectrum 
(a 'plateau' in the i{i + l)Ce vs. i plot). 

Inflationary theories predict tensor metric fluctuations, in addition to (and statis- 
tically independent of) scalar modes. As a consequence, the total multipole moment 
is actually the sum of two orthogonalQ contributions, which are ordinarily noted as 
C^^'' and C^^^^ for tensor and scalar, respectively. Tilted (n < 1) models, as are pre- 
dicted by e.g., power-law inflation, some versions of intermediate inflation, natural 
inflation and a few other models, imply more power (relative to a flat spectrum) on 
large scales and consequently enhance the contribution from long wavelength tensor 
modes (gravitational waves)P^. This is of relevance for the correct normalisation of the 
spectrum on large scales, as it is the actual case after the detection of anisotropies in 
the COB£^DMR maps. 



^''Here it is useful to remember the handy relation between the comoving wavenumber k and the 
index £: £~ (6, 000h-^Mpc)k. 

^^This correspondence applies for small £ (i.e., large scales). For smaller scales, however, it turns 
out that 9 ~ 60°/^ is more appropriate as a rule-of-thumb. 

^^To first order in the perturbations, scalar and tensor modes are uncoupled, and thus evolve 
independently from each other. 

^•^The exception being natural inflation where the level of gravitational wave contribution in negli- 
gible [Freese et al, 1990]. 



Chapter 5 

Primordial non— Gaussian features 



It has been realised that inflationary models also predict small deviations from Gaus- 
sianity. Every sensible inflation model will produce some small but non-negligible 
non-Gaussian effects, both by the self-interaction of the inflaton field, and by the lo- 
cal back-reaction of its self-gravity. It is our interest here to quantify this prediction, 
with the largest possible generality, so that observational results on the CMB angular 
three-point function can be used as a further test on the nature of the primordial 
perturbation process. 

Unfortunately, we will find that single-field inflation models generally imply mean 
values for the skewness which are well below the cosmic rms skewness of a Gaussian 
field, which confirms and generalises earlier results based on a simple toy-model [Falk, 
Rangarajan & Srednicki, 1993; Srednicki, 1993]. 

Falk et al. [1993] first gave a quantitative estimate of the size of non-Gaussian 
effects through a calculation of the three-point CMB correlation function from pertur- 
bations generated in a simple model, where the inflaton has cubic self-interactions. 

Here the problem is considered in a totally general and self-consistent way. We 
use the stochastic approach to inflation [e.g., Starobinskii, 1986; Goncharov, Linde 



& Mukhanov, 1987; see also § p.l.7|| , as a technical tool to self-consistently account 
for all second-order effects in the generation of scalar field fluctuations during infla- 
tion and their giving rise to curvature perturbations. We also properly account for 
the non-linear relation between the inflaton fluctuation and the peculiar gravitational 
potential. Our derivation moreover removes a non-realistic restriction to purely de 
Sitter background made by Falk et al. [1993], which is especially important when 
non-Zel'dovich perturbation spectra are considered. 
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The technical derivation in the frame of stochastic inflation of the connected two- 
and three-point function (or its Fourier counterpart, the bispectrum) , for the inflaton 
field first and the local gravitational potential next, is reported in § |5.1| . From this 
we compute (§p.2D the CMB temperature anisotropy bispectrum and apply it to the 



collapsed three-point function. § |5.3| deals with the zero-lag limit of the three-point 
function: the skewness, for which we provide a universal inflationary expression, which 
we then specialise to some of the most popular inflaton potentials: exponential [Lucchin 
& Matarrese, 1985], quartic and quadratic [Linde, 1983, 1985] as well as a simple 
potential for hybrid inflation [Linde, 1994]. Interestingly enough, this last model allows 
for a spectral index greater than one.Q In § |5.4| we discuss the effects of the cosmic 
variance on the possibility of observing a non-zero mean value for the skewness. As a 
by-product of this analysis we compute the 'dimensionless' rms skewness for a large 
range of the primordial spectral index of density fluctuations n. Finally, we provide a 
brief discussion in 50| 



5.1 Stochastic inflation and the statistics of the 
gravitational potential 

In order to take into account all the effects contributing to a non-vanishing primordial 
three-point correlation function of $, we will perform the computation in two steps. 
In § p.l.l| we compute the three-point function for the inflaton field perturbation 50. 
The most convenient way to perform this calculation is in the frame of the stochastic 
approach to inflation [Starobinskii, 1986; Goncharov et al, 1987], which naturally takes 
into account all the multiplicative effects in the inflaton dynamics that are responsible 
for the non-Gaussian features. Then, in § ^.1.2| we compute the extra-contribution to 
the three-point function of the gravitational potential that arises due to the non-linear 
relation between $ and This effect has been previously noticed by Barrow & Coles 
[1990] and Yi & Vishniac [1993]. 



^Note that inflation is able to produce perturbation spectra with essentially all values of n around 
unity [e.g., MoUerach, Matarrese & Lucchin 1994], as it can be needed to match the COBE-DMR 
data [Smoot et ai, 1992] with observations on smaller scales. 
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5.1.1 The inflaton bispectrum 

To study the dynamics of the inflaton, we will apply the stochastic approach. This is 
based on defining a coarse-grained inflaton field a), obtained by suitable smooth- 
ing of the original quantum field over a scale larger than the Hubble radius size, whose 
dynamics is described by a multiplicative Langevin-type equation (cf. § |2.1.7|) . This is 



obtained by adding to the classical equation of motion a Gaussian noise term whose 
amplitude is fixed by the rms fluctuation of the scalar field at Hubble radius crossing, 

da ~ 87r\/(0) ^ 277 ' > ' ^ > 

where V^(0) is the inflaton potential, primes denote differentiation with respect to and 
mp is the Planck mass. The Hubble parameter here should be consistently calculated 
from the local energy density of the coarse-grained inflaton. The noise term rj has zero 
mean and autocorrelation function [e.g., Mollerach et ai., 1991] 

{r]{x,a)r]{x\a)) = io{qs{a)\x - x'\)5{a- a) . (5.2) 

The use of the time variable a = In (a/a*) in this equation has been motivated by 
Starobinskii [1986], who noticed that a accounts for the possible time dependence of 
the Hubble parameter. This is particularly relevant when general, i.e. non-de Sit- 
ter, inflation is studied. We also defined the coarse-grained domain size through 
the comoving wave-number qs{a) = eH{a)a{a), with e a number smaller than one, 
H{a) = H {(pciid)) , with 4>d{ci) the homogeneous classical solution of the Langevin 
equation (i.e., that obtained with the noise term 'switched' off). Finally the scale 
factor a{a) is obtained by integration of H{a). 

The stochastic dynamics of the coarse-grained field within a single coarse-graining 
domain (i.e., for x = x') can be studied in terms of the Fokker-Planck equation for the 
probability distribution function of (j). In our case, instead, since we are interested also 
in spatial correlations of the field, we will solve directly the Langevin equation above. 
To the aim of computing the three-point function of 0, a second-order perturbative 
expansion around the classical solution is enough. We will require that the potential 
V{(j)) is a smooth function of its argument, which translates into requiring well defined 
values for the steepness of the potential X{a) = X(0ci(a)) = nipV {(pci) /V {(pci) and 
its derivatives [Turner, 1993] throughout the range of relevant scales. Apart from this 
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requirement we keep the analysis general; only at the end we will apply our results 
to some specific infiationary potentials (§ §5.3.1| - |5.3.2D . We first expand V{(f)) around 
(pel, up to second order in a) = (j){x,a) — (pd{c(), = ^{(pd) + V{(pd)S(p + 



'^V"{(f)ci)S(p'^ + ■ ■ ■■ Replacing this into the Langevin equation we obtain 



d6(f){x, a) 
da 



A{a)5(l){x, a) + B[a)5(j)^{x, a) + [Di[a) + D2{o!)5(l){x, a)] r]{x, a) , (5.3) 



where we have used dcpci/da = — mp^'(0c/)/(87rV^(0c/))- In Eq. (pl3|) we defined 
mp 



A 



-X' : B 



in 



^^-X" ■ Di = ^^^^ ■ D2 - ^^^^^ 



167r 



27r 



(5.4) 



Let us now split the field perturbation as 50 = 50i + 5(f)2 of first and second order, 
respectively. We can also define a rescaled variable fj = {H{a)/2n)ri. We then find 



50i(x, a) = X(a) / da'X {a')7](x,a') 
Jo 



6(j)2(x, a) = X(a) / da' 
Jo 



B(a') 



X(«') ' 2mp 

Let us now calculate the connected three-point correlation function of 6(f). The lowest 
order non-vanishing contribution reads 

{S(f){xi,ai)5(f){x2,a2)5(i)ixs,as)) = 

X{a3) da —— (50i(xi,ai)50i(x3,aO)('^0i(^2, "2)^01(^3, a')) + [^i^^2] 
Jo [X [a') 

X{a3) r^^ 



x,a) + 



1 



-501 (x, a')?7(x, a') 



(5.5) 



(5.6) 



2mp JO 
+ 2x2 terms . 



da' [(50i(xi, ai)(50i(x3, a')) ('^0i(^2, "2)^7(^3, + [^i ^^^2]] 



(5.7) 



The term proportional to X{a-i) /2mp in the r.h.s. of this equation can be recast in 
the form 

r'"'Vx-^(a')^^^(^0i(xi,«i)^0i(x3,«')) Jo{qs{a')\x2 - x,\) (5.8) 
Imp Jo [^TT) 

where we defined amin = min[a3, 02]. We need now to compute the 50 autocorrelation 
function. To this aim, recalling that a is the time when the perturbation wavelength 
^ a/q equals the size of the coarse-graining domain, we change the integration variable 
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in Eq. (|575|) ) from a to g = y 8-KV{a)/3 g*e"/ H^nip. Tlie subscript * denotes quantities 
evaluated at tlie time wlien we start to solve the Langevin equation; this is chosen in 
such a way that the patch of the universe where we live is homogeneous on a scale 
slightly above our present horizon [see, e.g., the discussion by MoUerach et ai., 1991]. 
We then find 

{6Mxua^)SMx3,a')) = — dqq'P{q) ^(^)^(^) Jo(g|xi - f3|) (5.9) 

where we defined P{q) = ^q^^H^{q) where a = a{q). Using this we can rewrite Eq. 
(pl) as 

X{q,)X{q2)X{qs) f d^q' f d\ P{q')P{q) 



2mp J {2Tif J (27r)3 X^{q) 

xW{q'- qrmniqs, q2))Wiq; g™„(gi, g'))e*'r-(xi-3)eir-(x2-x3) (5,10) 

where we defined the filter function W{q; qi) = Q{q — q*) — ^{q — qi) (© is the Heaviside 
function). 

A similar analysis can be performed for the term proportional to B{q) in Eq. ( p.?]) . 

We get 

X{q,)X[q,)X[q,) j^^ -jB(<l)X{<l) j j (^^^^^ 

xW{q"'- g™„(g2, q))W{q"- q^^Qu q'))e^'<-^~-^)e^'^"<^--^^^ . (5.11) 

So far we have been working in configuration space. In order to obtain the gravita- 
tional potential at Hubble radius crossing during the Friedmann era, it is convenient 
to Fourier transform the coarse-grained infiaton fluctuation. One has 5(j){x,a{q)) = 
(27r)^^/ d^k 6(j){k)Q{q — A;)e*^'^, where 6(j){k) denotes the Fourier transform of the full 
scalar field at the time a{q). This follows from the fact that at the time a{q) the only 
modes k that contribute to the coarse-grained variable are those which have already 
left the inflationary horizon, namely k < q. We can then obtain the Fourier transform 

— * — * 

5(j){k, a{q)) = 6(j){k)Q{q — k), and, in the limit k ^ q~ (or equivalently q — > k~^, that 
is, when we consider the horizon-crossing time of the given scale), we simply have 
6(j){k,a{k)) = 6(j){k). In other words, at the time a{k) the Fourier transform of the 
coarse-grained variable coincides with that of the full field. Using these results in Eq. 
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( p.7| ) we finally obtain the inflaton bispectrum through 

^X\k2)e{k2 - k^) + X\k^)Q{k^ -k2) , , /■'^i dq' 



X{k,) 

^X{k2)X{k 



3) 



2m p Jkt q 



+ 2 r ^B{q')X{q') 
Jk q 



^{h^k2} + {ki^k^] . (5.12) 

5.1.2 The three— point function of the gravitational potential 

We want now to compute the three-point function of the peculiar gravitational po- 
tential, ($(ro7i)'^'('"o72)'^'('^o73))- The approximate constancy of the gauge-invariant 
quantity C outside the horizon [e.g., Bardeen, Steinhardt & Turner, 1983] allows to 
obtain the gravitational potential during the matter-dominated era, given the value of 
Scj) during inflation. During inflation one has C,{x,a) ~ —S(j){x,a)/{d(f)/da), which is 
usually interpreted as a linear relation between ( and 6(f). However, when calculating 
the three-point function of $ one cannot disregard the second-order effects coming 
from the fluctuations of d(j)/da. 
Recalling that 

= - 87TW - 87 + ^ ("^^'^l ^'-''^ 

one gets 

C(f) = ^^^^^ (^50(f , a) - §^S<p\x, a)^ . (5.14) 

This equation is expressed in configuration space; the Fourier transform of the first 
term inside the brackets is just 6(j){k,a). For the second term we get a convolution 
of the type / d^p5(p{'p)5(t){k — p). For the scales of interest one has (^{k) ~ — 3C(A;)/5. 
Then, adding the two above contributions and evaluating the expression at the horizon 
crossing time we get 



5mpX{a{k)) 
From this equation we calculate 



"^^^^ " J^^IWf^'^^^' «(fc))50(fc-p, a{k)) 



. (5.15) 
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X{ki)X{k2)X{k^) 



X^{k,)X{k2)X{ks) J (27r)3 
+{h^k2} + {ki^h} (5.16) 

Using the fact that, at horizon crossing, {6(j){ki)6(j){k2)) = {2n)^P{ki)6^{ki + k2), we 
can write 

mhMk2)) = i2nfS%h + k2)f\a{k,))P{k,) , (5.17) 
with f{a{k)) = 247r/5mpX(a(A;)). Finally, using Eqs. (|5J^ ) and (§A^ we find 



^ ^ ^ / z47r \ ^ ^ ^ 

(<l>(A;i)<f (A;2)<f (fca)) = i27:f6'\k^ + ^2 + ^3 



X 



'X2(A;2)X2(A;3) 

■x2(A;2)e(A;2 - ^3) + X2(A;3)e(A;3 - k2) 2X'{k,)X{k2)X{k^ 



+2 ^5(g')X(g0| + {^1^4} + {^1-^4} (5.1^ 



2mp X^ki) 

dq' 
IK q' 

We are interested in considering perturbation modes that left the horizon about 60 
e-foldings before the end of the infiationary epoch. In the explicit examples below 
X{k) turns out to be a slowly varying function of k. We therefore approximate the 
steepness X{k) ~ Xqq in what follows. Then, to the lowest non-vanishing order, the 
three-point correlation function is 

"Xeo 2X^0 , 2 f'^eodq 



Jeo 



2 /■'^so cLa 
+ / -^B{a{q))X{a{q)) 



_2mp Xgo Xeo q 
X {2Txf5\k, + k2 + h) [P^{ki)P^{k2) + P^{k2)P^{k,) + P^{k,)P^{k^)l5.19) 

with /eo = 247r/5mpX6o. Using the expression for P{q) we obtain the power-spectrum 
for the peculiar gravitational potential 

P^ik) = 27r'fi,k-'H\aik))/An' ^ ^fl,Hl,k-\k/k^r~\ (5.20) 

up to possible presence of small logarithmic corrections. In this equation we used the 
fact that the variance per In k of the perturbation in the ( variable, = A;3p^/(27r^) oc 
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f^n-i -g gjygj^ horizon crossing by ~ (27r)^^/J^/0^. By using the equation of motion 
for the inflaton field, 2H ~ — 87r0^/mp we then have H^/cf^ oc H'^/H oc oc H"^ 
(with a ocfP and p > 1 during inflation), and so oc /c"^^, as was in fact used in the 
second equality in Eq. ( |5.20| ) above. 



5.2 The CMB angular bispectrum 

Now we can obtain the two- and three-point correlation functions for $ in configuration 



space by inverse Fourier transforming Eqs. ( |5.171 ) and (|5.19| ) . We then get the following 
general expressions for the mean two- and three-point functions of the primordial 
gravitational potential 

9 



($(ro7i)$(ro72)) = E(2^ + " 72)^^ 



(5.21) 



and 



(<l'(ro7i)<l'(ro72)$(ro73)) 



(An) 



,$3 5:(2j + l)(2£+l)C,Q 

j/>o 



X [Pj(7i ■ 73)P£(7i ■ 72) + Pj(72 ■ 71)^^(72 ■ 73) + Pj{l3 ■ li)Pi{73 ■ 72)] ,(5-22) 



where now the rms quadrupole (cf. Eq. ( [4.20| )) for the multipole coefficients is given 
by 



and 



$3 = — 



60 



2 _ Sti'^HIq r(3 - n)r(3/2 + n/2) 
5mlXio r2(2 - n/2)r(9/2 - n/2) 



fi 



X, 



60 



2m I 



X, 



+ 



60 



X, 



60 



^^60 da 
k, q 



(5.23) 



(5.24) 



In order to obtain Eqs. ( ^.21| ) and (|5.22| ), we had to use the definition of C^. The 



one given in §^]4| differs from the results obtained here because of the window function 
W{k), appearing in the inflationary expressions; this can be however neglected if we 
account for the oscillating behaviour of je for large arguments and for the fact that 
ji for small arguments (that helps in cancelling the lower divergence). 

The equations above allow us to recover Eq. (|4.18|) for the angular spectrum, and 
the following expression for the angular bispectrum. 



(5.25) 
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Figure 5.1: Collapsed three-point function (with and without the quadrupole contri- 
bution), normalised to the skewness, as a function of the angular separation a (left 
panel). The different curves refer to three typical values of the primordial spectral 
index: n = 0.8, 1, 1.2. The symbols in the left panel are the same as in the right one, 
which shows a vertical expansion of the same plot. 



Replacing the latter expression into Eq. (|4.15| ) we obtain the general form of the 
mean three-point correlation function. Some simplifications occur for the collapsed 
three-point function, for which we obtain 

3 

——$3 J2 PeAcosa)J^e,e2e3 
x{2i, + 1)(2£2 + 1)(24 + l)[CeA2 + CeA^ + ^3^,] W,, W,, >V,3 , (5.26) 



where the coefficients Ti^i^^.^ = (47r)^^/ dfl^ J dQyPi-^^^ ■ ^')Pe^{'y ■ ^')Pe^{'y ■ 7') may 
be suitably expressed in terms of products of factorials of ii, £2 and £3, using standard 
relations for Clebsch-Gordan coefficients, by noting that J^kim = (0 o*) ["^f- Messiah, 
1976]. The CMB mean skewness (C3(0)) immediately follows from the above equation 
for a = 0. A plot of the angular dependence of the collapsed three-point function 
above, normalised to the skewness, is reported in Figures ^.1| , for typical values of the 
spectral index n. 
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5.3 The CMB skewness 

In what follows we will restrict our analysis to the mean skewness as given by Eq. 
( p.2(j| ), for a = 0. In the numerical calculation below we will assume that the higher 
multipoles are weighted by a 7°5 FWHM beam, resulting in a = 3°2 [e.g., Wright et 
al, 1992]. To normalise our predictions we will consider the rms quadrupole obtained 
by [Seljak & Bertschinger, 1993] through a Maximum-Likelihood analysis, namely 
Qrms-ps = (15.7 ± 2.6) exp[0.46(l — n)]fiK [see also Scaramella & Vittorio, 1993 and 
Smoot et al, 1994] .0 For the models with n < 1, considered below, this expression 
should be multiplied by an extra factor [(3— n) / (14— 12n)]^/^ to account for the effective 
decrease in the estimated value of Qrms-PS due to the contribution of gravitational 
waves [Lucchin, Matarrese & Mollerach, 1992]. For the mean temperature we take the 
FIRAS determination Tq = 2.726 ± O.OIK [Mather et al, 1994]. 

To estimate the amplitude of the non-Gaussian character of the fluctuations we 
will consider the 'dimensionless' skewness Si = (C3(0))/(C2(0))^^^. Alternatively, if 
we want our results to be independent of the normalisation, we may also define the 
ratio ^2 = (C3(0))/(C2(0))^, as suggested by the hierarchical aspect of our expression 
for the skewness. By writing the multipole coefficients as Ce = (i.e., we factor out 
the amphtude as given by the rms quadrupole), we obtain 



Si 



327r2 



60 



^60 



z 



3/2 



n 



(5.27) 



and 



So 



15 
167r 



X, 



60 



-^60 



(5.28) 



with Q = AmpX^Q J^^°{dq/q)B{a{q))X{a{q)) and where we also defined the spectral 
index-dependent geometrical factor 



Ip{n) : 



|E^i/,/3(24 + l)(24+l)(24+l)[Q,Q,+Q,Q3+C,3C,jW,,W,,W,3J-,, 



E^(2/ + l)CeWi 



(5.29) 



^Note that this value of Qrms-PS assumes a multivariate Gaussian distribution function, account- 
ing for both the signal and the noise; while in principle one should repeat the analysis consistently with 
the assumed statistics of the temperature perturbations, the quasi-Gaussian nature of our fluctuation 
field allows to extrapolate this Maximum-Likelihood estimate without sensible corrections. 
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Figure 5.2: Geometrical factors Ip{n) as a function of n. Crosses are for X^i2{n), 
where the quadrupole term is removed. Diamonds correspond to X^j2{n) including the 
quadrupole contribution. Open boxes correspond to X2(n). For the latter case the 
difference in including or removing the quadrupole contribution is O(10~^), which is 
not perceivable in this graph. 



where the exponent p in the denominator takes values 3/2 and 2 for Si and S2, re- 
spectively. We see clearly from Eqs. (|5.27| ) and ( |5.28| ) that ^2 is independent of the 



normalisation, whereas Si it is not. The numerical factors 1p{n) in Eq. ( |5.29|) are 



plotted in Figure p.2| for different values of the primordial index. 

5.3.1 Worked examples: inflationary models 

Let us now specialise our general expressions to some simple inflationary models. 

Exponential potential 

Let us first consider power-law infiation driven by the exponential potential V{(f)) = 
Vq exp (^—Xy/8TTG(f)j , with A < -\/2 [Lucchin & Matarrese, 1985]. In this case the power- 
spectrum is an exact power-law with n = 1 — 2A^/(2 — A^). We note in passing that 
the right spectral dependence of the perturbations can be recovered using the above 
stochastic approach [Mollerach et al., 1991]. For this model we find X = —^/8n\, 
whose constant value implies A = B = 0. We then have 

r(3-n)r(3/2 + n/2) ^ 



3A Hqo 

Si = - l3/2{n) 

4 mp 



r2(2-n/2)r(9/2-n/2) 



1 5 

52 = y A^ Un) . (5.30) 



94 



Chapter 5. PRIMORDIAL NON-GAUSSIAN FEATURES 



The COBE-DMR results constrain the amphtude of Hqq. For the case n = 0.8 we 
have Hao/mp = 1.8 x 10"^ This gives Si = 9.7 x 10"^ and Si = 1.1 x 10"^ without 
and with the quadrupole contribution respectively, while ^2 = 1.3 in both cases. 



Quart ic potential 



Consider now the potential jXcj)'^ for chaotic inflation [Linde, 1983, 1985]. For 
this model X = Amp/cj) and therefore mpB ~ rn^iX" <^ X. Thus we can take 
^ ~ in the coefficients of Eqs. ( 5.28 ) and ( 5.27| ) for the dimensionless skew- 
ness parameters. By using the slow-roll solution for the infiaton we may write 
X ~ Am 



PV60 



1 + (^^J^mp/HQoj ln(fc/fc6o) where the logarithmic correction to the 
scale invariant power-spectrum is small. Inflation ends when the infiaton takes the 



value (pend - Y2/37r mp implying 
n ~ 1. We find 



^60 



60 /tt mp. In this case the spectral index is 



Si = \ — — X3/2(l) 
V 71 mp 



S. = -(^) 1,(1). 



TT 



P60 



(5.31) 



COBE-BMR constrains the value of A to be A ^ 1.4 x lO'^^ We get Si = 5.2 x 10^^ 
and Si = 6.0 x 10^^ without and with the quadrupole contribution respectively, while 
^2 = 0.5. 



Quadratic potential 



Another simple potential for chaotic infiation is [Linde, 1983, 1985]. In this 

case X ~ 2mp(j)^Q [1 + (mp/47r0gQ) \ia{k/kQo)], (pQQ ~ -y/sO/vr mp and n ^ 1. We find 

5. = ^ ^ ; 5. = lif;^)^.(l). (5.32) 

2 Vtt mp An \ 06o / 

with m^/mp ~ 1.1x10"^ from the COBE-DMR normalisation. We get Si = 3.9x 10~^ 
and Si = 4.5 x 10^^ without and with the quadrupole contribution respectively, while 
^2 = 0.4. 



Hybrid infiation model 
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Finally, let us consider a model of hybrid inflation recently proposed by Linde 
[1994]. Inflation happens in this model during the slow-roll evolution of the inflaton 
in the effective potential V{(()) = Vq + \rn'^4>'^-, where Vq — M^/4 is a cosmological- 
constant-like term. At a given time, when the inflaton fleld takes the value (j) — M, 
its coupling with a second scalar fleld triggers a second-order phase transition of the 
latter (whose vacuum energy density is responsible for the cosmological constant term) , 
which makes inflation end. An interesting prediction of this model is that the spectral 
index of density fluctuations for wavelengths which left the horizon while > m^0^/2 
is larger than unity, namely n ~ 1 + 2m^w?p/ {^tiVq) [MoUerach et al, 1994]. For this 
potential we flnd X — mpm'^ip/Vo, implying Q ~ 0. We then obtain 



r(3-n)r(3/2 + n/2) 
r2(2-n/2)r(9/2-n/2) 



1/2 



(5.33) 



15 mImVL/. 

^60 



1^7^{^-^Tir\Mn) (5.34) 



The slow-roll solution for the inflaton gives 06o = Mexp[60 m?m^p/ (SttVo)]. Choosing 
a value n = 1.1, the COBE-DMK results constrain the free parameters to be M ~ 
1.3 X lO^^mp and m ~ lO^^nip. In this case wc find Si = —1.1 x 10^^ and Si = 
—1.3 X 10~^ without and with the quadrupole contribution, while ^2 = —1.6. Note 
that these results suggest some correlation between the sign of the skewness and the 
value of the spectral index n. 



5.3.2 On the correlation between n and the sign of the skew- 
ness 

Let us now explore somewhat more this correlation between the sign of the skewness 
and the value of the spectral index n. We are here mainly concerned with single- 
field infiationary potentials. Some of the most popular models are characterised by 
their simplicity and universality (such as the quadratic and quartic chaotic potentials 
we already treated), by their being exact solutions of the equations of motion for 
the inflaton (like power-law and intermediate infiation), or by their particle physics 
motivation (as natural infiation with an axion-like potential). 

In contrast with these simpler models where one has just one relevant parameter 
more general potentials, with more freedom, were also considered in the literature. One 
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example of this is the polynomial potential [Hodges et al, 1990] which for an adequate 
choice of the parameters was found to lead to broken scale invariant spectra on a wide 
range of scales with interesting consequences for large scale structure. 

One should also worry about initial conditions [Goldwirth & Piran, 1992]. While 
for single-field models the only effect of kinetic terms consists in slightly changing 
the initial value of (leaving invariant the phase space of initial field values leading 
to sufficient inflation), for models with more than one scalar field initial conditions 
can be important, e.g., for double infiation [Starobinskii, 1985; Kofman et al, 1985; 
Silk & Turner, 1987] (with two stages of infiation, each one dominated by a different 
infiaton field) leading to primordial non-Gaussian perturbations on cosmological in- 
teresting scales. Within the latter models, however, the question of how probable it is 
that a certain initial configuration will be realised in our neighbouring universe should 
be addressed. More recently other examples of interesting multiple-field models with 
broken scale invariance have also been considered [see, e.g., Salopek, 1992]. Here all 
scalar fields contribute to the energy density and non-Gaussian features are produced 
when the scalar fields pass over the interfaces of continuity of the potential. Exten- 
sion of the single-field stochastic approach developed in § ^.1| for calculating the CMB 
angular bispectrum generated through Sachs-Wolfe effect from primordial curvature 
perturbations in the inffaton in order to include many scalar fields is therefore needed. 

It is also of interest here to calculate the accurate form of the spectrum of primor- 
dial perturbations (e.g., by finding the value of the spectral index) at the moment the 
scales relevant for our study left the Hubble radius. Our expression Si for the dimen- 
sionless skewness is accurate up to second order in perturbation theory (the first non 
vanishing order) and therefore we should calculate n at least to the same order. We will 
borrow the notation for the slow-roll expansion parameters from [Kolb & Vadas, 1994]. 
These are definedQ as e(0) = {H'{(t))/H{(j))f /An, 7^(0) = H "{(/))/ {An H{(f))) and 
^(0) = 771^3 H {(p) / {Att H' {(j))) . In the slow- roll approximation e and r] are less than one. 
The same is not true in general for ^ and this may cause consistency problems when 
this term is incorrectly neglected (see the discussion in Ref. [ |113|| ). We will see below 
examples of this. 

Let Q| (Qt) be the contribution of the scalar (tensor) perturbation to the variance 



•^It is easy to relate these parameters to the steepness and its derivatives. In fact e = X^/(167r), 
e + mpX' /{8tt) and £ = e + 3mpX' /{Sir) + w?pX" / {At: X). O f course, the value for e coincides 



with that defined in Eq. ( 2.14 ), whereas rj = —e + x-, cf- Eq. ( 2.15 ). 
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of the quadrupole temperature anisotropy. The complete second-order expressions 
for the tensor to scalar ratioQ and for the spectral index are given respectively by 
R = Qt/QI - 14e [1 - 2C{7] - e)] and 



n=l-2e 



+ 4(C + l)e- (5C + 3)r/ + Ce(l + 2(C + l)e-Cr/) , (5.35) 



evaluated at ~ ^go- In this equation C = — 2 + ln2 + 7~ —0.7296 and 7 = 0.577 is 
the Euler-Mascheroni constant [Kolb & Vadas, 1994; see also Stewart & Lyth, 1993]. 

In the rest of this subsection we will treat three different realisations of the infla- 
tionary idea. After a brief introduction to each of these models^ we compute the CMB 
skewness (cf. Eq. ( p.27| )) and the relation to the spectral index. Let us now turn to 
the examples. 

Natural inflation 

To begin with let us consider the natural inflationary scenario. First introduced 
in [Freese, Frieman & Olinto, 1990] this model borrows speculative ideas from axion 
particle physics [Kolb & Turner, 1990]. Here the existence of disparate mass scales 
leads to the explanation of why it is physically attainable to have potentials with a 
height many orders of magnitude below its width [Adams, Freese & Guth, 1991], as 
required for successful inflation where usually self coupling constants are fine-tuned to 
very small values. 

This model considers a Nambu-Goldstone (NG) boson, as arising from a sponta- 
neous symmetry breakdown of a global symmetry at energy scale / ~ mp, playing the 
role of the inflaton. Assuming there is an additional explicit symmetry breaking phase 
at mass scale A ~ mcuT these particles become pseudo NG bosons and a periodic po- 
tential due to instanton effects arises. The simple potential (for temperatures T < A) 
is of the form ^(0) = A^[l + cos(0//)]. 

For us this V{(f)) constitutes an axion-like model with the scales / and A as free 

parameters. It is convenient to split the parameter space into two regions. In the 

f » TUp zone the whole inflationary period happens in the neighbourhood of the 

minimum of the potential, as may be clearly seen from the slow-rolling equation 

^If the relations were derived for density fluctuation amplitudes rather than for the temperature 
anisotropies there would be a 12.5 instead of the 14 in the expression for R below. 

The models presented here were not described before in this thesis. In this sense, the present 



5 



subsection can be regarded also as a complement to the survey we performed in §2.1.6 
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[Steinhardt & Turner, 1984] |cos(0//)| < 247r/(247r + {nip/ff) which is only vio- 
lated near 0end — tt/, and from the small value of the steepness |X| = ^ tan(^) for 
(p smaller than (pend Q- Thus by expanding around the minimum it is easy to see 
the equivalence between this potential and the quadratic one V ~ m^(0 — 7r/)^/2 with 
= A^/ p. We have already studied the latter in § |5.3.1| and we will add nothing else 
here. On the other hand, let us consider the other regime, where / < mp. Reheating- 
temperature considerations place a lower limit on the width / of the potential. For 
typical values of the model parameters involved, a temperature T^eh ^ lO'^GeV is 
attained. GUT baryogenesis via the usual out-of-equilibrium decay of X-bosons ne- 
cessitates instead roughly T^eh ~ lO^^GeV (the mass of the gauge bosons) for successful 
reheating [Steinhardt & Turner, 1984]. Thus the final temperature is not high enough 
to create them from the thermal bath. Baryon-violating decays of the field and its 
products could be an alternative to generate the observed asymmetry if taking place 
at Treh > lOOGeV, the electroweak scale. This yields the constraint / > O.Smp 
[Freese, 1994], implying n > 0.6 (see below). Attractive features of this model include 
the possibility of having a density fluctuation spectrum with extra power on large 
scales. Actually for / < 0.75 mp the spectral index may be accurately expressed as 
n ^ 1 — mp/87r/^.[] This tilt in the spectrum as well as the negligible gravitational 
wave mode contribution to the CMB anisotropy might lead to important implications 
for large scale structure [Adams et al, 1993]. 

Let us take now / ~ 0.446 nip corresponding to ri ~ 0.8. Slow-rolling requirements 
are satisfied provided the accelerated expansion ends by (j)end — 2.78/, very near the 
minimum of the potential. Furthermore, the slow-rolling solution of the field equations 
yields the value of the scalar field 60 e-foldings before the end of inflation sin(06o/2/) — 



method for reconstructing the inflationary potential (appHed in particular to our present model) 
was developed in [Turner, 1993]; see also [Liddlc & Turner, 1994], for a second-order generalisation 
of this reconstruction process. 

^Recently Kolb & Vadas [1994] have given second-order accurate expressions for 1 — n, ut and R 
(ratio of the tensor to scalar contributions to density perturbations). Applied to our present model 
this importantly changes the relation between R and 1 — n. However the second-order correction to 
n (~ 0.2(mp/87r/^)^ ~ 6 x 10"'^ for our choice of /) is not relevant. 
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exp(— 15mp/47r/^) where we approximated (pend — tt/. We findQ (cf. Eq.( |5.27| )) 

X^o - 4mpX^o = imp/ff[2 + sm\(f)eo/2f)]{l - sm\<Peo/2f))-' - 2{mp/ff. (5.36) 

Two years of data by COBE-DMR [Bennett et al, 1994] are not yet enough to sepa- 
rately constrain the amphtude of the quadrupole and the spectral index. A maximum 
likelihood analysis yields Qrms-ps = 17.6 exp[0. 58(1 — n)]fiK. By making use of Eq. 
( ^.23| ) for the rms quadrupole, COBE-DMR results constrain the value of the free 
parameter A ~ 1.41 x 10~^mp. We find Si ^ 3.9 x 10~^, a rather small signal for the 
non-Gaussian amplitude of the fluctuations. 

Intermediate inflation 

We will now study a class of universe models where the scale factor increases at 
a rate intermediate between power-law inflation -as produced by a scalar field with 
exponential potential - and the standard de Sitter inflation. Barrow [1990] shows that 
it is possible to parameterise these solutions by an equation of state with pressure p 
and energy density p related hjp + p = 7p^, with 7 and A constants. The standard 
perfect fluid relation is recovered for A = 1 leading to the a{t) ~ e^*"^* (-ffm/ constant 
during inflation) solution of the dynamical equations when the spatial curvature k = 
and 7 = 0, while a{t) ~ t^/^'^ for < 7 < 2/3. This non-linear equation of state (and 
consequently the two limiting accelerated expansion behaviours) can be derived from a 
scalar field with potential V = Vq exp(— -^370). On the other hand for A > 1 we have 
a{t) oc exp{At^) (intermediate inflation) with A> and < / = 2(1-A)/(1-2A) < 1 
and again in this last case it is possible to mimic the matter source with that produced 
by a scalar field 0, this time with potential 

''w = (^tmf Its) bm * ) 



In Eq. (5-27) an additional contribution to the squared brackets of the form Q 



imp J^,'^° {dq/q)B{a{q))X{a{q)) was found (the scale fc* signals the time when we start to solve the 
Langevin equation, corresponding to a patch of the universe homogeneous on a scale slightly above 
our present horizon). As we saw, Q turns out to be always negligible or zero for the specific models 
so far considered. As we will see below also here this term yields a value much smaller than the 



others explicitly written in (5.27); e.g., within the frame of natural inflation Q can be worked out 
easily yielding \Q\ ^ 1.1 x lO^^mp//)*, much smaller than Xqq — AuipXqq ~ 2(mp//)^ (/ ^ nip, 
see below). 
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Figure 5.3: (a) Dimensionless skewness Si (in units of the quadrupole Q) as a function 
of d = (fP' — /3^/167r, the value of the squared of the field beyond the minimum allowed, 
for different choices of (3. Curves from bottom to top correspond to (3 from 1 to 8 (left 
panel), (b) Spectral index as a function of d. Now curves from top to bottom are those 
corresponding to (3 from 1 to 8 (right panel). 

with (3 = 4(/-i - 1). 

The equations of motion for the field in a A; = FRW universe may be expressed by 
3H^ = 87r(V^(0) + 0^/2)/mp and + SHcj) = —V. Exact solutions for these equations 
with the potential of Eq. (|5.37| ) are of the form [Barrow, 1990; Barrow & Saich, 1990; 
Barrow & Liddle, 1993] 

H{(f)) = Af{Al3/4nf/\(f)/mp)-^/^ ; = {AP/Any/V^^mp . (5.38) 

Solutions are found for all > but only for 0^ > (/?^/167r)mp we get a > (i.e., 
inflation). In addition /5 > 1 is required to ensure that the accelerated expansion 
occurs while the scalar field rolls (not necessarily slowly) down the potential, in the 
region to the right of the maximum (as it is generally the case). 

From the full potential ( ^.37| ) we may compute the value of the dimensionless 
skewness Si (for convenience we will be taking the field normalised in Planck mass 
units from here on) 



Si = 0.17Q 



(3((3 - 4) 4:P\p - 1) - 1927r/32(/? - 6)0^ 



(5.39) 



02(487r02-/52)2 

evaluated at ~ 06o- In Eq. ( p. 39] ) we tookX3/2(^) — 4.5 (from Eq. ( |5.27| )), that is the 



case for the specific examples we discuss below. A plot of this quantity as a function 
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of (i = 0^ — /?^/167r > (the value of the squared of the field beyond the minimum 
allowed) for different values of (3 is given in Figure |5.3| (left panel). Note that both 
positive and negative values of Si are therefore allowed just by modifying the choice 
of f3. Another generic feature is the rapid decrease of the non-Gaussian amplitude 
for increasing values of the field beyond P/y/lGn (i.e., d > 0). Clearly this is because 
for large (p we approach the slow-roll region where the steepness becomes increasingly 
small. 

Similar calculations may be done for the spectral index. The explicit expression of n 
calculated from ( p.35| ) for our potential ( p).37| ) is complicated and not very illuminating. 
Figure |5.3| (right panel) illustrates the variation of n as a function of d (i.e., the scale 
dependence of the spectral index) for different values of the parameter jS. 

These two figures show that for acceptable values of the spectral index very small 
amplitudes for Si are generally predicted. As an example we consider /3 = 1.2. This 
ansatz yields a negative Si being — 0.37 the value of the field that maximises 
|iSi| ~ 2.25Q. We show in Figure ^ the form of the infiaton potential ( p.37| ) for 
this particular /?. We see that for the scales that exit the Hubble radius 60 e-foldings 
before the end of inflation^ the value of the infiaton, 0eo, is located in the steep region 
beyond the maximum of the potential. The value of the spectral index associated with 
this choice of /5 is n ~ 1.29 (a 3% below the first-order result n ~ 1.34). This value in 
excess of unity for the scales under consideration yields a spectrum with less power on 
large scales (compared with a Harrison-Zel'dovich one) making the long wave length 
gravitational wave contribution to the estimated quadrupole subdominant. The slow- 
roll parameters for this scale are e = 0.16, rj = 0.37 and ^ = 1.010. While the first two 
are smaller than one, ^ is not and so cannot be considered an expansion variable on 
the same footing as e and rj. Terms proportional to ^ (and therefore non negligible) in 
Eq. ( ^.35|) are those not included in the second-order analysis done for the first time 
in Ref.[g03]. Taking the C05^DMR normalisation we finally get Si ~ -4.4 x 10"^ 

Let us now consider (3 = 7. Now the potential falls to zero as a power-law much 



more rapidly than in the previous case. If we take n ^ 0.8 we see from Figure 5.3 



^This model provides no natural end to the inflationary expansion. Therefore we could invoke 
arguments similar to those found in the literature [e.g., Barrow & Maeda, 1990] within which suitable 
modifications of the potential or bubble nucleation (e.g., within an extended inflation model) would 
be responsible for the end of inflation. 

^^Although the slow-rolling assumption is not necessary for finding solutions of the field equations, 
it is instructive to focus also on this limit as it represents the attractor solution for large enough time. 
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Figure 5.4: InHaton potential V{(j)) for the parameter choice {3 = 1.2 (f = 0.77). The 
Geld is taken in units ofmp, while the potential is normalised in units of (lO^^mpA^^^y . 



(right panel) that d ~ 6.05 {(f)QQ ~ 2.65) corresponding to the slow-rolhng region of the 
potentiaQ. For the parameters we find the following values: e = 0.13, t] = 0.17 and 
C, = 0.27. From these we see that the first-order result {n ^ 0.82) is 2.5 % above the full 
second-order one. In this case we get Si ~ 0.50Q. Now gravitational waves contribute 
substantially to the detected quadrupole. Actually we have ^ 1.99 [Liddle & 

Lyth, 1993; Barrow & Liddle, 1993; Kolb & Vadas, 1994] (while we would have had 
~ 1.89 up to first order). Thus the estimated quadrupole should be multiplied by a 
factor (1 + 1.99)^^/^ to correctly account for the tensor mode contribution. Finally, we 
get Si ~ 7.5 X 10"^ 



Polynomial potential 



We are interested in considering a potential of the form 

V{<P) = A(10^ + + f 02) + Vo (5.40) 

where for convenience is written in Plack mass units and A, a and (3 are dimensionless 
parameters. Translation invariance allows us to omit the linear contribution to V. 
A detailed analysis of a potential of the form ( ^.40|) was done by Hodges et al. [1990]. 



Parameter space diagrams were constructed and regions where non-scale invariance 
was expected were isolated. Here we will just summarise what is necessary for our 
study. 



^^In this region we can effectively neglect the term cx ^ in the parenthesis of Eq. (5.37) leading 
to a simplified form for V. 
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Taking f3 > 8a'^/9 ensures that = is the global minimum and therefore Vq = 0. 
If we further require (3 > a"^ then no false vacua are present. Scalar curvature perturba- 
tions are conveniently expressed in terms of the gauge-invariant variable ( [Bardeen, 
Steinhardt & Turner, 1983; Liddle & Lyth, 1993]. The power spectrum associated 
with it, assuming slow-roll evolution of the scalar field in the relevant region of the 
potential, is given by V^^'^{k) oc evaluated at horizon crossing time. Equivalently 

V^^'^{k) oc V^^"^ / {m^V) which suggests that regions of the parameter space where the 
slope of the potential goes through a minimum or a maximum will be of interest as 
far as broken scale invariance is concerned. The presence of these extrema in V is 
guaranteed by taking a"^ < (3 < 4q;^/3, with a < in order for the scalar field to roll 
down the potential from the right. We will thus concentrate our study in the vicin- 
ity of an inflection point 0j, i.e., near the curve /5 = (but still (3 > a"^). In this 
limit we have 0/ — —a/2 and the number of e-foldings taking place in the region of 
approximately constant slope about 0/ is given by iV = —7ia^/^3{j3 — a^) [Hodges et 
al, 1990]. Fixing = 60 as the number required for achieving sufficient inflation, the 
interesting parameters ought to lie close to the curve /3 ~ -|- vr^a^/ (3A^^) ~ -|- ya^ 
with y ~ 9.1 X 10~^. 

We flnd X|o — 4Xgo = 192 y {9y — a~'^)/[a'^{6y + a^^)^] where, again, we are 
normalising the fleld in Planck mass units, y ^ 9.1 x 10~^ and we have evaluated the 
fleld for 060 = 0/ — —a/2. 

Variation of a in the allowed range results in both positive and negative values for 
Si. Clearly /5 > 10a V9 yields >Si > 0. Ahhough in this region of the parameter space 
the value of a that makes Si maximal corresponds to a = —5.44, this value conflicts 
with the requirement (3 < 4a^/3 for the existence of an inflection point. We take 
instead a = —3.69 which makes Si the largest possible one and at the same time agrees 
within a few percent with approximating 0j ~ —a/2. Then /5 ~ + ya^ ~ 15.92. 
This guarantees we are effectively exploring the neighbourhood of the curve f3 = a"^ in 
parameter space. A plot of the potential for these particular values is given in Figure 
|5.5| (left panel). The slow-roll parameters in this case have values: e = 5.94 x 10~^, 
7] = 5.84 X 10-3 and ^ = 1.32 x 10"^ By using Eq. j ^M ) we get n ^ 0.99. 

Let us consider now the case where a"^ < (3 < 10a;^/9. This choice of potential 
parameters leads to Si < 0. The non-Gaussian signal gets maximised for a = —2.24 
and (3 = 5.13. Figure |5.5| (right panel) shows the potential in this case. Note the 
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Figure 5.5: Polynomial potential V{(f)) as a function of the scalar field ((j) is taken in 
Planck mass units). The overall normalisation parameter A is taken to be one; (a) for 
a = -3.69 and /3 = 15.92 (left panel); (b) for a = -2.24 and (3 = 5.13 (right panel). 



resemblance between this form of the potential and that of the hybrid [Linde, 1994] 
model Vo + m^0^/2 for Vq dominating and in the vicinity of the flat region. In that case 
the same sign of the dimensionless skewness (cf. § 5.3.1 ) and 'blue' spectra [MoUerach 
et al, 1994], n > 1, were predicted. For this value of a we get: e = 9.33 x 10~^, 
7] = 6.76 X 10-3 and ^ = 2.76. We get n ~ 1.01. 

We can now use Eq. ( p.23| ) for the quadrupole to find the overall normalisation 
constant A of the potential. Bennett et al. [1994] get a best fit Qrms-ps = (17.6 ± 
1.5)fiK for n fixed to one, as it is our present case to a very good approximation. Thus, 
for a = -3.69 we get A = 2.87 x 10"^^ ^nd 5i ~ 1.1 x 10"^ while for a = -2.24 we 
get A = 5.88 X 10^^^ ^nd Si ~ -1.9 x 10^^. 

The above two values for n show that the departure from scale invariance is actually 
very small (in fact, this is because we are exploring a very narrow range of scales). Note 
also the relatively large value of ^ in the last case compared with e and r]. This tells 
us that the terms proportional to ^ are non negligible in general. Also the rather 
small amplitudes for Si agree with previous numerical analyses [Hodges et al, 1990] 
in which adherence to the correct level of anisotropics in the CMB radiation under 
the simplest assumptions of inflation, like slow-rolling down with potential ( ^.40|) , 
practically precludes any observable non-Gaussian signal. 
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Figure 5.6: Normalised rms skewness of a Gaussian temperature fluctuation field as a 
function of the spectral index n, both including and removing the quadrupole contri- 
bution. 

5.4 The r.m.s. skewness of a Gaussian field 

It was first realised by Scaramella & Vittorio [1991] that detecting a non-zero three- 
point function or skewness for temperature fluctuations in the sky cannot be directly 
interpreted cLS cL SI gnal for intrinsically non-Gaussian perturbations. In fact, even a 
Gaussian perturbation field has non-zero chance to produce a non-Gaussian sky pat- 
tern. This problem is related to what is presently known as cosmic variance, and is 
particularly relevant for fluctuations on large angular scales, i.e., for low-order multi- 
poles of the temperature fluctuation field. One way to quantify this effect is through 
the rms skewness of a Gaussian field {C^{0))Gauss- ^^^Y 



Scaramella & Vittorio [1991] and, more recently, Srednicki [1993], focused on the most 
popular case of a scale-invariant spectrum, n = 1 and the COBE-DMR. window func- 
tion, corresponding to a Gaussian with dispersion a = 3°2 [e.g., Wright et al, 1992]. 
We will be interested here in the same quantity, but for various values of n. In Figure 



of the spectral index, both including and removing the quadrupole: in both cases this 
ratio is in the range 0.1 — 0.3 for interesting values of n. The values obtained for n = 1, 
both with and without the £ = 2 contribution, are identical to those given by Srednicki 





5]6| we have plotted the normalised rms skewness (Cf (0))^{^„g^/ (C2(0)) 



as a function 
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[1993], who adopted the same smoothing angle (the shghtly different definition of win- 
dow function cannot affect our dimensionless skewness ratios). As a rough criterion, 
we can conclude that, in order to detect primordial non-Gaussian signatures, {C^i^O)) 
must be at least of the same order as (Cf (0))^''f„g^. 



5.5 Discussion 

The results reported in the previous sections seem to preclude any chance of actually 
obtaining observable non-Gaussian signals at least in the frame of inflation, unless one 
resorts to more complicated multiple-field models [e.g., Allen, Grinstein & Wise, 1987; 
Kofman et al, 1990; Salopek, 1992]. However, it should be stressed that all the results 



reported in § |5.2| would apply to any large-scale anisotropy where the temperature 
fluctuation can be obtained by a local perturbative calculation (this is not the case, for 
instance, for secondary anisotropies, which are to be ascribed to an integrated effect - 
see Chapter |]): under these conditions the hierarchical form of Eq. ( 5.22 ), where the 
bispectrum is a sum of products of two power-spectra, holds. 

As we have seen, single-scalar-field inflationary models generally lead to very small 
skewness ratios Si = (C3(0))/(C2(0))^'^^, so that their non-Gaussian features cannot 
be distinguished from the cosmic rms skewness. Both for the exponential potential and 
for the quadratic and quartic ones the skewness can be estimated by taking just the 
first term in the brackets of Eqs. ( p.27| ) and (|5.28| ) . An upper limit on its magnitude is 



then obtained by requiring that they give rise to an accelerated universe expansion (i.e., 
inflation), which provides a constraint on the steepness of the potential, < 247r. 
This corresponds to Si < 10~^ and S2 ^ 20. 

Notwithstanding the non-zero value of the skewness, any actual detection of this 
signal is hardly distinguishable from the cosmic variance noise in which it is embedded: 
the intrinsic limitation induced by our impossibility of making measurements in more 
than one universe. In fact, the overall coefficient of Si is generically much smaller 
than the dimensionless rms skewness calculated from an underlying Gaussian density 
field. Removing the quadrupole contribution to reduce the cosmic rms skewness does 
not change this main conclusion, because the predicted mean skewness would also be 
reduced by a comparable factor. Things get even worse if we take into account the 
'sample variance' due to the galactic cut in the maps [Scott et al, 1994]. Therefore, any 
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possible skewness detection on the COBE scale is most probably due to these statistical 
effects than to any primordial non-Gaussian feature. In this sense the quasi-Gaussian 
inflationary predictions for the CMB anisotropics are in full agreement with the recent 
analysis of the three-point function and the skewness from COBE-DMR data [Hinshaw 
et al, 1994, 1995; Smoot et al, 1994]. 
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Chapter 6 
Integrated effects 



6.1 Introduction 

After the detection by COBE-DMR [e.g., Smoot et al, 1992; Bennett et al, 1994] 
of anisotropies in the cosmic microwave background, several other measurements of 
anisotropies at different angular scales have been announced. This new set of obser- 
vations has provided an extra tool to improve our understanding of the structure for- 
mation process. As we have already mentioned in previous sections, the usual analysis 
in terms of the two-point correlation function of the COBE data fixes the rms am- 
plitude of the primordial perturbations. A comparison with measurements at smaller 
angular scales encodes information about the matter content of the universe and other 
cosmological parameters to which these are sensitive. Furthermore, an analysis of the 
three-point correlation function of the data, as it has recently been performed for the 
COBE^MK first- and two-year data [Hinshaw et ai., 1994, 1995], can provide useful 
clues about the statistical properties of cosmological perturbations. Further analyses 
probing the statistical nature of temperature fluctuations have been performed on the 
COB£^DMR maps by Smoot et ai.[1994], Kogut et a].[1994] and Torres et ai.[1995]. 
As we already pointed out, these studies are of particular relevance as they could help 
to distinguish among inflation and topological defects as the source of the primordial 
fluctuations. 

We saw in previous Chapters that inflationary models predict in general a quasi- 
Gaussian distribution of density perturbations. This was due to the fact that they 
are originated by the quantum fluctuations of a very weakly coupled scalar field. The 
effect of the small non-linearities in the inflaton dynamics (cf. § |5.1j ) and the mean 
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value of the resulting imprints on the three-point function of the gravitational po- 
tential [Falk et al, 1993; Gangui et al, 1994; Gangui, 1994] were several orders of 
magnitude smaller than the typical values expected for a particular realisation of an 
ensemble of Gaussian universes [Scaramella & Vittorio, 1991; Srednicki, 1993]. As the 
anisotropies observed at large angular scales are essentially determined by the fluc- 
tuations in the gravitational potential, the CMB three-point correlation function at 
the scales probed by COBE is also expected to be Gaussian. An analysis of the rms 
skewness of inflationary models giving rise to isocurvature baryon perturbations, by 
Yamamoto and Sasaki [1994], yields values smaller than those obtained for curvature 
perturbation models. On the other hand, topological defects are the typical example of 
non-Gaussian distributed perturbations. However, for many observations, the relevant 
object is not the individual effect of each particular defect, but the superposition of 
many of them, which results in a nearly Gaussian pattern [e.g., Gott et ai., 1990; Coul- 
son et ai., 1994; Gangui & Perivolaropoulos, 1995]. We will discuss this in Chapter 
1^ when we will consider an analytic model for the generation of CMB kurtosis from 
cosmic strings. 

Moreover, it has been shown [Scherrer & Schaefer, 1994] that, for a wide class of 
models leading to non-Gaussian density perturbations, the corresponding temperature 
fluctuation field, as induced by the Sachs- Wolfe effect, is nearly Gaussian thanks to 
the Central Limit Theorem. Thus, a detailed analysis of the predictions for various 
models, taking into account all the relevant effects, is necessary to find out the best 
tools to distinguish among them. 

The analysis of the three-point correlations in the COBE-DVlYi two-year 
anisotropy maps shows evidence for a non-vanishing signal in the data but at a level 
consistent with a superposition of instrumental noise and Gaussian CMB fluctuations 
[Hinshaw et ai., 1995]. As the noise level will diminish rapidly with additional data, 
in the four-year map the sensitivity is expected to be limited by the cosmic variance. 

Even for the case of primordial Gaussian curvature fluctuations, the non-linear 
gravitational evolution gives rise to a non-vanishing three-point correlation function 
of the CMB. It has been argued by Luo & Schramm [1993] that the amplitude of this 
effect is several orders of magnitude larger than that predicted by Falk et ai. [1993] 
for a cubic self-interacting inflaton model, and with a similar angular dependence. A 
more recent estimate of this effect on the skewness of the CMB by Munshi, Souradeep 
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& Starobinskii [1995] finds that the amplitude of the gravitational and inflationary 
non-hnearity contributions is comparable for typical inflation models. We will provide 
in this Chapter a detailed and more general analysis of the possible observational 
consequences of the non-linear gravitational growth of initially Gaussian perturbations 
on the CMB three-point function. 

6.2 Contribution to the CMB three— point correla- 
tion function from the Rees— Sciama effect 

At large angular scales the anisotropics in the CMB are given by [e.g., Martinez- 
Gonzalez, Sanz & Silk, 1990; see also Chapter ^ 



where the first term represents the well-known Sachs-Wolfe effect [Sachs & Wolfe, 
1967], while the second one corresponds to the Rees-Sciama, or integrated Sachs- 
Wolfe effect [Rees & Sciama, 1968]. In the previous formula 7 denotes a direction in 
the sky and 77 is the conformal time, with rjo and rjr the present and recombination 
times, respectively. The contribution of the Rees-Sciama effect to the total anisotropy 
is small in a flat matter-dominated universe because the gravitational potential keeps 
constant in time within the linear regime. It gives however a non-vanishing contri- 
bution when the non-linear evolution of perturbations is taken into account. This 
non-linear contribution will generate a non-vanishing three-point function, even for 
primordial Gaussian perturbations in the energy density. 

Let us start off from the Poisson equation, which in comoving coordinates reads 
V^$(x, ?]) = A7cGpa'^6{'x,r]). We will perform a perturbative analysis and therefore 
write, e.g., 5(x, rj) = 5i(x, rj) + 52(x, rj) + . . .. Let us now outhne how to get the first 
two contributions to the gravitational potential (i.e., $(x, 77) = $i(x) + $2(x, 77)). We 
will normalise the scale factor such as ao = 1 (today) and the distance to the horizon 
is given by rjo = 3to = 2/-ffo, in a flat matter domination era (0(77) = {rj/rjoY). Thus, 
the background energy density is p = {?>Hl/S'nG)a~^ , and the Poisson equation may 
be cast as V^$(x, 77) = (3ifo/2a)5(x, r/). Fourier transforming this we readily get 




(6.1) 



$i(k) 



3i/g<5i(k,r/) 
2P a{ri) 



(6.2) 
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to lowest order. In linear perturbation theory we have 6i{k,r]) oc a{r]) and so we get 
the well-known result that $i is time independent^ 

The second-order quantities may be expressed in terms of first -order ones and so 
we have 

52(x, T]) = pi + ^V5, ■ V<l>i + f (^)2$i_^^.$i.^^ (6.3) 

After Fourier transforming and straightforward algebra we get (upon use is made once 
again of the Poisson equation) the second-order gravitational potential [Peebles, 1980; 
Fry, 1984] 

a(r/) 

The three-point correlation function for points at three arbitrary angular separa- 
tions a, 13 and 7 is given by the average product of temperature fluctuations in all 
possible three directions with those angular separations among them. The general 
expression is given in § [4.3| . In this section, for simplicity, we will restrict ourselves to 
the collapsed case, corresponding to the choice a = (3 and 7 = 0, that is one of the 
cases analysed for the COBEyDMH data by Hinshaw et ai. [1994, 1995] (the other is 
the equilateral one, a = /? = 7). The collapsed three-point correlation function of the 
CMB is given by 

C^{x] «) = y J ^^^wi)^ (72)'5(7i ■ 72 - cosa). (6.5) 

For a = 0, we recover the well-known expression for the skewness, 6*3(0). By expanding 
the temperature fluctuations in spherical harmonics, cf. Eq. ( [4.10] ), we can write the 
collapsed three-point function as in equation ( |4.17 ). Then, by taking the mean value 
we get 

The mean angular bispectrum predicted by a given model can be obtained from^ 

(«r«.?«r*) = / rf%rf%rf%>^r*(7i)>^r*(72)>^r(73) 

(^(71)^(72)^(73)) . (6.7) 



^As we have already mentioned, this conclusion also applies in a radiation dominated regime 
^For a general discussion of the properties of the angular bispectrum see, e.g., Luo [1994]. 
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The contribution from the Sachs- Wolfe term has been computed by Falk et al. [1993], 
Gangui et al. [1994] and Gangui [1994], accounting for primordial non-linearities due 
to inflaton self-interactions. The leading contribution coming from the Rees-Sciama 
effect is obtained from 

AT AT AT 

( — (7i) — (72) — (73)) = 



fVO 

J rjr 



) + (7i ^ 73) + (72 ^ 73)- (6.8) 

x=73{»?o-''?) 



At the same order in perturbation theory there are also non-vanishing contributions 
of the type ^{^i{'yir]o)^i{j2rio)^2{%Vo^Vr))] compared to the non-local Rees-Sciama 
terms, these terms are however suppressed by about two orders of magnitude, because 
of both the ry-dependence of $2 and the different numerical factors.^ Using the analytic 
expression for $2 from Eq. (|6.4|), a straightforward computation leads to 



AT AT AT 

(—(71)^(72)^(73)) 



L f oi^ki f d^k2 Z"^" J iki-7i7?0p«k2-72r;0p-«(ki+k2)-73(??0-r?) 

189 J (27r)3 J (27r)3 "''I'l^ ^ ^ 

[m + kl)+l{k, + k2f - (((ki + k2)-k0^ + ((ki + k2)-k2) 



P^{ki)P^{k2) + (71 ^ 73) + (72 ^ 73), (6.9) 

with Pq>{k) the gravitational potential power spectrum. Using this expression and 
well-known integral relations for spherical harmonics, the three angular integrations 
in Eq. ( |6.7| ) can be performed and we obtain for the bispectrum 



J J J Tlr 



■^One might think that contributions of the kind ~ i{ji'rio)^i{j2'r]o) ^2(^3110, Vr)) should, on the 
contrary, dominate the signal since, as no integration over the photon's path is involved, the three 
spots on the last scattering surface -among which we compute correlations- would be closer and 
therefore 'more correlated' than in the situation in which two spots (whose positions are, say, 71 770 
and 72770) la-y on the surface and the third one (73(770 — ??)) sweeps the trajectory of the photon towards 
us (and therefore would be further away from the other two), resulting in a 'loss of correlation' that 
would eventually suppress the signal. However, this is not the case for we should also remember 
the time dependence in $2 (cf- Eq. (^^)). It is precisely the increase in the scale factor (of order 
~ 1000 for standard scenarios) what counter-balances the above mentioned effect and (together with 
the factor 6 difference in numerical coefficients: ^ ^ |) niakes the term in Eq. ( |6.8| ) the dominant 
one. 
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Y,r{k)P^{h)P^{h) , (6.10) 

where je are spherical Bessel functions of order i. Now it is possible to make the 
integrations in dQ^ and perform the summation over rrii in Eq. ( |6.6D , using relations 
among Clebsh-Gordan coefficients^ to obtain for the collapsed three-point function 

DdTT 6Zr]Q ^^^^^^ 

roo roo pi 

/ ciwi<-^ / dw2W^-^ / dz ^ J ^{wlz) 

Jo Jo Jrir/vo ^1+2 ^1+2 



,2\2 



■J, A^2)J, A^2Z) lOK + wl) + 5 ^"^^ In 

^2+2 V W1W2 



Wi - W2 



Wi + W2 



(6.11: 



One can readily see from the last expression the advantages of our approach: any 
primordial spectrum may be studied, as we have P<s>{k) = A{kriQ)"'~~'^ . To find the 
amplitude A we just use Eq. ( p.20| ) and substitute Eq. ( [5.23[ ) in it, to get 

PM - ^ ("AV" '-^O^ r^(2-n/2)r(9/2-n/2) 

5 [kj Tr^^^-"-' r(3-n)r(3/2 + n/2) ' ^^'^^^ 



where we also used Q = v ^i^Qrms-Ps/To- Now recall that the wave number k^ refers to 
the largest scales crossing the horizon today, and so k^ ~ aoH^. As we are taking Oq = 1 
today and r/o = 3to = 2/ifo for niatter dominated regime, we have fc* = 2/rio. Thus, 
we are normalising the amplitude of the spectrum to the Qrms-ps value determined 
by COBE-DMR through the relation 

. _ 36 2'~-A7rQl^,^ps r^(2 - n/2)r(9/2 - n/2) 

5T2 r(3-n)r(3/2 + n/2) ' ^^'^^^ 

Moreover, the expansion in multipoles allows us to specify any chosen window function 
We as well as to subtract the desired multipole contributions from our expression, 
and thus match the different settings of various experiments (provided that we stay 
far enough from the Doppler peak). The integrals in Eg. ( |6.1l|) can be numerically 
evaluated to obtain the prediction for the collapsed three-point correlation function 
produced through the Rees-Sciama effect, to be measured by the particular experiment. 

''in particular, we use the identity J2m=-e(^ '^l^ ^ (^^ + l)<5fco- 
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6.3 The r.m.s. collapsed three— point function 

We turn now to the computation of the cosmic variance associated to the collapsed 
three-point correlation function. If the fluctuations in the CMB temperature were 
Gaussian, the mean value of the three-point correlation function over the ensemble of 
observers would be zero. However, as we are able to perform observations in just one 
particular sky, this prediction comes with a theoretical error bar, or 'cosmic variance', 
that indicates the typical values expected for particular realisations of the Gaussian 
process. Only non-vanishing values larger than this can be interpreted cls cl si gnal 
of intrinsic non-random phases in the distribution of CMB temperature fluctuations. 
The expected amphtude of the cosmic variance (C|(a)) at a particular angular scale 
can be computed from Eq. (|6.5|) 



A-n J 2-n J An J 2tt 
AT, ,AT^ .AT, ,Ar^ 

( — (7i)^ (72)^(73)^ (74)) . (6.14) 

We assume that the multipole coefficients are Gaussian distributed random variables 
with angular spectrum given by Eqs. (^4.18|) and ([4.20|) as it follows from Pq,{k) oc A;""'^. 
Using standard combinatorial properties, we can compute {Cl{a)) Gauss from Eq. (|6.14|) 



{Cl{a))Gauss = E C,A-AM?^ly^lmi + 1)(2^2 + 1) 

Pe,{cosa){Pe,{cosa) + Pe^{cosa) + Pe^{cosa)){ii£20 Ol^sO)^. (6.15) 

This last expression takes into account only the effect of the cosmic variance on the 
collapsed three-point function; one can easily modify it to allow for the experimental 
noise, by adding its contribution to the angular spectrum. 



6.4 Results for the COBE-BMR experiment 

We can now apply our formalism to the COBE-DMR. measurements. We consider a 
Gaussian window function. We ~ exp{—^i{i + l)cr^), with dispersion of the antenna- 
beam profile cr = 3°. 2 [Wright et al, 1992]. Figure ^]T] shows the expected collapsed 
three-point function arising from the Rees-Sciama effect, for three different values of 



116 



Chapter 6. INTEGRATED EFFECTS 




Figure 6.1: The collapsed three-point function C3(a) = Tq(C3(q;)) (where {C^^a)) 
is given in the text), in fj,K^ units, as predicted by the Rees-Sciama effect, vs. the 
angular scale a. The dot-dashed line curves refer to the n = 0.7 spectral index, the 
solid line ones to n = 1 and the dashed line ones to n = 1.3. The top, middle and 
bottom panels represent the cases Imin = 2, 4, 10 respectively. 

the primordial spectral index n = 0.7, 1, 1.3. Note that all three types of perturba- 
tion spectra can be easily generated in the frame of inflation models [e.g., Mollerach, 
Matarrese & Lucchin, 1994, and references therein]. In all cases the perturbation am- 
plitude has been normalised to two-year COBE-DMR. data, using the i = 9 multipole 
amphtude, = 8 /iK, according to the procedure proposed by Gorski et al. [1994], 
which leads to Qrms-ps = 24.2, 19.5, 15.9 /iK, for n = 0.7, 1, 1.3, respectively. The 
top panel is obtained by subtracting from the map only the dipole contribution. In the 
central panel also the quadrupole and the octopole have been removed. In the bottom 
panel all the multipoles up to £ = 9 have been subtracted. This procedure allows easy 
comparison with the recent analysis of the two-year data by the COBE team [Hinshaw 
et al, 1995]. 

Figure |6.2| contains the corresponding plots for the rms collapsed three-point func- 
tion obtained from our initially Gaussian fluctuation field, (C*3 

So far we have assumed full-sky coverage; the effect of partial coverage, due to the 
cut at Galactic latitude |6| < 20°, increases the cosmic variance and can be approx- 
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Figure 6.2: The rms collapsed three-point function {CKa))^^^^^ = Tj^ {C^^{o'))oauss 
(where {Cl{ct))Gauss ■'^ given in the text), in 10^ IJ,K^ units, as arising from the cosmic 
variance of Gaussian fluctuations. Symbols are as in the previous figure. Note that 
the top, middle and bottom panels, representing the cases Imin — 2, 4, 10 respectively, 
have different vertical scales. 
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imately taken into account by multiplying it by a factor 1.56 [Hinshaw et al, 1994]. 
An analytic estimate [Srednicki, 1993; Scott, Srednicki & White, 1994] gives a some- 
what smaller value, 1.23. Note that the subtraction of low order multipoles leads to 
a strong decrease both in the signal and in the cosmic variance. For all considered 
multipole subtractions and angular separations, the expected signal stays typically 
three orders of magnitude below the cosmic variance, which makes the effect unde- 
tectable. The results are found to be weakly dependent on the considered spectral 
index. For comparison, the amplitude of the collapsed three-point function produced 
by non-linearities in the inflaton dynamics [see Gangui et al., 1994] is about a factor 
of ten smaller than that arising from the Rees-Sciama effect. We can also evaluate the 
skewness parameter ^2 = ((AT/T)^)/((AT/T)^)^, which has the advantage of being 
normalisation (i.e., Qrms-ps) independent for the considered effect (cf. §|5.3|); in the 
case, e.g., of only dipole subtraction and n = 1, ^2 ~ —5. Note that, in the same 
case, the heuristic estimate by Luo & Schramm [1993] leads to ^2 ~ 18, while the 
method by Munshi et al. [1995] yielded 1S2 ~ — 1 (slight underestimate - now revised 
to ^2 ~ —5 (—2) with (without) Ofwhp = 10° smoothing). The corresponding rms 
skewness parameter, arising from the cosmic variance of a scale-invariant Gaussian 
field, is S^"^' ^ 1.3 x 10^ {Qrms-ps/l9-5 /^K)-\ 

One may try to increase the amplitude of the signal-to-cosmic variance ratio, by 
referring to experimental settings, such as FIRS [e.g., Ganga et al, 1993] and Tenerife 
[e.g., Hancock et al., 1994], which are more sensitive to intermediate angular scales. 
However, the little increase of this ratio is not enough to make it detectable; further- 
more, one has to bear in mind that the Tenerife measurement is further affected by a 
large sample variance. 

Finally, let us stress that our calculations, reported in [Mollerach, Gangui, Lucchin 
& Matarrese, 1995], apply to measurements of CMB anisotropics on scales large enough 
to be unaffected by the Doppler peak. The possible extension of this effect to smaller 
scales is discussed by Munshi et al. [1995], who however argue that also on those scales 
the cosmic variance largely overcomes the intrinsic non-Gaussian signal. 



Chapter 7 



On the CMB kurtosis from cosmic 
strings 

7.1 Introduction 

We have seen in Chapter ^ how the strong hnk between cosmology and particle physics 
led to the generation of mainly two classes of theories which attempt to provide phys- 
ically motivated solutions to the problem of the origin of structure in the universe. 
According to the one class of theories, based on inflation [e.g., Linde, 1990], primordial 
inhomogeneities in the energy density arose from zero-point quantum fluctuations of 
a scalar field during an epoch of superluminal expansion of the universe. These fluctu- 
ations may be shown to obey Gaussian statistics to a very high degree and to lead to 
an approximately scale invariant power spectrum of density perturbations. According 
to the second class of theories, those based on topological defects [e.g., Vilenkin & 
Shellard, 1994; Hindmarsh & Kibble, 1994], primordial fluctuations were produced by 
a superposition of seeds made of localised energy density trapped during a symmetry 
breaking phase transition in the early universe. Topological defects with linear geom- 
etry are known as cosmic strings and may be shown to be consistent with standard 
cosmology, unlike their pointlike (monopoles) and planar (domain walls) counterparts 
which require dilution by inflation to avoid overdosing the universe. Cosmic strings 
are predicted to form during a phase transition in the early universe by many but not 
all grand unified theories. 

Cosmic strings can explain the formation of large scale filaments and sheets [Vachas- 
pati, 1986; Stebbins et al, 1987; Perivolaropoulos et al, 1990; Vachaspati & Vilenkin, 
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1991; VoUick, 1992; Hara & Miyoshi, 1993], galaxy formation at epochs 2; ~ 2 — 3 [Bran- 
denberger et aJ., 1987] and galactic magnetic fields [Vachaspati, 1992b]. They can also 
account for peculiar velocities on large scales [Vachaspati, 1992a; Perivolaropoulos Sz 
Vachaspati, 1994], and are consistent with the amplitude, spectral index [Bouchet et 
aJ., 1988; Bennett, Stebbins & Bouchet, 1992; Perivolaropoulos, 1993a; Hindmarsh, 
1994] and the statistics [Gott et ai., 1990; Perivolaropoulos, 1993b; Moessner et ai., 
1994; Coulson et ai., 1994; Luo, 1994; Gangui & Perivolaropoulos, 1995] of the CMB 
anisotropics measured by COBE on angular scales of order 9 ~ 10°. Strings may 
also leave their imprint on the CMB mainly in three different ways. The best studied 
mechanism for producing temperature fiuctuations on the CMB by cosmic strings is 
the Kaiser-Stebbins effect [Kaiser & Stebbins, 1984; Gott, 1985]. According to this 
effect, moving long strings present between the time of last scattering tis and today 
produce (due to their deficit angle [Vilenkin, 1981]) discontinuities in the CMB tem- 
perature between photons reaching the observer through opposite sides of the string. 
Other mechanisms through which cosmic strings may produce CMB fiuctuations are 
based, e.g., on potential fiuctuations on the last scattering surface, by inducing den- 
sity and velocity fiuctuations to the matter surrounding moving strings; and on the 
Doppler effect, by producing peculiar velocities in the ambient matter formed by the 
photon's last scatterers. 

It was recently shown [Perivolaropoulos, 1994] how, by superposing the effects of 
these three mechanisms at all times from ti^ to today, the power spectrum of the total 
temperature perturbation may be obtained. It turns out that (assuming standard 
recombination) both Doppler and potential fluctuations at the last scattering surface 
dominate over post-recombination effects on angular scales below 2°. However this is 
not the case for very large scales (where we will be focusing in the present Chapter) 
and this justifies our neglecting the former two sources of CMB anisotropics. The 
main effect of these neglected perturbations is an increase of the Gaussian character of 
the fluctuations on small angular scales. Thus, in what follows wc will be considering 
exclusively anisotropics generated through the Kaiser-Stebbins effect. 

The main assumptions of the model were explained in [Perivolaropoulos, 1993a]. 
The magnitude of the discontinuity is proportional not only to the deficit angle but 
also to the string velocity Vg and depends on the relative orientation between the unit 
vector along the string s and the unit photon wave-vector 7. It is given by [Stebbins, 



7.2. THE FOUR-POINT TEMPERATURE CORRELATION FUNCTION 121 



1988] 

AT 

— = ±47rG/i7,f , 7 ■ {vs x s), (7.1) 

where 7s is the relativistic Lorentz factor and the sign changes when the string is 
crossed. Also, long strings within each horizon have random velocities, positions and 
orientations. In Eq. ( |7. 1] ) , G is Newton's constant and units are chosen such that 
h = c = 1, implying G = nip'^; /i ~ roQUj^ is the effective mass per unit length of the 
wiggly string. G^ ^ 10"^ is needed for consistency of the string scenario with both a 
physically realisable GUT phase transition and a correct level of CMB anisotropies. 

In what follows we will be mainly concerned with the computation of higher order 
correlation functions for CMB temperature anisotropies. In this Chapter, in addition 
to calculating the excess kurtosis parameter for the string-induced CMB temperature 
perturbations, we also show (with explicit formulae) how to estimate the intrinsic un- 
certainties with which the non-Gaussian signal ought to be confronted, when worked 
out at very large angular scales. Although the final result is somewhat disappointing 
(in the fact that the signal cannot be resolved) we think it is interesting to show ex- 
plicit computations (even under simplifying assumptions within an analytic approach) 
of these string-induced higher order correlation functions. In the following section 
we show how to construct the general q-point function of CMB anisotropies at large 
angular scales produced through the Kaiser-Stebbins effect. Explicit calculations are 
performed for the four-point function and its zero-lag limit, the kurtosis. Next, we 
calculate the (cosmic) variance for the kurtosis assuming Gaussian statistics for arbi- 
trary value of the spectral index and compare it with the string predicted value (§ [7.3|) . 



Finally, in §7.4 we briefly discuss our results. 



7.2 The four— point temperature correlation func- 
tion 

We discretise the time between tu and today by a set of N Hubble time-steps tj 
such that tj+i = 2tj, i.e., the horizon gets doubled in each time-step. For a redshift 
zis ~ 1400 we have ^ log2[(1400)^/^] ~ 16. In the frame of the multiple impulse 
approximation [Vachaspati, 1992a] the effect of the string network on a photon-beam is 
just the linear superposition of the individual effects, taking into account compensation 
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[Traschen et al, 1986; Veeraraghavan & Stebbins, 1990; Magueijo, 1992], that is, only 
those strings within a horizon distance from the beam inflict perturbations to the 
photons. According to the previous description of the Kaiser-Stebbins effect, the total 
CMB temperature shift in the 7 direction may be approximately written as the sum 
of all single impulses on the microwave photons caused by each of the strings, over all 
Hubble time-steps N between the time of last scattering and today [Perivolaropoulos, 
1993a] 

Arp N M 

— (7) = 47rG/i7.t;. E E /?"^"(7), (7.2) 

n=l m=l 

where /3'""(7) gives us information about the velocity f™" and orientation s"^" of the 
mth string at the nth Hubble time-step and may be cast as /5™"(7) = 7 ■ with 
pmn ^ ^mn ^ ^mn ^ (sin 6'™" COS f/)™'", siu 6'™" siu 0™", COS ^^™") a uuit vector whose 
direction varies randomly according to the also random orientations and velocities in 
the string network. In Eq. (|7.2|) , M denotes the mean number of strings per horizon 
scale, obtained from simulations [Allen & Shellard, 1990; Bennett & Bouchet, 1988] to 
be of order M ~ 10. 

Given an arbitrary angle on the sky A6, we will see now how to compute the 
number of perturbations (impulses) common to a pair of photon beams (see Figure 



7?TD . We are focusing here in the stringy-perturbations inflicted on photons after the 
time of last scattering, on their way to us. Recall we have discretised the time such 
that the horizon increases by a factor 2 with each time-step. Thus, at a particular 
step p (with l<p<A^~16)we have tp = 2^t/s- For a spatially fiat, matter 
dominated universe (at = a{t) oc t^^^) and in the small angle limit, we may write 
9t^ 1.5°(z(tp)/1400)~^/^, for the apparent angular size of the horizon at tp. Then, 
the ratio of the apparent horizon at time tp to that at last scattering is given by 



2^/3. (7.3) 



Let us now take an arbitrary angle ai2 on the sky. Consider a pair of photons with 
angular separation ai2 = arccos(7i -72) greater than Ot^, for a particular value of p 
[here 71 (72) is the line-of-sight direction of the first (second) photon]. Both of these 
photons will suffer temperature perturbations due to the presence of ~ M strings per 
Hubble time-step. These perturbations, however, cannot be causally correlated at tp, 
since the size of the causal regions surrounding each of these photons is smaller than 
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Figure 7.1: Two photon-beam paths separated by an angle A6, starting off their trip 
at the time of last scattering (noted tree in the plot) and being detected by an observer 
O at to- The horizon in three Hubble time-steps is shown. The effects of strings during 
the first two time-steps are uncorrelated for the particular angular beam separation 
shown. 

the angle au between them. Thus, existing no overlapping of their horizons, each 
photon will be perturbed by different strings, and therefore the perturbations will be 
uncorrelated. 

Given the angle au it is easy to compute the number of steps Nuncor-r during which 
the perturbations on these photons are uncorrelated: just substitute Ot^ by ai2 in 
Eq. (fOP to get Nuncorr{,oii2) = 3 logs (ai2/6'is) , which is valid for au > Ois ^ 1.5°. 
Obviously, the larger the angular scale ai2, the more steps Nuncor-r will be needed for 
the angular size of the horizon to overtake this scale. If we now define — such 
that ai2 = a'i2 + ^is^ subtract Nuncorr from the maximum number of steps of the 
discretisation of the time along the photon path, and finally multiply the result by the 
scaling number of strings per Hubble time-step M, we get 

Korr{ai2) = M[N - N„] = M[N - 31og2(l + ^)], (7.4) 

^Is 

where N'corr{oii2) is the number of correlated 'kicks' inflicted on a pair of photons on 
a scale q;i2, after this scale has entered the horizon and up to the present epoch. In 
this equation, ai2 > and primes were dropped. In other words: A/'corr(«i2) counts 
the total number of correlated kicks onto the pair of photons after the time at which 
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the angle au between their hnes-of-sight fits in the horizon scale. 

Much in the same way, kicks inflicted on three photon-beams will be uncorrelated 
at time tp/ if any one of the three angles between any two directions (say, ai2, 023, dis) 
is greater than 6t^,, the size of the apparent horizon at time tp'. So, in this case, we 
will be summing (cf. Eq.( |7.2|) ) over those Hubble time-steps n greater than p, where 
p (> p') is the time-step when the condition 9tp = max[ai2, ^23, ctis] is satisfied. The 
same argument could be extended to any number of photon-beams (and therefore for 
the computation of the q-point function of temperature anisotropies)Q. 

Let us now study the correlations in temperature anisotropics and put the above 
considerations on more quantitative grounds, by writing the mean q-point correlation 
function. Using Eq. ( [7.2|) we may express it as 

AT AT N M 

(^(7i)...^(7.)) = e E E (/?r"^-/?r"^) (7.5) 

ni,...,nq=l mi,...,mq=l 

where 71... % are unit vectors denoting directions in the sky and where = 7-,- • 

R'^j'^i ^ with 7j = (sin 6*-,- cos 0j, sin sin 0j, cos6'j). In the previous equation we called 
^ = AirGfi'ysVs- We are here defining the ensemble average in the usual way; namely, 
(^(71)^(72)) represents the average of ^(71)^(72), for fixed directions (71, 72) 
over a large number of realisations of the CMB sky. By assuming ergodicity, this is 
equivalent to performing a spatial average over the whole sky sphere, with fixed angle 
ai2 = arccos (71-72), for a single CMB realisation. Furthermore, we may supplement our 
definition with the property that also (■) = during those time-steps when the impulses 
on different photon-beams are uncorrelated. Thus, (■) accounts also for compensation 
effects, according to which at sufficiently early times no string could have perturbed 
both photon-beams whose causal regions were much smaller than the angle between 
them. For simplicity we may always choose a coordinate system such that 61 = 0, 
01 = (i.e., 7i lies on the z axis) and 02 = 7r/2. The seemingly comphcated sum ( fTS]) 
is in practice fairly simple to calculate because of the large number of terms that vanish 
due to lack of correlation, after the average is taken. The calculation proceeds by first 
splitting the product into all possible sub-products that correspond to 

general expression (suitable to any angular scale and to any source of temperature fluctuations) 
for the three-point correlation function was given in [Gangui, Lucchin, Matarrese & MoUerach, 1994]. 
Although of much more complexity, similar analysis may be done for the four-point function and a 
completely general expression in terms of the angular trispectrum may be found. 
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correlated kicks at each expansion time-step (i.e., fj^^j^^'s with the same pair {mj,nj) 
corresponding to the same string perturbing both beams), and then evaluating the 
ensemble average of each sub-product by integration over all directions of EJ^i^i . To 
illustrate this technique we evaluate the two-, three- and four-point functions below. 

Let us start off with the two-point function. Having a correlated pair of beams 
in 7i and 72 directions from a particular time-step p onwards means simply that 
pmim _ pm2n2 ^n = ni = n2>p.,m = mi = 1712', otherwise the ^"^j^-j's remain 
uncorrelated and there is no contribution to the mean two-point function. Therefore 
we will have 

AT AT ^ *^ . 

^(71)^(72) = E E (71 ■ Rn{i2 ■ nn, (7.6) 

n=pm=l 

where we wrote just the correlated part on an angular scale ai2 = arccos(7i-72). The 
uncorrelated part on this scale will vanish when performing the ensemble average (■), 
as we explained before. Thus the mean two-point function may be written as 

AT AT 

(^(71)^(72)) = eW2)Korr{au), (7.7) 

with J\fcorr{cii2) giveu in ( |7.4| ). Since we may always take 

A = (0,0,1)-^; /32 = (sinai2,0,cosai2) , (7.8) 

with R = = (sin 6^ cos 0, sin 6* sin 0, 0056*), (^(71)^(72)) may be calculated by 

integrating over {6, (p) and dividing by An. The result is 

i^T . , AT,„ ,, ^r,cosai2 .r / n /^, 

( — (71)^(72)) = e^^^fc.rMl2). (7.9) 

which may be shown [Perivolaropoulos, 1993a] to lead to a slightly tilted scale invariant 
spectrum on large angular scales. 

The three-point function may be obtained in a similar way. In our analytic approach 
however, the superimposed kernels of the distribution turn out to be symmetric with 
respect to positive and negative perturbations and therefore no mean value for the 
three-point correlation function arises (in particular also the skewness is zero). Similar 
conclusions were derived from exact lattice solutions to the Nambu equations, for string 
evolution and reconnection according to the Smith & Vilenkin [1987] scheme, adapted 
for evolving the network in an expanding universe [e.g., Coulson et ai., 1994]. Notice 
also, that the effects of cosmic string loops that are ignored in the present study are 
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expected to introduce a small skewness on scales of a few arcseconds or smaller (the 
typical angular scale of a string loop on the last scattering surface is a fraction of an 
arcsecond). On the other hand, the four-point function is easily found and a mean 
value for the excess kurtosis parameter [Gangui, 1995] may be predicted, as we show 
below. 

In the case of the four-point function we could find terms where (for a particular 
time-step) the photon-beams in directions 71, 72 and 73 are all correlated amongst 
themselves, but not with the beam 74, which may be taken sufficiently far apart from 
the other three directions. In such a case we have 

,AT,, ,AT,^ xAT,, ,Ar,^ ,AT,, .AT,, ,AT,, ,AT,^ , 

{—{%)—{%)—{%)—{%)) (—(71)^(72)^(73)) (^(74)) (7.10) 

and (cf. the discussion in the previous paragraph) this contribution vanishes. An- 
other possible configuration we could encounter is the one in which the beams 
are correlated two by two, but no correlation exists between the pairs (for 
one particular time-step). This yields three distinct possible outcomes, e.g., 
(^(71)^(72)) {^{13)^(14)) J ^iid the other two obvious combinations. The last 
possibility is having all four photon-beams fully correlated amongst themselves and 
this yields (^(7i)^(72)^(73)-^(74))^, where the subscript c stands for the con- 
nected part. 

In a way completely analogous to that for the two-point function, the correlated 
part for the combination of four beams gives 

AT AT AT AT ^ ^ ^ 

^(71)^(72)^(73)^(74) = E E {irRnii2-Rmi3-Rni%-Rn (7.11) 

n=pm=l 

Now we are ready to write the full mean four-point function as 



AT AT AT AT 

(—(71)^(72)^(73)^(74)) 

" AT AT AT AT 
(— (7i)^(72))(^(73)^(74))+2terms 



AT AT AT AT 
+(^(71)^(72)^(73)^(74)), 



-eKorr(0)cOs(e) 



,AT,. , AT,. , AT,. , AT,. ,, 
+ (— (7i)^(72)^(73)^(74))e (7.12) 



where for simplicity we wrote this equation using the same scale 6 for all directions on 
the sky (first term after the equality sign); after all, we will be interested in the zero-lag 
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limit, in which case we will have 0. By using Mathematica [Wolfram, 1991] it is 
simple to find that the second term includes just the sum of twenty-one combinations 
of trigonometric functions depending on 6'j, 0j, with i = 1,2,3,4, the spherical angles 
for the directions 7j on the sky (cf. Eq. ( |7.11D ). These terms are the only ones which 
contribute non-vanishingly after the integration over the angles 6'™", 0™" is performed 
(assuming ergodicity). 

When we take the zero-lag limit (i.e., aiming for the kurtosis) the above ex- 
pression gets largely simplified. After normalizing by the squared of the variance 
'^^ ~ [|^^-M::orr(0)]^ and Subtracting the disconnected part we get the excess kurtosis 
parameter 

1 AT^ rio 

^ = -l(^ (7i)) -3 ^1.125 



, , ^ , ,,,, _ X 10~ . (7.13) 

Note that for values of the scaling solution M increasing (large number of seeds) the 
actually small non-Gaussian signal /C gets further depressed, as it could be expected 
from the Central Limit Theorem. 



7.3 The r.m.s. excess kurtosis of a Gaussian field 

The previous section was devoted to the computation of /C, the excess kurtosis param- 
eter, as predicted by a simple analytic model in the framework of the cosmic string 
scenario. We might ask whether this particular non-Gaussian signal has any chance of 
actually being unveiled by current anisotropy experiments. Needless to say, this could 
provide a significant probe of the structure of the relic radiation, and furthermore give 
us a hint on the possible sources of primordial perturbations. 

However, as it was realised some time ago [Scaramella & Vittorio, 1991; see also 
e.g., Abbott & Wise, 1984, Srednicki, 1993], the mere detection of a non-zero higher 
order correlation function (e.g., the four-point function) or its zero-lag limit (e.g., 
the kurtosis) cannot be directly interpreted as a signal for intrinsically non-Gaussian 
perturbations. In order to tell whether or not a particular measured value for the 
kurtosis constitutes a significant evidence of a departure from Gaussian statistics, we 
need to know the amplitude of the non-Gaussian pattern produced by a Gaussian 
perturbation field. Namely, we need to know the rms excess kurtosis of a Gaussian 
field. 
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Let us begin with some basics. Let us denote the kurtosis K = J ^^^"^(7) 
and assume that the underlying statistics is Gaussian, namely, that the multipole 
coefficients are Gaussian distributed. Thus, the ensemble average for the kurtosis 
is given by the well-known formula: (K) = 3o"^, where cx^ is the mean two-point 
function at zero-lag, i.e., the CMB variance as given by 



a' ^ {cm = ^ E(2^ + 1) ' (7-14) 



where the Ce coefficients are also dependent on the value for the primordial spectral 
index of density fluctuations n and are given by the usual expression in terms of 
Gamma functions, cf. Eq. ( 4.201 ). In the previous expression represents the 



window function of the specific experiment. In the particular case of the COBE-DMR 
experimental setup we have, for i > 2, Wi ~ exp — + 1)(3.2°)^ , where 3.2° is the 
dispersion of the antenna-beam profile. 

However, (K) is just the mean value of the distribution and therefore we cannot 
know its real value but within some error bars. In order to find out how probable it is to 
get this value after a set of experiments (observations) is performed, we need to know 
the variance of the distribution for K. In other words, we ought to know how peaked the 
distribution is around its mean value (K) . The width of the distribution is commonly 
parameterised by what is called the cosmic variance of the kurtosis: cr^y = {K"^) — (K)"^- 
It is precisely this quantity what attaches theoretical error bars to the actual value for 
the kurtosis. Therefore, we may heuristically express the effect of a^y on K as follows: 
K ~ (K) ± acv, at one sigma level (a good approximation in the case of a narrow 
peak). Rearranging factors we may write this expression in a way convenient for 
comparing it with /C: ICcv = K/a'^ — 3 ~ idcy/cr"^, where ICcv is the excess kurtosis 
parameter (assuming Gaussian statistics) purely due to the cosmic variance. Not only 
is K,cv in general non-zero, but its magnitude increases with the theoretical uncertainty 
(Jcv- This gives a fundamental threshold that must be overcome by any measurable 
kurtosis parameter. Namely, unless our predicted value for K, exceeds ICcv "we will 
not be able to tell confidently that any measured value of /C is due to string-induced 
non-Gaussianities. 

Let us now compute the cosmic variance of the kurtosis explicitly. We begin by 
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calculating (K^) as follows 



71 



i)^ (72))- 



(7.15) 



By assuming Gaussian statistics for the temperature perturbations ^ (7) we may make 
use of standard combinatoric relations, and get 



AT^ AT'' 

(71)— (72)) 



T 



+ 72( 



T 
AT 2 



9(— (71)) (— (72)) 



T 



AT^ AT 

(71)) (72)) (^(71 



T 

.AT 



AT AT 

y (72))'+24(— (71) — (72))'.(7.16) 



As we know, the ensemble averages are rotationally invariant and therefore (^^(71)) = 
(C2(0)) = cr^ is independent of the direction 71. Plug this last equation into Eq. (|7.15|) 
and we get 



{K^) = 9(7^ + 36(7^ l^^d cos a {C2{a)f + 12 J^^dcosa {C2{a)f. 



(7.17) 



The above integrals may be solved numerically. Then, using this result in the expression 
for a^y we get JCcv, the value for the excess kurtosis parameter of a Gaussian field. 



It is also instructive to look at Eq. ( |7.17| ) in some more detail, so that the actual 
dependence on the spectral index becomes clear. By expanding the mean two-point 
correlation functions within this expression in terms of spherical harmonics and after 
some long but otherwise straightforward algebra we find 



/C, 



cv 



72 



+ 24 



,7712 1^3 +m4, njm^ ,1714 , —7773 —7(14 
h, L /T-^3, £4, L 



[E.(2£ + 1)C,W|]4 



1/2 



(7.18) 



where the coefficients 'hC^'^^l!^^ 



! dn^Y^^{^)Yj;^^{^)Y^-'{^) are those defined pre- 
viously in Eq. ( |4.16| ). In the above equation the n-dependence is hidden inside the 
multipole coefficients Cn (cf. Eq.( ^.20| )). 

Eq. ( |7.18D shows an analytic expression for computing Kcv which, in turn, rep- 
resents a fundamental threshold for any given non-Gaussian signal. For interesting 
values of the spectral index (say, between 0.8 < n < 1.3) we find no important vari- 
ation in ICcvi being its value consistent with the Monte-Carlo simulations performed 
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Figure 7.2: Excess Kurtosis parameter of a Gaussian temperature fluctuation field as 
function of the spectral index; imax = 15 and the quadrupole contribution is subtracted. 



by Scaramella & Vittorio [1991]. These authors concentrated on a Harrison-Zel'dovich 
spectrum, considered also the quadrupole contribution and took a slightly different 
dispersion width for the window function (3.0° in their simulations). 

We solved Eq. ( [7. 181 ) numerically. We subtracted the quadrupole contribution 
and performed the sums up to a £max = 15; the exponential suppression of the Wi 
window functions makes higher i contribution to the sums negligible. Our analytic 
analysis expresses JCcv as a ratio of averages, rather than the average of a ratio as in 
[Scaramella & Vittorio, 1991], and therefore exact agreement between their and our 
result should not be expected. Nevertheless, we get the same order of magnitude for 
spectral index n = 1. We also checked that a different n does not change very much 
the value of JCcv- We explored the cosmologically interesting range 0.8 < n < 1.3. 
In Figure |7]2| we plot JCcv versus the spectral index n, for i^nax = 15 and quadrupole 
contribution subtracted; the inclusion of the quadrupole in our results would do no 
more than slightly rise the whole set of points. This is not a new feature and simply 
reflects the intrinsic theoretical uncertainty of the lowest order multipoles. Note the 
small rate of variation of )Ccv with n and that, as expected, JCcv takes larger values 
for smaller spectral indexes. This is clearly due to the fact that a small n generates 
more power on large scales (i.e., small i), and precisely these scales are the ones that 
contribute the most to the cosmic variance of the kurtosis field. 
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7.4 Discussion 

In the foregoing sections we showed how to implement the multiple impulse approxi- 
mation for perturbations on a photon beam (stemming from the effect of the cosmic 
string network) in the actual construction of higher order correlations for the CMB 
anisotropies. We then focused on the four-point function and on the excess kurtosis 
parameter, finding for the latter a value /C ~ 10~^. 

We also found explicit formulae for calculating the rms excess kurtosis K-cv Pre- 
dicted to exist even for a Gaussian underlying field, and showed its dependence on the 
primordial spectral index of density fluctuations. This constitutes the main source of 
theoretical uncertainty at large angular scales, but not the only one: also the 'sample 
variance' [e.g., Scott et al., 1994] makes the signal hard to detect, since we should note 
that no complete sampling of the CMB sky is available, not even the COBE one, where 
galactic latitudes \b\ < 20° are cut out from the maps. These uncertainties completely 
hide the cosmic string signature: this is clearly seen from the results of §7.3| , where 
it was shown that JCcv is two orders of magnitude larger than JC. The results of our 
analytic approach presented in fully agree with those of [Coulson et al, 1994], 
where they found little evidence for non-Gaussianity in their simulated CMB maps, 
for a variety of smoothing angles. We conclude therefore that statistical effects of the 
kind just mentioned would be the most probable responsible of any kurtosis signal 
detected on COBE scales. The results presented in this Chapter have been published 
in [Gangui & Perivolaropoulos, 1995]. 
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Chapter 8 

Analytic modeling of textures 



The analytic modeling of topological defects, like cosmic strings and textures, captures 
the essential features of their cosmological predictions. Analytic models are com- 
plementary to full numerical simulations of defect evolution and CMB temperature 
anisotropy generation, and they rely on these simulations for fixing some free parame- 
ters present in the models. However, large-scale anisotropics alone can only very barely 
discriminate between competing models of structure formation. In the present Chap- 
ter we include work in progress towards studying to what extent CMB non-Gaussian 
features further constrain topological defects models. We derive an analytic expression 
for the collapsed three-point correlation function in the CMB temperature anisotropics 
for textures, as well as for the rms collapsed three-point function (cosmic variance) as- 
sociated with it. Applying our results to the COB£^DMR two-year data, preliminary 
estimates seem to indicate that the CMB non-Gaussianities predicted by this model 
are consistent with the theoretical uncertainties present in the maps, and therefore 
large scale anisotropics alone can not place further constraints on these models. 

8.1 Introduction 

One of the most compelling issues calling for thorough study regards the possibility of 
discriminating the real source of density perturbations leading to large-scale structure 
formation. Competing models are inflationary scenarios and topological defects (cf. 
Chapter |[), these latter formed as the consequence of a symmetry breaking phase 
transition in the early evolution of our universe. 

The CMB radiation carries a valuable information of that early epoch when mat- 
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ter and photons were tightly coupled. A useful way of trying to discriminate between 
models is given by the structure in the relic radiation predicted by them. However, pro- 
duction of anisotropy maps from numerical simulations of the 'defect' field equations 
turns out to be a fabulous problem, and therefore some analytic insight is desirable. 
This is the reason why people turned to consider analytic (although simplified) models 
for anisotropy generation by defects, wherein the physics is more transparent and com- 
putations can be pursued straightforwardly [Perivolaropoulos, 1993; Magueijo, 1995; 
Gangui & Perivolaropoulos, 1995]. Although it is clear that full-range numerical sim- 
ulations will have the last word regarding these and other observable signatures of 
defects, in the meantime progress with these simplified models will hint many of their 
important features. 

Defect induced large-scale structure formation has been well studied recently [see, 
e.g., Vilenkin & Shellard, 1994; Pen et al., 1994] and their predictions confront success- 
fully against a large bulk of observations. The amplitude of the anisotropics is directly 
related to the symmetry breaking energy scale, and normalisation to the COBE-DMR 
detection makes these seeds compatible with a physically realisable GUT phase tran- 
sition. 

However, more and valuable information can be readily extracted from the CMB 
maps. Recently Hinshaw et al. [1995] analysed the two-year COBE-DMR data and 
found evidence for a departure from Gaussian statistics, although consistent with the 
level of cosmic variance. Many are the models wherein non-Gaussian features can 
arise: in non-standard inflationary models, for example, primordial non-Gaussian 
features are easily generated [Falk et al, 1994; Gangui et al, 1994; Gangui, 1994; see 
also Chapter Furthermore, the post recombination mildly non-linear evolution of 
perturbations also predicts a non-Gaussian signature through the integrated Sachs- 
Wolfe effect [Luo & Schramm, 1993; Munshi et al, 1995; MoUerach et al, 1995; see 
also Chapter ^]. Here we will focus on the three-point correlation function of the CMB 
anisotropics. In particular, it turns out that a useful tool to quantify the departure 
from Gaussian statistics is given by the collapsed three-point function (cf. § |4.3|) , and 
we will work exclusively with this latter in the present Chapter. 

One might think that, due to the high level of theoretical noise inherent to the lowest 
order multipoles probed by COBE, very little can be said regarding non-Gaussian 
features in the CMB radiation on large scales. However, based on the analysis of 
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the equilateral and collapsed three-point functions, Hinshaw et al. [1995] were able 
to rule out the primordial isocurvature baryon (PIB) power-law inflationary model 
considered by Yamamoto & Sasaki [1994]: in this model the predicted mean skewness 
vanishes and the rms skewness (cosmic variance) is several orders of magnitude below 
the corresponding level for standard adiabatic models. These authors suggest that the 
cosmic variance of the full three-point function would be tiny as well. This would make 
the theoretical error bars for the three-point function predicted by this PIB model even 
smaller that the actual data from the maps, hence running into conflict with the result 
of the COBE-DMR analysis. 

Likewise, it is interesting to check whether the non-Gaussian signatures predicted 
by some recently proposed analytic defect models are in agreement with the analyses 
of the maps, and if not what the constraints are. This is the main aim of the present 



study, which, we reiterate, is in progress ||70[ . Let us now give a couple of definitions 
to set up the basic framework. 

The three-point correlation function for points at three arbitrary angular separa- 
tions a, P and 7 is given by the average product of temperature fluctuations in all 
possible three directions with those angular separations among them. The general ex- 
pression is given in Eq. ( [4.14| ). Here, for simplicity, we will restrict ourselves to the 
collapsed case, corresponding to the choice a = j3 and 7 = (one of the cases analysed 
for the COBE-DMR data by Hinshaw et al. [1994, 1995]). The collapsed three-point 
correlation function of the CMB is given by 

Csia) = J _2i_(^,)_ (72)^(71 ■ 72 - cos«). (8.1) 

For a = 0, we recover the well-known expression for the skewness, 6*3(0). By expanding 
the temperature fluctuations in spherical harmonics -^(7) = Z)^,™ a-™^^™(7), we can 
write the collapsed three-point function as 

where as usual Wc represents the window function of the particular experiment and 



we use the notation of Eq. ( ^.16| ). 

Predictions for different models usually come as expressions for the angular bis- 
pectrum (a^^a^^a^^). As we saw in Chapter 0, within the analytic model for cosmic 
strings considered there, the bispectrum simply vanishes, and so does the three-point 
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correlation function. In that case we are left just with the amplitude of the rms col- 
lapsed three-point function. For textures the situation is different, and it is one of 
the aims of the present study to show the variation of the collapsed function with the 
angle separation, as well as that of the rms collapsed three-point function. 

8.2 Texture— spot statistics 

One of the most serious complications in predicting CMB fluctuations in defect theories 
is the necessity of evolving the simulations over a large number of expansion times. 
Scaling helps in simplifying this issue in that if the defect is assumed to scale (and 
textures do scale) then the defect network looks statistically the same at any time, 
once its characteristic length is normalised to the horizon. Then, by dividing the sky 
cone into cells corresponding to expansion times and horizon volumes, one can restrict 
the study to the CMB pattern induced in horizon-size boxes during one Hubble time. 
Such studies have been performed by, e.g., Borrill et ai. [1994] and Durrer et al. [1994]; 
in the latter paper, in particular, a 'scaling-spot-throwing' process was implemented 
numerically, with CMB spots derived from a self-similar and spherically symmetric 
model of texture collapses [Turok & Spergel, 1990; Notzold, 1991; see also § |2.2|| . 

The implementation (within the texture scenario) of an analytic approach for the 
computation of the multipole coefficients Ce and the full two-point function, as well 
as for their cosmic variances, was performed in [Magueijo, 1995]. The model allows 
to study texture-induced spots of arbitrary shapes and relies on a hand-full of 'free' 
parameters to be specified by numerical simulations, like the number density of spots, 
n, the scaling size, dg, and the brightness factor of the particular spot, a„, which tells 
us about its temperature relative to the mean sky temperature (in particular, whether 
it is a hot or a cold spot). In the following we will briefiy review the necessary basic 
features of the model and then move on directly to the specific study of interest to 
us, namely the collapsed three-point function predicted by the model, and its cosmic 
variance. We make a short discussion at the end where we summarise what was found 
so far. 

Let us start off discussing the spots statistics. By dividing the whole volume in 
boxes of surface dS and thickness di we may write the volume probability density of 
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spot generation in the sky as follows 

dP = -^dSdidt , (8.3) 
H 

where n is the mean number of CMB spots expected in a horizon-size box volume and 
in a Hubble time. Now, we are interested in a two dimensional distribution of spots 
(since observations map the two dimensional sky sphere), and thus we may integrate 
the thickness d£ in an interval of order to get the surface probability density 

dP - -^dSdt . (8.4) 

Taking into account the expansion of the universe we may express dS = r^{t)dn, with 
r(t) = 3t((to/i)^^^ — 1), where t denotes proper time and to stands for its value today. 

We discretise the time between last scattering, t;^, and today by a set of Hubble 
time-steps t„ such that tn+i = 2tn, i.e., the horizon gets doubled in each time-step. 
Thus, it turns out to be convenient to change time variable to yniin) = log2(to/'^n)- 
For a redshift zis ~ 1400 we have yis — log2[(1400)^/^] ~ 16. In terms of yn the surface 
probability density may be cast as 

dP = N{yr,)dydn, (8.5) 

with 

my.) = -'-^ [2-" - if . (8.6) 

Within the model, CMB anisotropies arise from texture-produced spots in the 
sky. Thus, we may also express the anisotropies in terms of the profiles Wn of these 
spots, generated at time yn, as follows: ^ = I]n '^n^n(^n) Z/n)- Here, a„ characterises 
the brightness of the hot/cold ra-th spot [Magueijo, 1995] and is interpreted as a 
random variable: it is non-Gaussian distributed and its possible mean values are to 
be computed from numerical simulations [e.g., Borrill et al, 1994]. 6n is the angle 
subtended by a scale characterising the size of the spotsQ in the particular Hubble time 
step yn, as measured by a given observer. 

The non-linear evolution of the global scalar field responsible for texture formation 
is the source of the CMB anisotropies. As long as the characteristic length scale of 

^We are here allowing for arbitrary profile shapes. For circular spots the characteristic length is of 
course given by the radius. 
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these defects is larger than the size of the horizon the gradients of the field are frozen 
in and cannot affect the surrounding matter components^. However, when the size 
of the texture becomes comparable to the Hubble radius (this latter growing linearly 
with cosmic time), microphysical processes can 'push' the global fields so that they 
finally acquire trivial field configurations, reduce their energy (which is radiated away 
as Goldstone bosons) and, in those cases where the field was wounded in knots, they 
unwind. In so doing they generate perturbations in the metric of spacetime and hence 
affect also the photon geodesies. These unwinding events are mainly characterised by 
one length scale, that of the horizon at the moment of collapse, implying that the 
outcome perturbations arise at a fixed rate per horizon volume and Hubble time. 

It is this 'scaling' behaviour what makes texture seeds a viable model for structure 
formation; being a generic outcome of any symmetry breaking phase transition involv- 
ing a non-Abelian GUT group, if they did not collapse not only would they be the 
storage of huge amounts of energy (e.g., in the form of gradients of the global fields) 
but they would also produce too much anisotropy in the CMB and would therefore be 
cosmologically intolerable. It is therefore natural to expect a similar scaling for the 
pattern of the spots they produce. A spot appearing at time ?/„ has typically a size 
e'ivn) ^ ds 0^'"'{yn), with e^°''{yn) the angular size of the horizon at ?/„, and where 

Spot anisotropies are generated by causal seeds, and therefore their angular size cannot 
exceed the size of the horizon at that time; thus we have dg < 1. The scaling hypothesis 
implies that the profiles satisfy Wn{9n,yn) = Wn{0n/9'^{yn))- 

Combining both previous expansions for ^ we easily find the expression for the 
multipole coefficients a'^ in terms of the brightness and the profile of the spots 



' 21 
where 



A-jr , . 

7i:^nWiYr{ln) . (8.8) 

-L nr, 



= J ^ Wn{x)Pe{x)dx . (8.9) 



2 



We are here simplifying the discus sion by considering the correlation length ^ < i to be of the 



order of the particle horizon. See Ref. 228 for a complete account of these matters 
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8.3 Spot correlations 

Having derived the coefficients of the expansion of ^ in spherical harmonics, we 
are able now to compute correlations, namely the angular spectrum and angular bis- 
pectrum predicted within this analytic model. From Eq. ( ^.8|) we may easily calculate 
the quantity 

^ E K\' = 7^^^,T.^nan'WX,mn-%,) . (8.10) 

m=—l \ ' nn' 

Now let us compute ensemble averages (■). Any two texture-spots located at dif- 
ferent 7n's are uncorrelated, hence we have 

{a^a^W'^W^,P,{^n ■ In')) = Snn'al {W^y . (8.11) 

In the present model, due to this lack of correlation, the angular part of the surface 
probability density (cf. Eq. ( ^.5|) ) yields just a kronecker delta and disappears from 
the expression of the angular power spectrum = Ci. In fact, from Eqs. ( 8.10 

and ( |8.11| ) we get 



The cumulative effect of the different spots on the Q's (as reflected from the profiles 
W^) is computed from the surface probability density (cf. Eq. (|8.5|) ). Going to the 
continuous limit in the time variable y the above equation may be cast as 

C, = {a^)(-P^yi',' (8.13) 



.2i 

where we defined 

I'jj-'^ ^ I dyN{y)W''{y) . . . W'-{y) . (8.14) 

In these equations and hereafter we drop the „ subindex and consider y in the contin- 
uous limit, (a^) in Eq. ( |8.13| ) is the mean squared value of the spot brightness, to be 
computed from the distribution {an} resulting from simulations. 

8.4 Textures and the three— point function 

Let us now compute the angular bispectrum. Following the same steps as in the 
previous section we may express it as 

{ai, a,^ ) - ^^^^ ^ ^^^^^^ + i)^2i, + 1)^' ^'^'^'^ ' ^^'^^^ 
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This time the angular part of Eq. (^) is the responsible for the appearance of the 
factor Ti^^^g™^. Now we are ready to study the mean collapsed three-point function 
predicted by this model. From Eq. (K^ we get 



{Cs{a))=4n{a^) ^ P,,{cos a)We,We,Wej'/''-' J^e,i,i, , (8.16) 

where Ti^t^io, = (q^ c?) are just the squared of the the 3j-symbols (cf. Messiah 
1976). The factor (a^) in Eqs. ( ^.16|) , ( p.l5|) deserves some comments: a brightness 
distribution {a„} symmetric in hot and cold spots leads to a non-vanishing angular 
spectrum (cf. Eq. ( ^.13|) ). On the other hand, and as we might have expected, such 
symmetric distribution implies (a^) = and therefore yields a vanishing three-point 
function. This is precisely what happens for spots generated by spherically symmetric 
self-similar (SSSS) texture unwinding events [Turok & Spergel, 1990]. 

However, the properties of the SSSS solution are not characteristic of those of 
more realistic random configurations. Borrill et al. [1994] have shown that randomly 
generated texture field configurations produce anisotropy patterns with very different 
properties to the exact SSSS solution. Furthermore, spots generated from concen- 
trations of gradient energy which do not lead to unwinding events can still produce 
anisotropics very similar to those generated by unwindings, and the peak anisotropy 
of the random configurations turns out to be smaller by 20 to 40 % with respect to 
the SSSS solution. Importantly enough, they also find an asymmetry between maxima 
and minima of the peaks for all studied configurations other than the SSSS. In our 
context this implies that (a^) ^ and we will use their results for the amplitude of 
the spots (the brightness) to estimate the amplitude of the mean collapsed three-point 
function below. 

There is no particular form inherent to the spot profiles obtained in the simu- 
lations. The outcome of the random field configurations (especially when no sym- 
metries are assumed) may yield spots of arbitrary distribution of brightness and 
shapes. A Gaussian distribution for the spot brightness, as given by the profile 
W{9/9s{y)) oc exp(— 6'^/6'^(y)) (cf. Eq. (p.7|)) looks as a good approximation to the 
bell-shaped spots obtained by Borrill et al. [1994]. One of the interesting conclusions 
of that work regards the number of spots one is expected to find: whereas in the SSSS 
model approximately only 4 out of 100 simulations showed unwinding events, the num- 
ber density of nonunwinding events (which as we mentioned above also lead to spots) 
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is much greater, typically of the order of one such event per simulation. This implies 
a mean number of CMB spots per horizon-size volume per Hubble time n ~ 1. 

Estimations of the size of the spot indicate that these are approximately 10% of 
the horizon when generated, hence the scaling size ds — 0.1. As we do not have at our 
disposal the maximum and minimum brightness of each of the spots found by Borrill 
et al. [1994] in their simulations, we will estimate (a^) and (a^) from their quoted 
(O'max) — 0.75 X 27re and {amm) — —0.65 x 27re, for the random collapse (and possible 
unwinding) of gradient energy on sub-horizon scales, in order to have a realistic working 
example (e = A-nGrf, cf. §^l2|). We have {a?) ~ 19.44e^ and {a?) ~ 18.26e^. 

We may now apply our formalism to the COBE-DMR measurements and normalise 

1/2 

the perturbations to the two-year data, using the £ = 9 multipole amplitude, Qig = 
8 /iK, according to the procedure proposed by Gorski et al. [1994]. The mean value 
for the CMB collapsed three-point function can hence be computed from Eq. ( p. 161) , 
as well as its dependence on the angular separation [Gangui & Mollerach, 1995]. 

Let us turn now to the computation of the cosmic variance associated to the col- 
lapsed three-point correlation function. We proceed by first calculating the ensemble 
average of the combination of six a™'s, namely, (a^^ . . .aJ^^). The actual expression 
is awkward to read and not very illuminating, and we are therefore not writing it ex- 
plicitly. From that it is then an easy matter to compute the cosmic variance for the 
collapsed three-point function, which we note as a^y^a) = (C|(a;)) — {C^la))'^. We 
get 



2 {Anfiay ^ P,,(a)(P,,(a) + P,,(a) + P,3(a)) T,,,,^,X{X 



\ ' 2£ + 1 



+ (47r)2 (a2)(a^) ^ [P^ (a) + 4P,,(a)P,,(a) + 4Pfe(a)P,,(a) 
£1/2 AAA 



2£i + l 
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+47r (a^) Yl Pi,{a)Pi,{a) Wi,Wi,Wi,Wi,Wf„We, 

xX^^^^'^'*'^'«^.,...3^£4.5^e • (8.17) 

The computation of {C^la)) ± (Jcv{(^) will be an estimate of the maximum band of 
signals predicted from our fluctuation field, and thus allowed by the model. 

8.5 Discussion 

So far we have assumed full-sky coverage in our analytic formulae. We have already 
discussed previously (e.g., in Chapter ^ for the contribution to the three-point function 
from the Rees-Sciama effect), the effect of partial coverage, due to the cut in the maps 
at Galactic latitude |6| < 20°. As estimated in [Hinshaw et al, 1994] this increases the 
cosmic variance by a factor ~ 1.56, and should also be taken into account to quantify 
the uncertainties. 

In Figure ^.1| we show the pseudo-coUapsed three-point correlation function, in 
thermodynamic temperature units, evaluated from the two-year 53+90 GHz average 
COBEI-DMIi map containing power from the quadrupole moment and up [from Hin- 
shaw et al, 1995]. The name 'pseudo-coUapsed' used by Hinshaw et al. [1994, 1995] 
comes from the following fact: whereas in the 'collapsed' case discussed here two out of 
the three directions 7i, 72, 73 on the sky are exactly parallel, in the case when analysing 
the actual data it turns out that taking, say, 71 exactly parallel to 72 is source of large 
errors. This is due to the fact that the temperature associated with each pixel of the 
map comes with a noise contribution. The signal is thus divided into a part with the 
'true' signal plus a bias term that is a cross-correlation of the true sky temperature 
with the average noise pattern. As one can prefigure, given the present noise level of the 
DMR detector and the smallness of the true signal one is after, the bias term dominates 
the behaviour of the statistics. Thus, the strategy of taking 71 'nearly' parallel to 72, 
while leaving 73 arbitrary, is adopted. This technique bypasses the above mentioned 
difficulties and at the same time provides an enhancement in the statistics: instead of 
taking 71 and 72 the same pixel, one takes 71 and the eight nearest neighbours of 72, 
and thus one is summing over roughly four times as many independent pixel triples as 
in the true collapsed case. Needless to say, this enhances the sensitivity [see Hinshaw 
et al, 1994 for details]. 
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Figure 8.1: The 'pseudo-coUapsed' three-point function as computed from the analysis 
of the two-year COBE-DMR data [from Hinshaw et al, 1995]. The quantity plotted 
is the analogue of our {Cs{a))TQ (see the text). Error bars represent instrument noise 
while the grey band represents the rms range of fluctuations due to a superposition of 
instrument noise and cosmic variance. Three different subtraction schemes are shown: 
the top, middle and bottom panels represent the cases Imin = 2, 4, 10 respectively. 



The idea is to confront our results for the mean {C^{a)) plus its cosmic variance 
'error bars' with the experimental data as given in Figure |8.1| . Roughly speaking, if the 
amplitude of our signal plus cosmic variance is larger than the actual data points plus 
the experimental noise given in the maps, then the analytic texture model considered 
in the present Chapter cannot be constrained by the experiment. If, on the other hand, 
our level of cosmic variance (on top of the collapsed three-point signal) lies below these 
points, then (based on the published analysis of the data) something might go wrong 
with the analytic model. Further study may then lead to modifications or even to 
the realisation that some of the assumptions is not fully met, thereby helping in the 
working out of a more realistic version of the model. 
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As we mentioned in the introduction, at the time of the writing up of this thesis, 
the work described here is still in progress. Preliminary (and conservative) estimates 
seem to indicate that the CMB mean three-point correlation function plus the cosmic 
variance predicted within the analytic texture model by Magueijo [1995] is consistent 
with the COBB-DMR two-year data reported in Figure ^]l]. However, no definitive 
statement can be made, since careful checks in our computations are still to be per- 
formed. 

A similar analysis based on a recently proposed analytic model for CMB anisotropy 
generation from cosmic strings [Perivolaropoulos, 1993; Gangui & Perivolaropoulos, 



1995] is also under way |]70|, and we hope to report on this work in the near future. 

The prospect is interesting in itself: it would give us a hint that, even on large 
angular scales, the analysis of non-Gaussian features in the structure of the cosmic 
microwave background radiation maps provides us with a tool to constrain models, or 
at least to make a consistency check of their predictions against experimental data. 



Chapter 9 

Global textures and the Doppler 
peaks 

In the present Chapter we review recent work aimed at showing how global topological 
defects influence the shape of the angular power spectrum of the CMB radiation on 
small scales. As we explained before, the CMB fluctuations on large angular scales 
(greater than a few degrees) are caused by perturbations in the gravitational field 
via the Sachs-Wolfe effect. However, on intermediate scales, 0.1° < ^ < 2°, the 
main role is played by coherent oscillations in the baryon radiation plasma before the 
era of recombination. Unless the universe is reionised at some redshift z > 50, these 
oscillations lead to the 'Doppler peaks' in the angular power spectrum. On even smaller 
scales the anisotropics are damped due to the finite thickness of the recombination 
shell [Bond & Efstathiou, 1984, 1987; Efstathiou, 1990], as well as by photon diffusion 
during recombination, due to the imperfect coupling between baryons and photons 
(Silk damping [Silk, 1968]). 

Inflation-based cold dark matter models predict the location of the first peak to be 
at £ ~ 220/-\/n^, with a height which is a few times the amphtude of the Sachs-Wolfe 
'plateau' on large scales [Steinhardt, 1995]. Here we present a corresponding study for 
perturbations induced by global textures. As we will see below, the first Doppler peak 
is reduced to an amphtude comparable to that of the Sachs-Wolfe contribution, and 
it is shifted to £ ~ 350. We also briefly comment on the relation between these results 
and what open models predict. 
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9.1 Acoustic oscillations and 



AT 
T 



The importance of the observations of anisotropies on angular scales ranging 0.1° < < 
2°, and smaller, cannot be over-emphasised. Presently a host of experiments scrutinise 
different regions of the sky trying to uncover the real amplitude and structure of the 
relic radiation. While the COBE-DMR data played undoubtedly an important role 
in determining the normalisation for the angular spectra predicted by the different 
models, it is clear now that the next important step towards discriminating between 
competing models lies on the degree scale. Whether it will be inflation [Steinhardt, 
1993] or defect models [Kibble, 1980] the 'final' theory, able to explain faithfully the 
experimental data, is something that hopefully will not make us wait so long. The 
CMB anisotropies are a powerful tool to discriminate among these models by purely 
linear analysis. As we saw before, they are parameterised in terms of Cgs, defined (cf. 
Chapter ^) as the coefficients in the expansion of the angular correlation function 

= -^^(2^+ l)QP,(cos^9). (9.1) 

For a scale invariant spectrum of perturbations i{i + l)Ci is constant on large angular 
scales, say i ~ 50. Both inflation and topological defect models lead to approximately 
scale invariant spectra and therefore their shape on these scales does not help in our 
aim of finding a signature to distinguish between models. 

On smaller scales the situation is more promising and we will see below how the 
predictions of global textures on the CMB angular spectrum profile depart from those 
expected in the standard CDM inflationary model. Disregarding Silk damping, gauge 
invariant linear perturbation analysis leads to [Durrer, 1994] 



6T 
~T 



_ i^M + K.y - ^ + $ 



1/ rf 



+ / - ^')dT , (9.2) 

J i 



where $ and \l/ are quantities describing the perturbations in the geometry (Bardeen 
potentials)[| and V denotes the peculiar velocity of the radiation fluid with respect 
to the overall Friedmann expansion. Primes stand for derivatives with respect to 



^It is well known that — $ is the natural variable for the generalisation of the gravitational potential. 
Since in our analysis we neglect anisotropic stresses in the matter components we have ^ = — Hence, 
unlike the notation we used in previous Chapters, we will here stick to the notation in [Kodama & 
Sasaki, 1984] and interpret ^ as the generalised gravitational potential. 
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conformal time. D^J'^ specifies the intrinsic density fluctuation in the radiation fluid. 
There are several gauge invariant variables which describe density fluctuations; they 
all differ on super-horizon scales but coincide inside the horizon. Below we use another 
variable, D^, for the radiation density fluctuation [Kodama & Sasaki, 1984]. Since the 
coherent oscillations giving rise to the Doppler peaks act only on sub-horizon scales, 
the choice of this variable is irrelevant for our calculation. 

$ , \I' and D'^^^ in Eq. ( |9.2D determine the anisotropics on large angular scales^, 
and have been calculated for both inflation and defect models [Scott et al, 1995; 
Bennett & Rhie, 1993; Pen et al, 1994; Durrer & Zhou, 1995]. Generically, a scale 
invariant spectrum is predicted and thus the Sachs-Wolfe calculations yield mainly 
a normalisation for the different models. In this Chapter we present a computation 
for the Doppler contribution from global topological defects; in particular we perform 
our analysis for vra-defects, textures [Turok, 1989], in a universe dominated by cold 
dark matter. Many of the features endowed by textures are also present in other global 
defects and we believe that the analysis we present below remains valid for these defects 
as well. 

The Doppler contribution to the CMB anisotropics is given by 

Doppler ^ 

where x^ec = 5 + 7^0- Iii the previous formula 7 denotes a direction in the sky and r] is 
the conformal time, with 7]^ and rj^ec the present and recombination times, respectively. 

9.2 Linear theory power spectra 

To evaluate Eq. ( |9.3| ) we calculate Dr and V at rjrec- We study a two-fluid system: 
baryons plus radiation, which prior to recombination are tightly coupled, and CDM. 

We start off by considering the evolution equations in each fluid component as 
are given in the classic paper by Kodama & Sasaki [1984]. The evolution of the 

^One might think that o'^g^ leads just to coherent oscillations of the baryon radiation fluid, but 
this is not the case. Note that, e.g., for adiabatic CDM models without source term one can derive 
(l/4)Z)g'-' = —(5/3)^' on super-horizon scales. Since for CDM perturbations $ = —5* and ^P' ~ 0, we 
see how the usual Sachs-Wolfe result 5T/T = (l/3)^'(Xrec, ^rec) is recovered. 
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perturbation variables in a flat background, f2 = 1, can be expressed as 

Dl-3w^Dr = (l + w)[3f^-3$'-A:V;-|(f) k-\l + '^)Vr] ^ 
D', = 3f^-3$'-A;K-|(f)'A;-i(l + ^)K, 

where subscripts r and c denote the baryon-radiation plasma and CDM, respectively; 
Z), V are density and velocity perturbations; w = Pr/Pr, = P'r/p'r P = Pr + Pc- 

Throughout this section we will normalise the scale factor as a{f]eq) = 1, where r]eq 
stands for the conformal time at the time of equality between matter and radiation 
(we will come back to the usual convention a(?7o) = 1 at the end, when expressing 
our results). For a flat universe we have {a' /of = STrGa^^prad + Pmat)/^, where both 
CDM and baryons are included in pmat- Taking into account now the different scaling 
behaviour of both forms of energy with a, and deflning r = (87rG'/3)^^/^(peq/2)^^/^, 
where peg is the total energy density at equality, we get a{r]) + 1 = {\/2 + (77 — r7eg)/2r)^, 
which automatically satisfies it should. By further requiring that a vanishes 

at the origin we get r = rjeq/ {2{\/2 — 1)) while the scale factor may be cast as 

"^^^ = 7(' + ir) • 

It is clear from this expression that for early (late) conformal times we recover the 
radiation (matter) dominated phase of expansion of the universe. 

Furthermore, as both the energy density of baryons and CDM go as oc a^^ we 
have the ratio Pbar/ Pcdm — const, and today yields VLb- From this we may calculate 
Pbar/Prad = {phar / Pcdm){pcdm/ prad) ^ ^sct- Applying this to w and of Eq. (|9]|) we 
get their time dependences as follows 

wc^^-il + nsa)-' ; cl^^-{l + lnBa)-K (9.6) 

The only place where the 'seeds' (texture-seeds or whatever) enter the system ( |9.4| ) 
is through the potentials ^ and $. These potentials may be split into a part coming 
from standard matter and radiation, and a part due to the seeds, e.g., \1' = {c,r) + '^ seed 
where seed-, ^seed are determined by the energy momentum tensor of the seed; the 
global texture in our case. Having said this, one may easily see how the seed source 
terms arise [Durrer, 1994]. 
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It is our aim now to get rid of the velocity perturbations and express the system 
( ^.4| ) as two second order equations for Dr and Dc. For the sake of simphcity in 
these intermediate steps let us neglect, for the time being, the source terms (we will 
incorporate them at the end). After straightforward algebra we find 

^ 12nGa'pr{l + w)Vr-kD'^ 

127rGa^Pr{l + w) + k^ ^ ' ^ 

and 

" k (A;2(127rG'a2)-i + p,(l + w) + p,) ' ^^'^^ 

which express the velocity perturbations in both fluids in terms of the intrinsic density 
perturbations and the background variables. After long (and tedious) computations 
in which we differentiate Vr and Vc from Eqs. ( p.7|) , ( |9.8|) and we substitute back into 
Eqs. ( p. 41) , including this time the seed source term, we finally derive two second order 
equations for and Dc, namely 

D';. + ^[1 + 3c^ -Qw + F-'p,]Dl - %cF-\l + w)D', 

+4'KGa^[pr{3w^ -8w + 6cl - 1) - 2F'^wpc{pr + Pc) 

+p,{9cl-7w) + -——cl]Dr-A7iGa^p,il + w)D, = il + w)S; (9.9) 
47rGa 

D'^ + ^[1 + (1 + w)F-^pr{l + 3cl)]D'^ - ^(1 + 3cl)F-^prDl 

-A7rGa^PcD^-AnGa^Pr{l + 3cl)[l-2{pr + Pc)F-^w]Dr = S, (9.10) 

where F = k'^{127TGa^)^^ + pr{l + w) + pc and S denotes a source term, which in 
general is given by S' = ATiGa?{p + ?)pY^^'^. In our case, where the seed is described by 
a global scalar field 0, we have S = 87rGa^(0)^. 

From numerical simulations one finds that the integral of a^|0p over a shell of 
radius k can be modeled by [Durrer & Zhou, 1995] 

a2(U|2\(A;,r/) = i^I^———, (9.11) 

with fj denoting the symmetry breaking scale of the phase transition leading to texture 
formation (cf. § p.2|) . The parameters in (|9.11|) are A ~ 3.3, a ~ — 0.7/(27r) and 
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(5 ~ 0.7/(27r)^. On super-horizon scales, where the source term is important, this fit 
is accurate to about 10%. Analytical estimates support this finding [^. On small 
scales the accuracy reduces to a factor of 2. By using this fit in the calculation of 
Dj. and Dc from Eqs. (|9.9D, (|9.10|) we effectively neglect the time evolution of phases 



of (0')^; the incoherent evolution of these phases may smear out subsequent Doppler 
peaks [Albrecht et ai., 1995], but will not affect substantially the height of the first 
peak. 

9.3 Choosing initial conditions 

In defect models of large-scale structure formation one usually assumes that the uni- 
verse begins in a hot, homogeneous state. All perturbation variables are zero, and the 
primordial fields (e.g., the global field (f) in our case) are in thermal equilibrium in the 
unbroken symmetry phase. When the phase transition occurs the field (f) leaves thermal 
equilibrium and enters the process of phase ordering. The large scales of interest to us 
get perturbed when crossing the horizon, long after the phase transition took place. 

Before the phase transition the defect fields have no long-range correlations. The 
dynamics of the transition as well as the field ordering are completely causal processes. 
Isotropy and the vanishing correlation function of the defect's stress-energy tensor on 
super-horizon scales lead to the power spectrum of the source term to take the form 
of white noise for kr] « 1. 

Following the phase transition, the distribution of the defect field enters a 'scaling' 
regime rapidly [Pen et ai., 1994], and so numerical simulations need to be performed 
only within a reduced dynamical range: from 'natural' initial conditions (namely, no 
correlations on super-horizon scales) and for a few expansion times. This is enough to 
produce a realistic field distribution on large scales. 

Going back now to our set of Eqs. (|9.9D, ( |9.10|) we will specify our initial conditions 



as follows: for a given scale k we choose the initial time ?7j„ such that the perturbation is 
super-horizon-sized and the universe is radiation dominated. In this limit the evolution 
equations reduce to 
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Figure 9.1: The dimensionless power spectra, k^\Dr\'^ (solid line) and k\D'j.\'^ (dashed 
line) in units of (Ae)^, are shown as functions of k. These are exactly the quantities 
which enter in the expression for the Ce's. We set h = 0.5 , Qb = 0.05 and z^ec = 1100. 



D';-^Dr = (9.12) 



' 3 ^ 

'-D'-^D'-^D = ^ 

c o r r> 9 r /— 

7] ZT] /r]^ y/r] 



D'c + zD'c-izD'r-ihDr = ^, (9.13) 



with particular solutions 



= -^eV^' ; = -Ua^"' . (9.14) 
15 7 

In the above equations we have introduced e = AirGfj'^, the only free parameter in 
the model. We consider perturbations seeded by the texture field, and therefore it 
is incorrect to add a homogeneous growing mode to the above solutions. With these 
initial conditions, Eqs. (|9.9| ), ( p.lOj ) are easily integrated numerically, leading to the 
spectra for Dr{k,rirec) and D[.{k,rirec) (see Figure ^A\j . 
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9.4 Angular power spectrum from global textures 



In this section we calculate the expression for Ci as function of the baryon-radiation 
power spectrum and its derivative, evaluated at r/^ec- 

Following a standard procedure [Durrer, 1994] we express the Fourier transform of 
the velocity perturbation in the radiation fluid as0 



V{k) 



(9.15) 



A;2(l + w) ■ 

Now, let us Fourier transform the Doppler contribution to the anisotropies as given by 
Eq. (U). Using ( ^[5]) we get 

- Doppler 



^(^,^0,/i) 



Jkriofi 



— Dr{k, Tjr 



+ 



ifiD'j.{k, rjr 
k{l + w) 



(9.16) 



where fi = A; ■ 7. Now, it will prove useful to perform an expansion in Legendre 
polynomials as follows 



6T 
~T 



{k,Vo,IJ') 



Doppler 



Y,z'Ae{k,r]o)P,{fi) . 



(9.17) 



After some algebra we get the coefficients in the previous expansion given by 
^eik, Vo) 



4 k[l + w) 



(9.18) 



(2^ + 1) 

where ji denotes the spherical Bessel function of order ^ and j[ stands for its derivative 
with respect to the argument. Now we claim that we may cast the Doppler contribution 
to the multipole coefficients Ci as follows 



p(Doppler) 



dk k^ 



(9.19) 



|A,(fc,r7o)p 
(2£+l)2 ■ 

To prove this we make use of the ergodic hypothesis and substitute the ensemble 
average by an integral over all possible observers. Thus we may write the left-hand 
side of Eq. (Ol) as follows 



AT AT* 

{ — {x,^)— (f,7)) 



T 



T 



(Tx 



dVL^dVL^, AT , ^ ,AT^ 



87r2 



T 



T 



[a;,7 



(9.20) 



•^In the previous section we have calculated the velocity perturbations in both (coupled) fluids. 
From them we found the intrinsic density perturbations. In Eq. (9.15) the velocity perturbation 
refers to that intrinsi c to r adiation at the last scattering surface. Given the uncertainties in the fit of 
the source term (cf. ( |9.11 )) this is a good approximation. See Rcf. [|2j for details. 
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with the constraint 7 ■ 7' = cos'd, and where * denotes complex conjugation. Upon 
using Parseval's relation to transform the / d^x integral into an integral on k we see 
that Eq. ( ^.2(J| ) can be cast as 

1 



(27r) 



(9.21) 



Now, the last integral in Eq. ( |9.21 ) may be written as 



^^^5(7 ■ 7' - cosmCk ■ l)Pe{k . 7'), 



(9.22) 



where we introduced a Dirac delta to enforce the constraint on the 7's. Standard 
algebra like the one used in Chapter ^ yields Eq. ( p. 20 ) in the form 



{C2m 



(9.23) 



(27r)3y {2i+l) 
If we assume now that |A£(A;,?7o)p depends only on the absolute value k, we can 



integrate out the directions k and, upon comparison with the expansion (|9.1| ) it is easy 
to see that the claimed result is achieved. Finally, knowing the spectra for Dr and 
we calculate the Doppler contribution to the Ce's according to 



dk 



16 



\Dr{k,r]rec)\ je{kr]o) 



w)-'\DUk,Vrec)njKkVo)y 



(9.24) 



where we drop the label (-^"pp'^'') from now on. As we mentioned above, the source 
term of equations (|9.9|) and ( |9.10|) was estimated from numerical simulations of the 
evolution of the global scalar field 0. These simulations give a fit for the absolute 
value of and therefore we can say nothing with regard to phases. This also implies 
that we cannot evaluate the crossed term missing in Eq. (|9.24|) , which we neglect as 



a first approximation to the problem. Integrating Eq. (|9.24|) , we obtain the Doppler 
contribution to the angular spectrum of the CMB anisotropics [Durrer, Gangui & 



Sakellariadou, 1995]. We plot the result in Figure ^]2 



9.5 Discussion 



As we see from the plot, the angular power spectrum i{i + l)Ci yields the Doppler 
peaks. For £ < 1000, we find three peaks located at £ = 365, £ = 720 and £ = 950. 
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Figure 9.2: The angular power spectrum for the Doppler contribution to the CMB 
anisotropies is shown in units of e^. We set cosmological parameters h = 0.5 , Qb = 
0.05 and Zrec = HOO. 




Silk damping, which we have not taken into account here, will decrease the relative 
amplitude of the third peak with respect to the second one; however it will not affect 
substantially the height of the first peak. Our most important results regard the 
amplitude and position of the first Doppler peak, for which we find 

i{i + l)C,\,^,,,=5e\ (9.25) 

It is interesting to notice that the position of the first peak is displaced by Ai ~ 150 
towards smaller angular scales than in standard inflationary models [Steinhardt, 1995]. 
This might be due to the difference in the growth of super-horizon perturbations, which 
is Dr oc r]^/^ in our case (cf. Eqs. ( |9.14| )), and Dr oc r^^ for inflationary models. 

One may understand the height of the first peak from the following analytic es- 
timate: matching the sub-horizon with the super-horizon solutions of Eq. ( |9.10| ), in 

the matter dominated era, one finds |Ae(A;/27r)^/2^2_ p^^^ ^q. we then 

obtain in this limit Dr ~ Aek~^^'^. Plugging this latter value into Eq. ( |9.24D we get 
roughly i{i + l)Ce ~ (Ae)"^, for the height of the first peak. 

Let us now compare our value for the Doppler peak with the level of the Sachs- 
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Wolfe plateau [Bennett & Rhie, 1993; Pen et al, 1994; Durrer & Zhou, 1995] 

e{i+l)Ce ,^~(6-14)e^ (9.26) 



It is apparent from Eqs. ( |9.25|) and (|9.26| ) that the Doppler contribution from textures 



is substantially smaller than for inflationary models. 

Normalising with the COBE-DMR experiment [Smoot et al, 1992; Gorski et al, 
1994] 

i{i + 1)C, I ~ 8 X 10-1° , 

one finds the symmetry breaking energy scale when textures are produced in the early 
universe, and from this e ~ 10~^. This value for e depends of course upon the numerical 



simulations, Refs. 123, 161, 54 



We believe that this result, stating that the first Doppler peak has a height compa- 
rable to the Sachs- Wolfe plateau, is valid for all global defects. This depends crucially 
on the 1/y^ behavior of (a0)^ = 0'^ on large scales (cf. Eq. ( |9.11| )), which is a generic 



feature of global defects: on super-horizon scales (p'^{k) represents white noise superim- 
posed on the average given by (0'^) ~ 4'''^{k = 0) oc l/?7^ (the usual scahng behaviour). 
The Fourier transform of 0'^ determines the fluctuations on this background on a given 
comoving scale A. A patch of size A consists of = (A/?])^ independent horizon-sized 
volumes. The fluctuation (f)''^{k) on this scale is thus proportional to (0'^) / oc 1/ y/rj. 
Notice that this argument does not apply for local cosmic strings. 

The displacement of the peaks towards smaller scales is reminiscent of open models 
[e.g., Kamionkowski, Spergel & Sugiyama, 1994], where the first Doppler peak is lo- 
cated approximately at ipeak — 220/ v^Ho. A value Qq ~ 0.3 would produce ipeak — 365 
for the first peak, as we actually find. However, the amplitude of the peaks within 
open inflationary models is still much larger than what we find from global defects. 

Similar results to those worked out in the present Chapter were independently 
derived in [Crittenden & Turok, 1995]. Following a different approach, these authors 
calculate the Doppler peaks from cosmic textures in the synchronous gauge. A main 
assumption of that analysis is that spatial gradients in the scalar field are frozen outside 
the horizon, and therefore time derivatives are negligible. This is not the case in the 
simulations for the texture source term from Durrer & Zhou [1995]. Apart from this 
we basically agree with the shape and position of their Doppler peaks. Regarding the 
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height of the peaks, however, these authors seem to obtain an amphtude similar to 
that predicted in standard inflation-based models, unlike what we find. 

Based on our analysis we conclude that if the existence of Doppler peaks is indeed 
confirmed (recall that not all the teams agree in their findings on the ^ ~ 1° scale) 
and the height of the first peak is larger than about twice the level of the Sachs-Wolfe 
plateau, namely 



and if the first peak is positioned at £ < 300, then global topological defects are 
ruled out. On the other hand, if the first Doppler peak is positioned at £ ~ 350 and 
its height is below the above value, global defects are strongly favored if compared 
to standard cold dark matter inflationary models. This is a clear Gngerprint, within 
present observational capabilities, to distinguish among these two competing models 
of structure formation. 



peak 



> 2 X lO"'^ , 
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